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1. INTRODUCTION 


GENERAL 


FORTRAN (from FORmula TRANslator) is a programming language designed for 
extensive use in mathematical, scientific, and technological areas. The advantages 
of FORTRAN are minimum programming time and cost, and maximum interchange- 
ability of FORTRAN programs on different FORTRAN processors. 


FORTRAN statements resemble English statements and the equations of elementary 
algebra. Therefore, FORTRAN statements are self-documenting, that is, the intended 
operation is apparent from the statement itself. For example, to find the average of 
two numbers, the programmer can write a statement such as: 


AVRGE = (A+B)/2.0 


Since the FORTRAN programmer uses a programming language that resembles the 
language ordinarily used for the solution of problems, relatively littie time is 
required to leam the language. As a result, programming effort can be devoted to 
the logic of the problem without being troubled by the intricacies of computer opera- 
tion. This self-documenting feature of FORTRAN reduces debugging time and- 
enables other programmers to readily grasp the logic of a program so that it can be 
modified or adapted to other purposes with minimal effort. 


Because FORTRAN is the first programming language to be generally accepted as a 
standard by the data processing community, a FORTRAN program written for a partic- 
ular FORTRAN processor can be accepted by many different FORTRAN processors 
with a minimum of change. The FORTRAN specifications described in this manual 
are those of the American National Standards Institute, Inc. (ANSI, formerly known 
as USA Standards Institute, Inc. (USASI), in FORTRAN, X3.9-—1966. Another ANSI 
standard in current use is Basic FORTRAN, X3.10—1966, which is a subset of 
FORTRAN. Differences between the two standards are described in Appendix A. 


Fundamentals of FORTRAN is designed to introduce FORTRAN to the novice 
programmer while providing sufficient depth of coverage for experienced FORTRAN 
programmers. The examples and descriptions of operation, therefore, will be partic- 
ularly valuable for the novice. 
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1.2, PROGRAMMING LANGUAGES 


The series of steps specified for the solution of a particular problem is called a 
source program and the notation that the programmer uses for specifying these steps 
is a programming language. Except in very few cases (where the programmer uses 
machine language) the computer cannot ‘‘understand’’ the programming language; 
the programming language must be translated into instructions that the computer can 
comprehend. Any programming language that resembles English or the language of 
mathematics must be translated into machine language before it can have any effect 
upon computer operation. The combination of the mechanism that accomplishes the 
translation from programming language to machine language and the data processing 
system is called the processor. 


The following paragraphs trace the development of programming languages from 
machine language to assembly language to FORTRAN and show the application of 
each to the problem of evaluating Z in the following equation: 


z -(R+S$~-T)Xx 
<3 


where the values of R, S, T, X, and Y are known. 
1.2.1. Machine Language 


The fundamental unit of information handled by a data processing system is the 
bit (from binary digit). In the strictest sense, machine language is a combination 
of bits interpreted by a data processing system as an instruction. Each bit has 
two mutually exclusive states, represented by 0 and 1. However, almost every 
data processing system has built-in facilities for accepting an abbreviated form 
of machine language. In an ‘‘octal’’ computer, each of the octal digits, 0 through 
7, represents a combination of three bits. For example, an octal 6 is the equiva- 
lent of binary 110. In a ‘‘hexadecimal’’ computer, the digits 0 through 9 and A, 
B, C, D, E, and F each represent four bits, e.g., the hexadecimal D represents 
the binary 1101. In a ‘‘decimal’’ computer, each of the digits, 0 through 9, repre- 
sents four bits. For example, the decimal 91 may be represented internally as 
1001 0001. The following descriptions of machine language refer to this abbreviated 
form, using a typical octal computer. 


Each computer has its own unique set of instructions and machine language. Using 
machine language, the programmer must know the operation code (indicating the 
operation to be performed), the location (in main storage) of at least one operand, 
and the location which is to contain the result. Typical machine language coding 
for the problem described in 1.2 is indicated after the following assumptions: 


The value of R is in location 10002. 
The value of S is in location 10010. 
The value of T is in location 10020. 
The value of X is in location 10030. 
The value of Y is in location 10040. 
The result, Z, is to be stored in location 10050. 
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@ A typical machine language set of coded instructions for performing the evaluation 
is: 


772110002 moves the value of R to an arithmetic register. (A register is a fixed 
location which stores an operand or the result of an arithmetic operation.) 


770110010 adds the value of S to the value in the register and retains the result in 
the register, giving R +S. 


770210020 subtracts the value of T from the value in the register and retains the 
result in the register, giving R +S — T. 


770310030 multiplies the value in the register by the value of X and retains the 
result in the register, giving (R + S — T) X. 


770510040 divides the value in the register by the value of Y and retains the result 
in the register, giving ((R + S — T)X]/ Y. 


. 


772510050 moves the contents of the register to the location specified for Z. 


The first four digits of each instruction in this example constitute the operation 
code; the rightmost five digits specify a memory location. 


Machine coding is tedious; it requires time for learning (since it is unique), the 
programmer must keep track of all memory references, and it requires many lines 

of source coding for a relatively simple operation. However, machine coding requires 
the least machine time since no language conversion is involved. 


1.2.2. Assembly Language 


& The next logical step from machine language is what is commonly known as assembly. 
language. Assembly. language requires a language translation program (an assembler). 
Using assembiy language, the programmer is permitted to use symbolic references to 
memory and specified mnemonic codes instead of numeric codes to designate the 
operation to be performed. Typical assembly coding for the problem in 1.2 is: 


ENT*R enters R into the register. 


ADD*S adds S to R in the register. 


uae 


SUB*T subtracts T from the contents of the A register and retains the result in 
the register. 


MUL*X multiplies the contents of the register by X and retains the result in the 
register. 


DIV*Y divides the contents of the register by Y and retains the result in the register. 


STR*Z_ stores the contents of the register in Z. 


Assembly language coding is less tedious than machine coding, and the operation 
code is more meaningful, making the language easier to learn. The assembler 
assigns memory locations to the symbolic memory references, relieving the program- 
mer of that chore, thus making the coding even more meaningful. However, it is 
essentially a one-to-one source coding process (i.e., one source language instruction 
generates one machine language instruction) when compared to machine coding and 
still requires many lines of source coding for a simple problem. 
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1.2.3. FORTRAN Language 


FORTRAN is one of many higher level languages than have evolved from assembly 
language. It is considered a higher level language because the translation of a 
FORTRAN statement may result in many machine language instructions. This 
conversion is performed by a program called the FORTRAN compiler. The design 
of a FORTRAN compiler is definitely machine-oriented and is not part of the 
FORTRAN language. 


Offsetting this relative complexity (and consequent increased overall computer time), 
FORTRAN is a self-documenting language that cuts down the cost and time required 
for learning, writing, debugging, and maintenance. In standard FORTRAN there are 
approximately 30 statements to be remembered by the programmer, whereas, in 
assembly language there may be from 50 to 350 mnemonics. The compactness and 
self-documentation of FORTRAN are apparent in Figure 1-1. For instance, a 
program that requires from 10 to 20 pages of assembly coding may.require only 

one page of FORTRAN coding. 


FORTRAN ASSEMBLER MACHINE 





DIV*Y 770510040 


iO ‘ O: 
O STR*Z 772510050 


Figure 1+1, FORTRAN—Assembler—Machine Coding 


OQ Z =( (R+8-T)*X)/Y ENT*R 772110002, © 
O | ADD*s | 770110010 © 
O | suet = | 770210020 ©} 
O | MUL*X | 770310030 =O 
| 
| 


1.3. SOURCE AND OBJECT PROGRAMS 


1.4, 


A FORTRAN program written by the programmer represents a series of logical steps 
for the solution of a particular problem. This program is the source program. 


A source program must be translated to machine language for a particular data 
processing system. Translating the FORTRAN source program to machine language 
is generally accomplished by a prewritten program, the FORTRAN compiler. (The 
compiler is furnished with the data processing system and is not the responsibility 
of the programmer.) The output of this translation process (compilation) is the 
object program. 


COMPILATION AND EXECUTION 
The complete compilation process follows the steps below: 


(1) The source program is keypunched onto cards to produce the source program 
deck. 


(2) The FORTRAN compiler is read from auxiliary storage into the computer’s 
memory (main storage) as shown in Figure 1-2. 
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FORTRAN : MAIN 
COMPILER STORAGE 


AUXILIARY STORAGE 





Figure 1-2. Compiler to Memory 


(3) The source deck is then read into main storage and control is turned over to the 


compiler which determines the operations to be performed and generates the 
required machine language instructions to create the object program. This trans- 
lation process includes the assignment of memory locations for variables and 
constants, and the utilization of routines stored in auxiliary storage when required. 
The compiler produces an object program deck of machine instructions and a list- 
ing of the source program as shown in Figure 1-3. 


SOURCE 
PROGRAM 


DECK 





INPUT 


MAIN STORAGE 














OUTPUT OUTPUT 












SOURCE FORTRAN OBJECT 
PROGRAM COMPILER PROGRAM 
LISTING DECK 


Figure 1-3. The Compilation Process 


(4) The object program and internal data are entered into main storage. The data 


processing system can then execute the program: fetching data from input and 
auxiliary storage devices as required, operating on data as directed by the 
object program, and producing required output. 
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Steps 1 through 4 can be performed in a fairly rapid sequence referred to as ‘‘the 
compile and execute process’’ shown in Figure 1—4. The object program is placed 
in main storage as it is created by the compiler (rather than being transmitted in the 
form of punched cards as described in step 3) and immediately executed. 





COMPILATION EXECUTION 
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Figure 1—4. The Compile and Execute Process 
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1.5. THE COMPUTER SYSTEM 


A computer system is made up of hardware components (equipment) and software 
(operating system). The operating system is a program usually furnished by the 
computer manufacturer made up of routines that coordinate hardware activity and 
furnish various services to the user. 


1.5.1. Computer Hardware 


A digital computer is a data processing system that processes data in accordance 
with a set of instructions (program) and produces useful results. The programmer 
may regard the computer as a complex of devices with the functions shown in 
Figure 1—5, 


STORAGE 


INPUT OEVICE MEMORY OUTPUT DEVICE 


CENTRAL 
PROCESSING 
UNIT 
(CPU) 





AUXILIARY 


STORAGE 





Figure 1-5. Elements of the Computer 
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1.5.1.2. 
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Input Devices 


An input device furnishes the program (and possibly data) to main storage. The 
central processing unit (CPU) fetches the program and data from main storage, 
processes the data, and stores the result in main storage from where it is sent to 
an output device. If the program requires, input may be obtained from auxiliary 
storage devices and the result stored in auxiliary storage for updating or for 
transmission to an output device at a later time. Similarly, the input program 

can be held in auxiliary storage to be processed on demand by the CPU. 


Input devices are varied in nature. Currently, the most widely used device is 
the 80-column punched card reader and all descriptions presented here are 
oriented toward such a device. Other commonly used input devices are the 
console typewriter and punched paper tape readers. 


Main Storage 


Main storage retains the program, intermediate results, and output in addressable 
memory storage units. The contents of these storage units are transferred to the 
CPU for processing, and the storage units receive results from the CPU. Since 
processing in the CPU is generally much faster than the rate at which data is 
obtained from, and entered into, main storage, the processing speed of a com- 
puter is determined by the time required to enter and retrieve information into 
and from main storage. 


Central Processing Unit 


The central processing unit (CPU) decodes the machine language instruction, 
performs the arithmetic and logical processing functions indicated, and supplies 
the timing and control signals which synchronize the actions of the other 
elements of the computer. Each machine language instruction must indicate 

the operation to be performed and the location of the operands. 


Internal representation of values depends upon the particular computer; however, 
number representation is divided into two classes: fixed-point representation 

and floating-point representation. All computers can represent fixed-point numbers; 
however, some may uSe preprogrammed routines for the representation of floating- 
point numbers. : 


1.5.1.3.1. Fixed-Point Representation 


In fixed-point number representation, any external string of digits (optionally 
signed) is represented internally as a signed string of digits. A decimal point 
(radix point) is implied to be a fixed number of positions (this fixed number may 
be zero) from the rightmost position, hence the term ‘‘fixed point.’’ When arith- 
metic operations are performed on fixed-point operands, it is the programmer’s 
responsibility to keep track of the decimal point to ensure that both the operands 
and the result are correctly aligned. This is usually done by scaling the operands 
(multiplying by a power of 2 or 10). However, when many numbers of widely dif- 
fering magnitudes are involved, floating-point numbers are used because fixed- 
point representation requires an extensive programming effort for the scaling 
operations. The advantage of fixed-point representation is that any external 
string of digits can be represented exactly provided it is within the limits of 
range prescribed by the computer. This range provides limited representation 
when compared to floating-point representation. 
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1.5.1.3.2. Floating-Point Representation 


In floating-point representation, any number can be represented as a fixed- 
point number (mantissa) multiplied by an integral power (exponent) of 2 or 10, 
The exponent part is also a fixed-point number. 


The problems that arise from the use of floating-point representation of values 
are truncations and roundoffs that result from arithmetic operations. For example, 
the division operation 1.0/3.0 has as its true result, 0.333...3. However, since 
the mantissa can only contain a limited number of significant digits, only an 
approximation to the true result can be stored. Because of this, a floating-point 
value actually represents an approximation to the true value of a number. These 
errors are propagated as successive operations are performed upon floating-point 
operands and their results. With some computers another approximation is required 
if the internal representation of the mantissa is in binary form, rather than in 
binary coded decimal, because a decimal fraction cannot always be expressed 
exactly in binary form. For these reasons, comparisons involving floating-point 
numbers require special attention from the programmer. 


1.5.1.4. Auxiliary Storage 


Auxiliary storage is storage in addition to main storage and includes devices such 
as tape (magnetic and paper), disc, drum, and card (magnetic and paper). The time 
required to access data contained in auxiliary storage is significantly longer than 
for main storage; however, the cost of auxiliary storageis lower and capacity is 
generally larger. Auxiliary storage devices may contain files of data required for 
a program, different modules of the operating system, and/or results of a program 
that will later be transcribed on an output device. 


1.5.1.5. Output Devices 


An output device usually furnishes a visual display of the results, as specified 

in the program. In addition, an output device can inform the operator of various 
conditions arising during operation. The most widely used devices are the console 
typewriter, printer, and CRT (cathode ray tube) display. 


The auxiliary storage devices described in the previous paragraph can also be 
considered as output devices, depending on their end use. ; 


1.5.2, Computer Software 


The entire complex of software (programs furnished with the computer) is called 
the operating system. The operating system may contain programs which control 
scheduling, input/output, compilation, debugging, storage assignment, linking, 
loading, assembly, and other necessary functions depending upon the individual 
computer system. 


1.6. SAMPLE PROGRAM 


The following paragraphs present a simple executable program (a self-contained 
computing procedure) for the purpose of introducing FORTRAN programming and 
terminology. The concepts and terminology used in the description of the program 
are described in detail in other sections of the manual. 
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This program calculates the average of a series of numbers, each of which is supplied 
by a punched card. The program is general enough to calculate the average no matter 
how many values are involved. The last card of the data deck contains a value known 
to be outside the range of values expected. This card is used as an end data card. 
After the end data card is detected, the average is calculated and printed, together 
with explanatory text. - 


Figure 1-6 is 2 flowchart, with explanatory text, outlining the program. A flowchart 
should be constructed for any extensive program. It is a convenient means for detect- 
ing logical errors and provides documentation for other programmers who may be able 
to use or modify the same program. If the flowchart provides enough detail, the actual 
writing of the program is greatly simplified. 












INITIALIZE 

RUNNING COUNT AND 

RUNNING TOTAL TO 
ZERO 





READ IN 
DATA VALUE 


-” 
. 


{nitialize a running court and a running tatal to zero. The run- 
ning count will indicate the number of values (one value per data 
Card) to be averaged. Each time a value is obtained, it will be 
added to the running total. 


2. Read in a value from a card. 







3. Test the vatue for end data indicator. if the value is greater than 
9.0(108) it indicates that the previous value read in was the iast 
value to be used in computing the average, and the next step is 
step 7. if the vaiue is less than or equal to 9.0(103}, proceed to step 4. 


UPDATE 
RUNNING COUNT AND 
RUNNING TOTAL 
4. Add the input vatue to the running total. 


5. Increase the running count by 1: 


6. Read in the next value and repeat the processing steps by returning 
to step 2. 


7. Divide the running total by the running count ta compute the average. 
8. Print this average with some explanatory text. 


9. Indicate that there are no more instructions to. be executed. 







PRINT 
AVERAGE 
AND 
EXPLANATORY 
TEXT 


Figure T}-—6. Sample Prohlem Flowchart 
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Figure 1-7 is the FORTRAN program written from the procedure outline in Figure 
1-6. It is not the only program that could have been written for the problem, nor is it 


terms of lines required. It does introduce FORTRAN nomenclature and 
— the most important of which is the programmed loop. Such a loop 


is shown in the flowchart portion of Figure 1~6. The purpose of a loop is to repeat 

a series. of operations without respecifying the steps for each repetition. In conjunction 
with every loop, there must be some test that will make it possible to leave the loop; 
otherwise, it will be performed without end, being limited only by the time allotted to 


the program. In 
from the loop. 


5 16 


7 10 


"S* FOR COMMENT 
[rae irontean STATEMENT —_________» 


q) Q 


on eee 
ay fo, cum ttt ayt.tyz,6, oR 
w Li. | [pataL KOU 
() ke oe.tlalin, .vjat.u.e 
6) |, 3,0] |rea.o 1,1 
(7) | 10] [FORMAT ( 
(8) FOR, END 
ol... | toecvyatue 
co) lc, upolalte RUN! 
ant... | fro rage 
(12) KO. UN)T 
aa |e, RePejat Loop 
(4) | lo,0, ,TI0, ,3 
(15) elu te A,V,E 
(16) | |x.«,0,uIN.T 
(17) XM 6, AIN 
(18) T, AVE RA 


ayt., | tr 1 TE, ,3,, 


| le ORMAT (C1 


TiE, JA, VER 


G 


-° (20) 
(21) C,A,T,E] , END 
(22) | Is TOP 
(23) iC ate NO 
aot. | leno 





the sample program, detection of the end data card provides the exit. . 


20 30 


TjO, ,B,€, ,R,£,A,9 


DATA, |C, AROS 


A,GE, OF! NUMBERS FR,OM 
UN NING COUNT AND, RUNNING TOTAL 


NT... T.O,TjA,L./,0,,,0,.,0,/ 


0), VAL Ue 
G10.0) 
DATA ies 
CRE Le 
NG TOTAL 
TOTAL 


K OUNT 


OE 8) 
AND 
VALUE 


+ 1 


G.O 
RUN 


TO 20 
ING COUNT, 


+ 


OR Hae OO PS A OS 


AGE FIT, E.R, LAST, WAL UE 

K,0,U,N,T ; 
TOTAL, j/, XK, OUNT l 
—E AND |TEXT 

40), XM/E AN 
7H AVERAGE VALUE = , F10.5) 

O.F, ,£,Xj£,C,U,T,1,0,N 
MORE LINES 


Figure 1~7. Sample Program 
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Line 1 of the program is a comment line, indicated by the character C in column lL. ; 
This line (including the C) is printed when the program is compiled, but does not @& 





affect execution of the program. Comment lines provide documentation for the 
programmer. 


Line 2 is also a comment line with all blank characters producing a blank line 
(except that the C is printed). 


Each of the remaining comment lines applies to the statement(s) immediately follow- 
ing it. Again, comment lines are not required, but are included to aid the reader in 
following the program. 


Line 4 is a DATA or initialization statement that sets an initial value of zero for 
KOUNT (symbolic name for the running count) and TOTAL (symbolic name for the 
running total). The DATA statement does this at compilation time rather than 
execution time, in order to reduce execution time for the program. KOUNT will be 
a fixed-point value; TOTAL, a floating-point value. 


Line 6 is a READ statement that instructs the computer to read a punched card for 
the number to be represented by VALUE. 


Line 7 is a FORMAT statement that indicates the required number is found in the 
first ten character positions of the punched card. 


Line 9 is a logical IF statement. It tests the number that is read in for VALUE, if 

the number is greater than 9.0 (108), program control is transferred to the statement 

with a statement label of 20 (in columns 1 through 5). If the number is less than or 

equal to 9.0 (108), the next executable statement which follows in physical order 

is to be executed. The logical IF provides the exit from the loop. @ 





Line 11 is an artthmetic assignment statement that updates the running total. It 
obtains the current value of TOTAL, adds to it the number just read for VALUE, 


and assigns this sum as the new value for TOTAL. 


Line 12 is also an arithmetic assignment statement that increases the value of 
KOUNT by 1 each time a number for VALUE is read. 


-" Line 14 is a GO TO statement that completes the loop. When executed, it transfers 
control to the statement with statement label 30, so that the next card can be read. 


Line 16 is executed only after the end data card, which contains a value greater than 
9.0 (108), has been read. Its effect is to convert KOUNT to floating-point form to be 
compatible with TOTAL so that one can be divided by another (line 17). 


Line 17 is an arithmetic assignment statement that specifies the division necessary 
to calculate the average, XMEAN. 


Line 19 is a WRITE statement which specifies that the value of XMEAN is to be 
printed. 


Line 20 contains a FORMAT statement which specifies that the characters AVERAGE 
VALUE = are to be printed, followed by 10 print positions for the value of XMEAN, 
which is to be printed with five digits to the right of the decimal point. 
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Line 22 is a program control statement, the STOP statement, that terminates execution 
of the program. 


Line 24 is an end line, terminating compilation of the program unit. 


Note that every statement (columns 7 through 72) except the assignment statements 
starts with a keyword. (Comment lines are not considered statements.) The keyword 
is an English word that describes the purpose of the statement. Every statement in 
FORTRAN, except statement functions and certain assignment statements, begins 
with a keyword. Keywords are not reserved words in FORTRAN; they may be used 
anywhere in the program as symbolic names. ; 


Figure 1-8 is the complete deck of punched cards for the program in Figure 1-7. 

This deck consists of the FORTRAN program cards followed by data cards and an 

end data card. It should be clearly understood that punched cards are not the only 
means of primary input. Punched cards are shown here because of their widespread 
use and because they are easily used for illustrative purposes in examples. At least 
one system control card is required preceding the program, another preceding the data, 
and possibly an end of file card after the last data card. Requirements of these control 
cards depend upon the computer being used, and they are not part of the standard 
FORTRAN language. 


9.1EB ENO DATA CARD 30 
231.53 F Cid 
0.9876 28 
~48, 27 
~122.45 26 
12345.678 28 
ENO 24 







@ /C INDICATE NO MORE LINES 23 
© INDICATE END OF EXECUTION 

40 FORMAT (17H AVERAGE VALUE =,F10.5) 

WAITE(3, 40) XMEAN 

iC PRINT AVERAGE ANO TEXT 

XMEAN = TOTAL, XKOUNT 


. 20 XKOUNT = KOUNT 16 
C COMPUTE AVERAGE AFTER LAST VALUE 














‘C UPDATE RUNNING TOTAL AND RUNNING COUNT 
1F(VALUE .GT. 9.0€8) GO TO 20 9 
. 'C TEST FOR END DATA r) 
. 10 FORMAT (G10.0) 7 
30 READ(1.10) VALUE 6 
C OBTAIN VALUE 3 
DATA KOUNT, TOTAL/0,0.0/ 7 
/CNITIALIZE RUNNING COUNT ANO RUNNING TOTAL, 










COMPUTE AVERAGE OF NUMBERS TO BE READ. FROM DATA CARDS 


aan 1 oenena | aa i 

Oooo PCreFoosooccsccgeascgosecccoorscccnsoceagacsocrgaoccacdooggcngednacoacaas 
agar 

mit 





CPOE TONGS GROAN ANGAANABT AMS ONS HSE AADTLOKNDY DHMNTOMBMTAENN NNT OA 
cat SEED SURO RERREEREESEBOSEER! SEROGRERL Of [EI DERE SURE CCRROREREROGSRERSEERS | 
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2 G3333933392399993393322239IG9393239999333392292993GI3F33993933799933993939931933 
Ved [e CUSSCCECP CESS EC) | PECCSECOCEOCOL Cosel CECE EESL) CES EC ECOSOC OCs ere ee re ne) 
SSSSSSSSHSSQUPSSSQSSSSHSSSESSSSSSSPSSESSSSHSSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSHSS555S 
SESUCSECCSESECEEESUMSC GES SC OSUSES CESS ECU CUGCECCECESSEEE EES SESCCSG SSS EEC EEGGESS 
PERRI PRRRRRRES) LRORRRRRRRORPRAPRRSRRORERAR PORE ERE PRE ee Pe Penere Reese eeeeee renee! 
| SSOLESASECESEAEATAASTAACS KS TF CERES TATE AKSHSERSRTATRAT ARTES ZEGTHAH ESET 
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@ Figure 1-8. Sample FORTRAN Problem Deck and Data 
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2. WRITING A FORTRAN 
PROGRAM 


2.1.. GENERAL 


é This section discusses the organization of FORTRAN programs and the rules for 
writing a source program, with particular attention to constants, variables, and arrays. 


2.2. ORGANIZATION 


Organization of a FORTRAN program necessitates discussion in two areas: the con- 
cept of subprograms and the organization of program units. 


2.2.1. FORTRAN Program 


A FORTRAN program is made up of one, and only one, main program and as many 
subprograms as required. The main program contains the steps required to solve 
a given problem, the subprograms are subordinate program units used by the main 
program. Both are referred to as program units. 












UP-7536 
Rev. 1 





2 










FUNDAMENTALS OF FORTRAN 





SECTION: PAGE: 


The various types of program units (see Figure 2-1) are as follows: 





a A main program is a series of comments and statements which does not contain 
a FUNCTION, SUBROUTINE, or BLOCK DATA statement and is terminated by an 
end line. 


z A function subprogram is a series of comments and statements starting with a 
FUNCTION statement and terminated by an end line. 


w A subroutine subprogram is a series of comments and statements starting with a 
SUBROUTINE statement and terminated by an end line. 


m A specification subprogram is a series of comments and specification statements 
starting with a BLOCK DATA statement and terminated by an end line. 


MAIN 
PROGRAM 


SPECIFICATION 
SUBPROGRAM 


BLOCK DATA 


$ 





EXTERNAL PROCEDURES 
ee ee 


PROCEDURE SUBPROGRAMS 


external procedures 
en ee written in other than 


FUNCTION SUBROUTINE FORTRAN language 
SUBPROGRAM SUBPROGRAM 


. FUNCTION... SUBROUTINE... 


END 





Figure 2~]. Program Units of FORTRAN Program 
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Each program unit is independently compiled and then linked together by the operating 
system to form an executable program starting with the main program unit. Program 
units may be written in languages other than FORTRAN but must conform to the 

tules for FORTRAN subprograms. Such program units and procedure subprograms 
(function and subroutine subprograms) are termed external procedures. 


Execution of a program always starts with the first executable statement of the 
main program. In. the first example shown in Figure 2—2, the main program proceeds 
until it encounters a reference (call) to the external procedure. The external pro- 
cedure assumes control until it encounters a RETURN statement, which sends 
control back to the calling program unit (in this case, the main program). The main 
program then continues processing until another reference transfers control to the 
external procedure. The external procedure assumes processing until it encounters 
a RETURN statement (not necessarily the same as the first RETURN statement) 
and transfers control back to the main program. The main program then resumes 
processing until it encounters the STOP statement, which transfers job control to 
the operating system. 


The second example in Figure 2-2 shows how a procedure subprogram can call upon 
another procedure subprogram during execution. 


A specification subprogram (with a BLOCK DATA header) consists entirely of 
nonexecutable statements and therefore never assumes control during execution. 


Subprograms are useful because they eliminate repetitive coding of procedures used 
many times in a program. In addition, a library of mathematical external procedures, 
called basic external functions, is present which contains debugged procedures for 
computation of mathematical functions such as square root, sine, etc. The main 
program of a large FORTRAN program can be coded as a logical skeleton consisting 
primarily of references to subprograms; these subprograms can be independently 
coded and compiled concurrently with the main program. 
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EXAMPLE | 





: PROCEDURE 
MAIN PROGRAM SUBPROGRAM 


headers 


irst executable statement first executable statement 


RETURN 


RETURN 





EXAMPLE 2 


PROCEDURE PROCEDURE 


MAIN P 
AIN PROGRAM SUBPROGRAM SUBPROGRAM 


first executable statement first executable statement 





NOTE: A header is either a FUNCTION or a SUBROUTINE statement. 


Figure 2-2. Controf Path During Execution 
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2.2.2. Program Unit Organization 


A program unit consists of comments, statements, and one end line. A FORTRAN 
statement falls into one of two categories: an executable statement or nonexecutable 
statement. An executable statement specifies an action; a nonexecutable statement 
describes the characteristics and arrangement of data, editing information, statement 
function definitions, and classification of program units. Nonexecutable statements 
are generally intended as instructions to the compiler; no executable machine language 
instructions are generated. Executable statements result in executable machine 
language instructions, effective at execution time. As an example, lines.6 and 7 of 
Figure 1~7 are as follows: 





The READ statement is an executable statement specifying an action; the FORMAT 
statement is a nonexecutable statement specifying the characteristics and arrange- 
ment of data. Table 2-1 is a precise guide for determining whether a specific 
FORTRAN statement is executable or nonexecutable. The order (sequence) of 
statements within each program unit is shown in Table 2-2. 


STATEMENT GENERAL CATEGORY 
EXECUTABLE STATEMENTS 


arithmetic assignment statement 
logical assignment statement 
GO TO assignment statement 












assignment statements 


GO TO statements 
'F statements 
CALL statements 
CONTINUE statement 
RETURN statement 
STOP statement 
PAUSE. statement 

OO statement 









control statements | 


READ statement 
WRITE statement 
REWIND statement 
BACKSPACE statement 
ENOFILE statement 









10 statements 


Table 2~1. FORTRAN Statements (Part 1 of 2) 
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STATEMENT 
NONEXECUTABLE STATEMENTS 


DIMENSION statement 
COMMON statement 
EQUIVALENCE statement 
EXTERNAL statement 
type-statements: 

INTEGER statement 

REAL statement 

DOUBLE PRECISION statement 

COMPLEX statement 

LOGICAL statement 


DATA statement 

FORMAT statement 
statement function definition 
FUNCTION statement 


SUBROUTINE statement 
BLOCK DATA statement 
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GENERAL CATEGORY 









specification statements 









data initialization statement 
format statement 


function defining statement 














subprogram statements 


NOTE: The end line is nonexecutable and is not considered a statement. 


PROGRAM UNIT 


MAIM 
PROGRAM 


PROCEDURE 
SUBPROGRAM 


SPECIFICATION 
SUBPROGRAM 


Table 2-1. FORTRAN Statements (Part 2 of 2) 


ORDER OF STATEMENTS 


specification statements and‘FORMAT statements, in any combination 


statement function definitions, DATA statements, and FORMAT state- 
ments, in any combination 


executable statements, FORMAT statements, and DATA statements, 
in any combination 


end line 


Minimum requirements: an executable statement and one end line. 


FUNCTION or SUBROUTINE statement 


specification statements and FORMAT statements, in any combination 


statement function definitions, DATA statements, and FORMAT state- 
ments, in any combination 


(4) executable statements, FORMAT statements, and DATA statements, 
in any combination 


(5) end fine 

Minimum requirements: for function subprogram — one FUNCTION statement, 
an executable statement, a RETURN statement, and an end !ine; for sub- 
routine subprograms — one SUBROUTINE statement, a RETURN statement, 
and an end line. 

(1) BLOCK DATA statement 

(2} specification statements (except EXTERNAL), in any combination 

(3) DATA statements 

(4) end line 


Minimum requirements: one BLOCK DATA statement, a COMMON statement, 
a DATA statement, and one end line. 





Table 2-2. Ordering of FORTRAN Statements 
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CHARACTER SET 


The FORTRAN character set (see Table 2—3) consists of all the characters required 
for writing FORTRAN statements. 


LETTERS A through Z (only upper case letters) 


ALPHANUMERIC { 


DIGITS Q through 9 


blank, represented by B in text and 
coding form 
equals 
pilus 
minus 
asterisk | 
SPECIAL CHARACTERS / slash 
left parenthesis 
right parenthesis 
comma 
decimal point 


currency symbol 


Table 2-3. FORTRAN Character Set 


In addition to this set, each processor has its own set of additional characters. The 
FORTRAN character set and any additional characters in the processor set are 
referred to collectively as the processor character set. 


FORTRAN PROGRAMMING FORM 


A typical FORTRAN programming form (for 80-column punched card input) is shown 
in Figure 2-3. A FORTRAN line uses only columns 1 through 72; the information in 
columns 73 through 80, shown only in the program listing, may be (depending upon 
the processor) used for documentation. (In Figure 1—7 these columns are used to 
serially identify the lines of the program.) Columns 1 through 5 are used for state- 
ment labels; column 6 is used to indicate the continuation of a statement; columns 
7 through 72 contain the statement (or its continuation) or an end line. 


In writing a FORTRAN program, dlank characters may be used freely to improve read- 
ability with these general exceptions: comment lines, the end line, and Hollerith 

data. Except where noted, blank characters have no meaning and are ignored. Exceptions 
to the general rule are explained in detail in the applicable paragraphs of the manual. 


Comment Line 


A comment line is printed only when a program unit is being compiled and is a 
convenience for the programmer. Only characters from the FORTRAN character set 
may be used. A comment line is indicated by the character C in column 1 (as shown 
in Figure 1-7). Any blank character in a comment line will be printed as such. A 
comment line may not be used between the initial line of a statement and its con- 
tinuation lines, and must always be followed by an initial line, another comment 
line, or an end line. No executable code is generated for a comment line. 
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2.4.2. End Line 





One end line must appear as the last line of each program unit to be compiled. Its 
purpose is to indicate the end of compilation for a particular program unit. Columns 
1 through 6 must be blank; columns 7 through 72 must contain the characters END. 
The characters may be interspersed with, preceded by, and followed by blank 
characters. Thus, starting from column 7, i 


END 
E N D 
EN D 


are all valid end lines. 


UNIVACG FORTRAN 


PROGRAMMING FORM 
PAGE. 


PROGRAM PROGRAMMER. DATE 











“c* FOR COMMENT 


srAmeca | HEORTRAM STATEMENT————__ : 
1 | a See tiie a SSS ee £ | een 

















Figure 2-3. FORTRAN Programming Form 
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2.4.3. Statements 

@ A statement consists of an initial line and, if required, up to 19 continuation lines 
which follow in sequence. Each line of the statement is written in columns 7 through 
72. An initial line must have either the digit 0 or a blank character in column 6 and 
must not have the character C in column 1. Each continuation iine must have a 
character other than the digit 0 or blank character in column 6, and must not have 
the character C in column 1. Columns 1 through 5 are available to the programmer 
for documentation (except for the character C in column 1). Each continuation line 
must be immediately preceded by an initial line or another continuation line. 


For example, line 20 in Figure 1-7 was written: 


"“C* FOR CONMENT 


"NUMBER: moet Jronrean STATEMENT» . ; 


Po. wo) [oR WALT OZ, A VIEIRA GE, VALUE! |= ,,F,1,0,.,5) 


It could be written as: 


| 


| 4 ol [FoR MIA a 
Ve cise hi AVERAGE VALUE te 


Pos (Bh F101. 5) 


However, it can not be written as 





because the 17H of the FORMAT statement means that the 17 characters immediately 
following the H are Hollerith data. As was pointed out in 2.4, blank characters are 
significant in Hollerith data. 


Frequently, the character in column 6 of a continuation line indicates the sequence 
of continuation lines (it is hoped that what is described as a common practice in 
this manual is not deemed mandatory by the programmer). For example, the first 
continuation line might contain the digit 1 in column 6; the second, the digit 2; and 
so on. Also, since the first five columns of continuation lines are available for 
documentation (using any characters from the FORTRAN set), these are sometimes 
used to contain the statement label of the initial line. 
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2.4.4. Statement Labels @ 
A statement label is an unsigned integer (1 through 99999) that identifies a FORTRAN 
Statement and is written in columns 1 through 5. Only the digits 0-9 and blank charac- 
ters may be used in a statement label. The same statement label cannot be used more 
than once in a program unit; the same statement label can appear in more than one 
program unit. The value of a statement label does not affect the order in which 
statements are executed. 


A statement label is meaningful only when used with the initial line of a statement; 
it is ignored when used with continuation lines. It enables other statements of the 
program unit to reference it; therefore, it is superfluous when used with a statement 
that is not referenced by another statement. Some processors will print out a warning 
diagnostic message that indicates an unreferenced statement label. The possibility 
of the programmer having forgotten to insert the reterenicing 3 statement is strongly 
indicated. 


A maximum of five digits may be used in a statement label. All blanks and leading 
zeros are neglected. 


For example, in the sequence 


"Cc" FOR COMMENT 
STATEMENT 


NUMBER a take AER 
7 


30 





— 
wn 
“NI 
ail 
u 
> 
+ 
o 
t 
a 
a 
9 


rae GREG ECR ecm 
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because control is being transferred to statement label 157, not 1570. 


2.5. FORTRAN DATA 


arithmetic, logical, or Hollerith. Arithmetic data is used in computations restricted 
to numbers. Logical data is used to indicate whether a specific condition is true or 

- false. Hollerith data is information to be used literally; it may contain any and all 
@ characters of the computer set of characters and is generally used for printing 


| 
Data is information manipulated by a program. An item of data is classified as 


messages, titles, and headings. 


2.5.1. Data Types 


The various forms in which data may appear are known as types. The term type, 
when applied to data in FORTRAN, has a special meaning. An item of data, if 
arithmetic, must be integer type, real type, double precision type, or complex 
type. If not arithmetic, it must be logical type or Hollerith type. The different data 
types are discussed in the remainder of this section. 


One important attribute of a data type is the number of storage units required for 
internal representation of a datum. This attribute is of special importance in the 
use of COMMON and EQUIVALENCE statements, and is listed in Table 2—4 for 
reference purposes. 


DATA TYPE STORAGE UNITS* REQUIRED 


integer 


real 


double precision 


compiex 


logical 





«The term ‘‘storage unit’ is defined differently for different computers. In one computer, it may 
@ mean a word; in another computer, it way mean five consecutive bytes, and so on. 


Table 2-4, Memory Requirements for Data Types 
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An item of data can be expressed as a constant or can be represented by a symbclic 
name. For example, in 





" FOR COMMENT 
STATEMENT 


NUMBER Jrowtean 2 aia eae A 
fi; 10 30 
A = B + 5,.,9 


A and B are symbolic names used as variables and 5.0 is a constant. In the DATA 
statement 


T NiALM SH TUT LE t rig Lod 


NAME is a variable; SHTITLE is a Hollerith constant representing the Hollerith 
datum TITLE. 


2.5.1.1. Integer Type 


An integer type datum is always the exact representation of an integer (a string 

digits without any decimal point). It may assume positive, negative, or zero 

values (the value zero is considered neither positive nor negative). Each compu- & 
ter places a limit on the number of digits (or the maximum absolute value) that 

may be contained in an integer type datum. An overflow condition occurs when 

this is exceeded. The term ‘‘exact representation’’ is used because an integer 

type datum is represented internally as a fixed-point number. 


2.5.1.2. Real Type 


A real type datum is the processor approximation to the value of a real number 
(as opposed to a complex number) which may or may not have a decimal point. 

-" It is limited to a specific number of significant digits and a specific range of 
values, both of which differ for different processors. The real type value may be 
positive, negative, or zero (the value zero is regarded as neither positive nor 
negative), When the range is exceeded fora value, an overflow condition occurs; 
when the value is less than the lower limit of the range and is not zero, an under- 
flow condition exists. A real type value is an approximation because it is repre- 
sented, intemally, in floating-point form. 
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2.5.1.3. Double Precision Type 


A double precision type datum, like a real type datum, is also the computer 
approximation of a real number, but the degree of approximation is greater than 
that of a real type datum because it can contain more significant digits. The 
number of significant digits is not necessarily double that of a real type value; 

it will be greater than that of a real type value, and this number must be specified 
for each processor. In addition, its range of values may be the same as or greater 
than, rarely less than, that of a real type value. This range must also be specified 
for each processor. A double precision type datum may be positive, negative, or 
zero. Here, as for all arithmetic values, the value zero is treated as neither 
positive nor negative. 


2.5.1.4. Complex Type 


A complex type datum is the processor approximation to the value of a complex 
number. It is represented as a pair of real type values; the first of the pair 
represents the real part of the complex value, the second of the pair represents 
the imaginary part of the complex value. Each part, therefore, has the same degree 
of approximation as a real type datum. 


2.5.1.5. Logical Type 


A logical type datum is an item that may assume only the values true or false, 
There is no standard internal representation for these values; internal representa- 
tion depends upon the processor. 


2.5.1.6. Hollerith Type 


A Hollerith type datum is a string of characters taken from the processor set 
(which includes the FORTRAN set) of characters. The blank character is a sig- 
nificant character in Hollerith type data. 


2.6. CONSTANTS 


A constant is a value (arithmetic, logical, or Hollerith) that is defined by its appear- 
ance in a DATA or an executable statement. Once defined, it must not be redefined 
during execution of the program. For example, in the arithmetic assignment statement 





the integer 300000 is an arithmetic constant. 


< 


Throughout this manual the term ‘“‘constant’’ alone, when applied to an arithmetic 
constant, means a constant without a sign. An integer, real, or double precision 
constant is said to be signed when it is written immediately after a plus or minus 
sign. An optionally signed constant may have the plus sign omitted; the minus sign 
cannot be omitted for a negative value. 
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In the following discussions it is arbitrarily assumed that an integer type datum can ~ 
consist of up to nine digits and a sign; a real type datum, up to eight significant digits; @ 
and a double precision type datum, up to 15 significant digits. It is also assumed that 

one storage unit can contain five Hollerith characters. These assumptions are made 

for use in examples. These specifications differ for each processor and the programmer 

should consult the reference manual for the particular processor being used. 


2.6.1. Integer Constant 


2.6.2. 


An integer constant is a nonempty string of decimal digits. No decimal point, comma, 
alphabetic, or special character (except blank characters) may appear anywhere in 
the string. Blank characters and leading zeros are ignored except where noted. The 
value 0 is a valid integer constant. 


Examples: 

0 valid integer constant 

0000 valid integer constant 

123000123 = valid integer constant 

123,456 invalid integer constant because of comma 

“0 invalid integer constant because of decimal point 

0. also invalid integer constant 9ecause of decimal point 
bbobbb invalid integer constant because the string must not be empty 


Real Constant 


A real constant may be written in any one of the following three ways: 





(1) As a basic real constant. A basic real constant consists of an integer part, @ 
a decimal point, and a decimal fraction part, in that order. Both the integer 
and the decimal fraction part are strings of decimal digits and (possibly) 
blank characters. Either part may be empty (consisting of only blank characters), 
but not both. Except for this restriction, any blank characters in the constant 
are ignored by the processor. 


Examples: 
.0 valid real constant 
0. valid real constant 
0.0 valid real constant 


123.45678 valid real constant 
000123.45678 valid real constant. Leading zeros are ignored. 


123.456789 valid real constant. The presupposed limit of eight significant 
digits is exceeded, but the processor can approximate the 
written constant by an appropriate real type value. 


bb. 66 invalid real constant because integer part and decimal fraction 
part are empty 
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ALGEBRAIC VALUE FORTRAN EQUIVALENT 


0.1234°1074 


~0.1234°104 


2 


SECTION: PAGE: 


As a basic real constant followed by a decimal exponent. A decimal exponent 
is the letter E followed by an optionally signed decimal integer. This decimal 
integer has a range unique to each processor, usually two or three digits, with 
leading zeros ignored. The decimal exponent-has the effect of multiplying the 
preceding constant by a power of 10, as specified by the optionally signed 
decimal integer. , 


Examples: 
.0E0 valid real constant 
1.23E2 valid real constant. It could also have been written as a basic 


real constant, 123.0 
123456.01E—2 valid real constant. It could also have been written as a basic 
real constant, 1234.5601 ; 
123456,.01E+05 valid real constant. It could also have been written as the 
basic real constant 12345601000.0. The zeros trailing the 
rightmost 1 are not handled as significant digits, but as 
positional digits. 


There is no assurance that equivalent mathematical versions of a value will yield 


identical processor values; these will be equivalently approximated by the processor. 


As an integer constant followed by a decimal exponent. 
Example: 


12345601E—4 valid real constant. It is equivalent to the real constant 
123456.01E—2 or to the real constant 1234.5601. 


Figure 2-4 shows different ways of stipulating the values 0.1234 - 10~4 and 
~—0.1234 + 104 as real constants in FORTRAN statements of a program. 


jO,.,1,2 3,4,E- 4 | = 
ut See a ee eee eee 
1142,3.4,E,-,0.8, , |]... , \ ee ly rey 
,o.- 000 0 12,34) : | , ' 

Bi IO es Ce ER OR OLAS OOO VO a 
1-,0.. 12,3,4,E.+,4 | 
Lo-11,2,3,4,E,+,0.4) | ' pt 
1-,1,2,3,4,E, 0 | 

-~ 1234. xy a | 


Figure 2-4. Real Constants in FORTRAN Storements 
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A double precision constant may be written in either of the following ways: 


(1) As a real part (of exactly the same form as a basic real constant) followed by 
a double precision exponent. A double precision exponent is similar to a 
decimal exponent, except that the letter D is used instead of the letter E. 


Examples: 


12345678.9012345D0 valid double precision constant 
12345678.9012345D—5 valid double precision constant 
123.0D3 valid double precision constant 

(2) An integer part followed by a double precision exponent. 


Example: 


123456789012345D-—5 valid double precision constant 


Complex Constant 


A complex constant is the processor approximation of a complex number. It is 
written as a pair of optionally signed real constants enclosed in parentheses. The 
first of the pair is. the real part. This is followed by a comma and then the second 
of the pair, which is the imaginary part. 





Examples: 
(1.2 , 3.4) valid complex constant. It represents 1.2 + 3.41. ©® 
(9.0 , 3.4) valid complex constant. It represents the imaginary 
number 3.47. 
(-34E-2, —4.5E3) valid complex constant. It represents ~0.34—4500:. 
(-3400, —4.5E3) invalid complex constant. Both constants of the pair 


- must be real type constants. 
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2.6.5. Logical Constant 


A logical constant is written as either . TRUE. or .FALSE. The periods are required, 
as shown. It represents a condition as either true or false, respectively. An analogy 
is a switch which is either on or off, with . TRUE. corresponding to on, and .FALSE. 
corresponding to off; then, in a program, the programmer can test the setting of the 
switch and condition the sequence of execution accordingly. In the following 
sequence, a logical IF statement is used for testing. 


Peers it OG)! ICAL, ,€,VE,N 
Berea li 


E.ViEWNI 1= «TRUE. 


| | UNF, CeWV,E.N,) 6,0 TO 1100 


The first statement is a type-statement which informs the compiler that the variable, 
EVEN identifies a logical value, that is, either true or false. The second statement 
is a logical assignment statement that sets EVEN to true. The third statement is a 
logical IF statement that means: IF EVEN is true go to statement 100; otherwise, 
go to the next statement (statement 10). 


1 


= 
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2.6.6. Hollerith Constant 


A Hollerith constant is written as the letter H preceded by an integer constant 
which specifies the number of characters following the H that are part of the 
Hollerith constant. A Hollerith constant may appear only in a DATA statement 
and in a CALL statement. 


- Examples: 

SHTITLE valid Hollerith constant 

6HIS NOT valid Hollerith constant. The space character must be counted. 

SHA=$55 valid Hollerith constant. All characters from the processor 
character set are permissible (see 2.3). 

NOTES: 


(1) Although Hollerith constants can appear only in a DATA or CALL statement, 
Hollerith data may appear in input/output data and in the FORMAT statement. 


(2) The internal representation of a Hollerith value may be different from the 
internal representation of an arithmetic value. For example, 1H1 need not 
have the same intemal representation as the arithmetic value 1, and, if they 
are compared, they will generally be treated as unequal. The FORTRAN 
language does not define any correspondence between Hollerith values and 
arithmetic or logical values. ; 


2.7. SYMBOLIC NAMES 


A symbolic name consists of one to six alphanumeric characters (any blank charac- 
ters are ignored, other symbols are prohibited), the first of which must be a letter. 
A symbolic name is followed by a subscript only when used as an array element 
reference. In a program unit, a symbolic name may identify an item in one (and 
usually only one) of the classes (array and array elements are considered the same 
class) shown in Table 2—5. The remainder of this section discusses the use of 
symbolic names for variables, arrays, and array elements. The programmer is 


entirely free in the choice of words for symbolic names. Keywords such as FORMAT, 
GO TO, READ, etc., are not considered reserved words; these may be used as symbolic 


names, or the first letters of a symbolic name anywhere in the program. 





PAGE: 
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A SYMBOLIC NAME IS CALLED TYPE TEXT 
THAT IDENTIFIES A(N) ASSOCIATIGN?* REFERENCE 
A(N) ; 
variable variable : es see 2.7.3 
array array name i see 2.7.4 


array element array element see 2.7.4 
reference 


statement function statement see 8.2 
function name 


intrinsic intrinsic see 8.3 
function function name 


subroutine subroutine name no see 8.6 


external function external function yes. see 8.5 
name 


external procedure external procedure not necessarily see 8.7 
name alone as an name 
actual argument 


block block name see 8.8 





*Type association means that the data type is either implied in the name or explicitly 
specified (see 2.7.2). 


Table 2~+5. Uses of Symbolic Names 


2.7.1. Uniqueness of Symbolic Names 


A symbolic name cannot be used in a program unit for more than one purpose except 
as noted in the following rules. 


A block name can be used as an array name, variable, or name of a statement 
function in the same program unit. 


In a function subprogram, the function name that follows FUNCTION in the header 
f the subprogram must also appear as a variable in the same program unit. 


Once a symbolic name is used as an external function name, subroutine name, 
procedure name, or block name in any program unit, it may not be used anywhere 
in the program for the same purpose except as originally used. 


A symbolic name that appears in an EXTERNAL statement of a program unit and 
is used only as an actual argument within that program unit must be the name of 
an external function or subroutine in the same program. However, that same name 
may be used in any other program unit for a valid purpose. 


A symbolic name used to identify an intrinsic function in one program unit may be 
used for any purpose in another program unit. 
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2.7.2. Typing of Symbolic Names 





It was indicated in Table 2—5 that certain uses of symbolic names have data type & 
association. This means that the compiler must be informed of the data type for 

the symbolic name, that is, integer type, real type, double precision type, complex 

_type, or logical type. There is no mechanism in FORTRAN for specifying Hollerith 

type; data of this type, other than Hollerith constants, are identified under the 

guise of a real, integer, or logical type variable or array element, as will be shown 

in one of the following paragraphs. 


The data type of a symbolic name can be explicitly declared or, if not explicitly 
declared, implied by the first letter of the symbolic name. A constant automatically 
defines its own type by the form of its appearance. 


2.7.2.1. Explicit Type Declaration 


Explicit type declaration means that the symbolic name is declared in a type- 
statement. The one exception for explicit type declaration is the function name 
that appears in the FUNCTION statement of a subprogram. In the function sub- 
program, this symbolic name can only be explicitly typed in the FUNCTION 
statement. 


“C* FOR CONMENT 


“ranee | JFORTRAN STATEMENT» 
Ho 5 (617 10 30 


P| feeb) OT AL, |b/ALMi8,0,A 


1D,O U BILE, PRE CIS IjO.N BETA, ., GAMMA 
| levowea ex ALP HLA A.2,4, | 


INTEGER, S|GMA 
LOGIICAL, XXX, .., 1Y, YY, KKK 





This group of five statements consists of type-statenents. The first statement 
specifies that IOTA and LAMBDA represent real type data, the second. that 
BETA and GAMMA represent double precision data; and so on. 


2.7,2.2. Implied Type Declaration 


Implied type declaration means that the data type associated with a symbolic 
name is implied by the first letter of its name. unless the data type of the 
symbolic name is declared explicitly. If the first letter of the symbolic name is 
any of the six letters I, J, K, L, M, or N, the compiler will assume integer data 
type for that name; if the first letter is any of the other alphabetic characters 
the compiler will assume real data type for that name. 
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This rule for implied data type association clearly indicates that symbolic names 
representing double precision, complex, or logical data must always be typed 
explicitly. An exception to this mle is that references to intrinsic or basic 
external functions (see Section 8) do not require explicit type declaration in the 
referencing program unit. The function names involved already have the data 
type of the function names known to the processor. 


Examples: 


a ee ee eee ee 
Preeak B PR (8) V2 V4 Vi6 
| | hoe ean TEST, ,. 


C,0,.M,P{L,E,X, .N.M,B,R,5,, |NM,B,R6 re Te Les ey eae 


| dK AD Dy =, NMB.R.2, +) NM BRA 


Po | detest VAL UL El= VALUE, s+ 1.19 
is sis | 


re =| ,Vi4: - Vi6 ' 
N MB, R16 = .(,12.,3,41,,4,5..6,7,),+N MBR 5 


KADD is the sum of two integer values (all implied type declarations). In the 
next statement, if the logical value of TEST (explicitly declared) is true; 
VALUE is increased by one (implied real type). In the statement that follows, 
the double precision V2 is formed by the subtraction of V6 from V4 (all explicitly 
declared as double precision). In the last statement, the complex value NMBR6 
is formed by adding 12.34 to the real part of NMBRS and 45.67 to the imaginary 
part of the complex type NMBRS (explicitly declared). 
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2.7.2.3. Hollerith Values 





It has been stated that there is no mechanism in FORTRAN for specifying 
Hollerith type, as such, and that it must be guised as a real, integer, or logical 
type. The following example shows how this can be done. 


S* FOR COMMENT 


Sarma ieee STATEMENT —————__p> 


a ee} _ rae oo) 


| loar, tre VALUE 6.6.0. / 54 N AMES, 4HN ONE, 3H YES / 


The symbolic names LGC, ITEM, and VALUE each can represent five Hollerith 
characters (from the original assumption in 2.6). The variable ITEM contains 

the characters NAMES; the variable VALUE, the characters NONE followed by 

a blank character; the variable LGC, the characters YES followed by two blank 
characters. It is the programmer’s responsibility to remember that LGC, ITEM, 

and VALUE contain Hollerith data; there is no way for the processor to ‘‘remem- 
ber’’ this. In FORTRAN, there is no provision for moving Hollerith data from one 
internal location to another, or otherwise processing the data (as in comparisons). 
The Hollerith data from the previous example can be used in the following sequence: 





! 
P| eas Muact) 
to} |r o.Rw TiGUuy, SCA 51%) 
aa el 


i 





This will cause a line to be printed as: 


“BN AMESBN ON EbBY ES55B 


The AS in the FORMAT statement informs the processor that the contents of 
ITEM, VALUE, and LGC are to be treated as Hollerith data. 


If the variables ITEM and VALUE are used in arithmetic statements, or the 


variable LGC in a iogical expression, the results are unpredictable in standard 
FORTRAN and will vary, depending upon processor implementation. 
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Variables 


A variable is a symbolic name that identifies a single value. The rules for de- 
claring the data type of a variable have already been discussed. Once the value 
for a variable has been defined, it may be redefined as many times as required 
by the program. It may be defined by a DATA statement (Figure 1-7, line 4), 

by an input/output statement (Figure 1-7, line 6), by an assignment statement 
(Figure 1-7, line 17), by its use as an argument in a subprogram reference 
(discussed in Section 8), by its use as a DO statement or DO-implied list control 
variable, or by association in a COMMON or EQUIVALENCE statement. A variable 
cannot be used unless its value has been defined at least once. For instance, in 
Figure 1~7, if the original value of KOUNT had not been defined as zero, its 
use as a running count would have been meaningless. The same is true for the 
variable TOTAL used as a running total. Failure to define a variable used for 
such purposes is quite common and leads to unpredictable results. 


Arrays 


An array is an ordered set of values, each of which is called an array element. The 
entire set of values is identified by a symbolic name called the array name. All 
elements of the array must have the same data type, which is determined by the 
data type of the array name. 


The use of array elements in FORTRAN corresponds to the use of subscripted 
variables in ordinary algebra. For instance, the algebraic expression a ~ b—c 
uses simple variables which correspond to the use of variables as defined in 
FORTRAN. However, the algebraic expression aj — ag — a3 uses subscripted 
variables with one subscript; the algebraic expression 141 - b1 2 - bi 3 uses 
what is called in FORTRAN two-dimensional array elements; the expression 

Cy 4,1 ~ S1,1,2 ~ £1,1,3 uses three-dimensional array elements. The first expres- 
sion is written in FORTRAN as A(1) — A(2) — A(3), for the one-dimensional array 
named A; the second is written as B(1,1) — B(1,2) — B(1,3), for the two-dimensional 
array aamed 8B; and the third is written as C(1,1,1) — C(1,1,2) — C(1,1,3), for the 
three-dimensional array named C. An array, therefore, is nothing more than an 
ordered list, any element of which can be identified by the array name followed by 
the appropriate subscript. : 


If it is assumed that the maximum value of each dimension is 3, then the order of 
the elements in array A is: 


A(1) 
A(2) 
A(3) 


The order of elements in array B is: 


B(1, 1) 
B(2, 1) 
B(3,1) 
B(1,2) 
B(2,2) 
B(3,2) 
B(1,3) 
B(2,3) 
B(3,3) 














UP-7536 awe 5 






FUNDAMENTALS OF FORTRAN 


SECTION: PAGE: 





The order of elements in array C is: 





Note that the leftmost integer of the subscript varies most rapidly in the order 
of progression. 


The total number of storage units required for an array is the sum of the storage 

units required for its elements. If the data type of array name A is integer, real, 

or logical, array A will require 3 storage units (see Table 2—4); if array name A 

is double precision or complex, it will require 6 storage units. Under the same 

conditions, array B will require 3 - 3 or 2(3 ° 3) storage units; array C will require a 
3°3.°3 or 2(3 3 - 3) storage units. 


In a simple situation, visualize a book with the title NAME consisting of only 
columns of statistics on each page. Let the leftmost subscript represent the row; 
the next subscript, the column; and the last subscript, the page. Then, to obtain a 
Single statistic, reference the desired statistic by title, row, column, and- page. An 
example is NAME (3,4,5), referring to the third row of the fourth column of the fifth 
page of the book called NAME. This type of reference is called an array element 

- reference (sometimes called a subscripted variable). 


2.7.4.1. Array Declaration 


Before an array or any of its elements can be referenced.in an executable state- 
ment.of a program unit, the array must be declared. This means that the name, 
data type, and maximum number of array elements required must be stipulated. 


An array can be declared in a DIMENSION statement, COMMON statement, or 
type-statement. The form of the declarator is 


v (1) 


where: vis a symbolic name called the array name, and 1, the declarator sub- 
script, is composed of 1, 2, or 3 expressions, each of which may be either an 
integer constant or an integer type variable, except that a COMMON statement 
does not permit variables. 
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The parentheses are required as shown. Each expression is separated from the 
next by a comma. If any expression is an integer type variable, the array is known 
. as an adjustable array. 


The data type of the array name is governed by the rules for implied and explicit 
ty ping of symbolic names. 


Examples: 


Q PRE CUS TOM ViCT R13, 14, 5,) T.R.X,(;2,,,3,),. 8,7, AF, XC, 3,5,),, MANY 
1M EIN, S)1,0,N 1S, T7120, 415.236.) 
O,Gi1 Ayh, ,b,1,8, 7.2 


| eos M10, MAW, Yi (123,134) 








D1 MEINS!,0,N, 14,7, GIR, S,(,K,2,.4,, K 5) 


M D YPLENIG : i 


The first statement shows array declaration with a type-statement. The arrays 
VCTR, MTRX, and STATX are double precision arrays: VCTR is a three-dimen- 
sional array requiring 120 storage units; MTRX is a two-dimensional array 
requiring 12 storage units; STATX is a one-dimensional array requiring 30 storage 
units. 


The next group of three lines shows two different ways of declaring an array. The 
first is array declaration with a DIMENSION statement and a type-statement; the 
second shows array declaration in a COMMON statement of the double precision 
array MANY. 


The last group of lines shows array declaration for a three-dimensional adjustable 
array of integer type, INTGRS. Such a declaration (with integer type names as 
subscripts) can appear only in a procedure subprogram. The values for K2 and Ké 
will be furnished at execution time when the subprogram is referenced. 


In a program unit, any appearance of an array name must be followed immediately 
by a subscript except, possibly, in any of the following cases: 


® in the I/O list of a READ or WRITE statement; 

@ in the list of dummy arguments of a FUNCTION or SUBROUTINE statement; 

m in the list of actual arguments supplied in the reference to an extemal procedure; 
@ in a COMMON statement; 

@ ina type-statement; 


@ in a formatted READ or WRITE statement, to designate the format specification. 
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2.7.4.2. Array Element Reference @ 


An array element is one of the items in the set that constitutes the array. An 
array element is identified by an array element reference of the form: 


array name (se) 
or 

array name (sel,se2) 
or 


array name (sel,se2,se3) 


where array name must be the same as that in the array declaration and each se 
is a subscript expression. Parentheses and commas are required as shown. 


A subscript expression is any of the following forms. In these forms, the single 
asterisk is the FORTRAN symbol for multiplication. 


+k 
-—k 
*v 
*vitk 


*y—k 


a 9 4 ¢ © 8S 94 


Each c and k is an unsigned integer constant; each v is an integer variable name. : Se 


Examples: 
13 
NMBRI1 
NMBRI1 + 13 
NMBRI — 13 
2 * NMBR1 + 13 
2 * NMBR1 — 13 


Note that the order of the different elements must follow the order described in 
the previous paragraph. As an example: 2 * NMBRI + 13 is a valid subscript 
expression; NMBR1 * 2 + 13 is not a valid subscript expression. 


Rules: 


(1) An array element may not be referenced unless the array has been declared 
previously. 


(2) Except in an EQUIVALENCE statement (where only constants are permitted), 
any and all forms of subscript expressions can be used. 


(3) Except in an EQUIVALENCE statement, the number of subscript expressions 
must correspond to the number of dimensions in the array declarator. 


(4) The value of a subscript expression must be greater than zero. 
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& Examples: 


(1) The following sequence shows how the elements of one array, LIST1, can be 
added to the corresponding elements of a second array, LIST2, to compute 
the values of the array elements in a third array, LIST3. 


"C* FOR COMMENT 


1S huMBer srortean STATEMENT——_________» 
éi7 10 20 


5 30 40 





| fet mei 5.1.9.8 Li.s.Tj 1 ¢,10,) L.t,$,T)2,(,1,0,) ,,,6,1,8,T13,(, 1,0,) 


eeeen a 

eres ieee 

eee Taare 

eerie = hI SIT. 1 OK +, ,L 1S, 7,2, (0, K 
eee K 


K 
| F K ...uT,. 11); ,6,0, ,T,0, ,1,0 
ed lie 
Foes Fae | 
Perna a 


(2) The following sequence shows how only the even-numbered elements of 
LIST1 can be added to the corresponding even-numbered elements of LIST2 
to form another array, LIST4. 


" 
+ 
— 


Geese t? EIN: Sit ON ol S$, TI, (, 1,0 L.t,$,72,(,1,9)) Lit s,T14,(,5,) 


ra) 
u 
~ 


lust os. ty4, (0K) =, ,b,1,8'T, 1,(,2, ,*, ,K) +, ,b,1,8,7,2,(,2, {*, ,«) 


va) 
a“ 
a 
+ 
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(3) The following sequence shows how, by means of a nested loop, the elements 


of a two-dimensional array can be printed out, one element per line. The 


operator .LT. in the logical IF statements means ‘‘less than.’’ 


*C* FOR COMMENT 


ae r 
NUMBER dirorre 
51617 10 


eerne Sear 


The elements will be printed in the following order: 


N.S! 


= 1 
1 


NARRAY (1,1) 
NARRAY (2,1) 
NARRAY (3,1) 
NARRAY (1,2) 
NARRAY (2,2) 
NARRAY (3,2) 
NARRAY (1,3) 
NARRAY (2,3) 
NARRAY (3,3) 
NARRAY (1,4) 
NARRAY (2,4) 
NARRAY (3,4) 


O.N 


0 


NAR RIALY 


O.RiMIALT 1,5,X%.,1!,210 
=, IM, i+ 1 
-L.T,., 4), 160 
=, ({K, .+ 1 | 
(|K, ,-,,T,. 5), {G90 


AN STATEMENT ———__—________» 


3 


T,0 


TO 


4) 


1,0 


3 0 


30 


TE, (,3,,,5,9,) NARRAY (CM, K) 
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Location of Elements Within Array 


Table 2-6 indicates how to find the relative location of an array element within 
an array. 


IF THE ARRAY THEN THE ARRAY REFERS TO THE Nth ELEMENT 
1S DECLARED ELEMENT WHERE N 
AS REFERENCE 1s 


ARRAY (A) ARRAY (a) a 
ARRAY (A,8) ARRAY (a,b) at+A-+(b— 1) 
ARRAY (A,B,C) ARRAY (a,b,c) atA-(b—1)+A-B (ce —1) 





NOTE: A,B,C are the integer values of the array deciarator; a,b,c are the integer values of 
the. subscript expressions in the array element reference. 


Table 2-6. Array Element Location in Array 


Examples: 


(1) If an array is declared as LIST1 (15), the array element reference LIST1 
(9) refers to the ninth element in the array LIST1. 


(2) If an array is declared as NARRAY (3,4), the array element reference 
NARRAY (2,3) refers to the Nth element where N is 2 + 3(3—1) or the eighth 
element. (Check with example (3) in 2.7.4.2.) 


(3) If an array is declared as INPUT (3,4,5), the array element reference INPUT 
(3,4,5) refers to the 60th element, since it is the last array element. By 
checking Table 2-6, it is found that N = 3 + 3(3) + (3) (4) (4) or 60. 


(4) If an array is declared as ARRAY (2,3,2), the array element reference 
ARRAY (1,3,2) refers to the 11th element. However, the array element refer- 
ence ARRAY (3,2,2) also refers to the 11th element. Using the expression 
in Table 2—6 fora three-dimensional array, the relative location within the 
array for the element is evaluated as: 3+2(1)+6(1) or 11. : ~<« 
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3. FORTRAN EXPRESSIONS 


GENERAL 


An expression is a group of one or more elements and operators which, at each 
execution, is evaluated as a single value. FORTRAN expressions are: arithmetic 
expressions, relational expressions, and logical expressions. 


Evaluation of expressions is governed by the priority of operators in the expression. 
Parentheses can be used to force the order of evaluation, regardless of the operators; 
innermost groups within parentheses are evaluated first. 


ARITHMETIC EXPRESSIONS 


An arithmetic expression is a group of one or more arithmetic elements and operators 
which, at execution time, is evaluated as a single arithmetic value. 


The rules for the formation of arithmetic expressions are much the same as those used 
in algebra, except that operations are restricted to exponentiation, multiplication, 
division, addition, and subtraction operations, each indicated by the proper FORTRAN 
operator. With these fundamental operations, more elaborate mathematical functions 

can be built up, such as trigonometric functions and definite integration. Some of these 
functions are supplied as basic external functions (see Section 8); others may be written 
by the programmer. 


Arithmetic Operators 


Table 3—1 is a list of the FORTRAN arithmetic operators and their meanings 
Associated with each operator is a priority number that determines its order of 
evaluation within an arithmetic expression. 


EXAMPLE : 
OPERATOR FUNCTION PRIORITY 
EXPRESSION WRITTEN AS 


exponentiation 
multiplication 
division 


unary minus 
(zero minus) 


unary plus 
(zero plus} 


addition 


subtraction 





Table 3~]. Arithmetic Operators 
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Formation of Arithmetic Expressions 


An arithmetic expression is formed in much the same fashion as in algebra (except 
for the FORTRAN operators and format requirements for array element references and 
function references) using (as primary operands) constants, variables, array element 
references, and function references. Blank characters anywhere in the expression 
have no significance and are ignored during execution; they may be used freely by 
the programmer. Parentheses also may be used to indicate grouping of operands and 
operators. No two-arithmetic operators may appear in succession (two asterisks, even 
if separated by blank characters, are interpreted as the exponentiation operator). 

The unary plus, which is redundant, and unary minus operators must be preceded by 
a left parenthesis except when either one is the first (leftmost) nonblank character 
of the arithmetic expression. 


A more rigorous description of an arithmetic expression with formal FORTRAN 
nomenclature is shown in Figure 3—1 and accompanying text. Note that this figure 
implies the order of evaluation: primary, factor, term, signed term, simple arithmetic 


expression, arithmetic expression. 


atithmetic expression 


signed term+tsimple signed term-simpie 
: arithmetic expression aithmetic expression 


wee) Cae] Cad 


arithmetic 
expression 


ewes simple arithmetic simple arithmetic 


exoressiontsimple expression-simpie 
aithwetic expression arithmetic expression 


ee ee ee ee 


| prieary _| Saimarye* prio: 
[_cvstane | | veiante _| tasetien artnet 
reference expression 
in 


parentheses 





Figure 3}. Structure of Arithmetic Expression 
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An arithmetic primary is an arithmetic constant, an arithmetic variable, an arithmetic 
@ array element reference, an arithmetic function reference, or an arithmetic expression 
in parentheses. 


Examples: 


1 2 3.4 4.5E-2 5.6D3 (S4E-2, 3.0) constants 

ALPHA BETA LIST1 NAME variables 

SUM (IOTA, 2*NMBR-l1, 5) array element reference 
SQRT (A**2 + B*¥*2) function reference 


An arithmetic factor is an arithmetic primary or a construction of the form: 
primary **primary. 


Examples: 
1 4.SE-2. K2**2  K3**NAME  (3E-2,0.04) * CMPLX 


An arithmetic term is an arithmetic factor or a construction of one of the forms: 
term/factor or term*term. 


Examples: 
1 4.5D-2  (3E-2, 0.04) K3**NAME SORT (4.5), 
& . K3**NAME/34 K3**NAME/34/LIST, 
5*6 3.2*SQRT(4.5) K3**NAME/34*2. ALPHA*BETA/GAMMA 


A simple arithmetic expression is a term or two simple arithmetic expressions 
separated by a+ ora ~. 


Examples: 

1 4.5D2 32.*SQRT(4.5) 

3E-2-ALPHA CMPLX + (3.2, 4E-1) 

3E-2-ALPHAtBETA 

A signed term is an arithmetic term preceded by - or +. 

Examples: 

+2.4 -2.4D2 +SQRT(A-B) {2.2,3.3) -ARRAY(1,2,3) 

An arithmetic expression is a simple arithmetic expression, a signed term, ora 


signed term followed by a + or - immediately followed by a simple arithmetic 
expression. 





-UP-7536 
Rev. 1 





3.2.3, 














FUNDAMENTALS OF FORTRAN 3 


SECTION: 


Example: 

( -B - SQRT(B**2 - 4.0 * A * C)) / (2.0 * A) which is one of the roots of a quadratic 
equation 

Type Rules for Arithmetic Expressions 


The data type of an arithmetic expression involving the exponentiation depends upon 
the data type of its operands, as shown in Table 3-2. 


primary 1**primary2 


PRIMARY 1 PRIMARY 2 
u u 


INTEGER integer 
















REAL real 










O.P. d.p. 





COMPLEX complex 


NOTES: O.P. means double precision. 
The letter u means that the result is not defined in standard FORTRAN and it depends 
upon processor implementation of exponentiation. 
\f primary I has a negative value, primary 2 must not be real or double precision. 


Table 3—2. Type Rules for Exponentiation 


The data type of an arithmetic expression involving an arithmetic operator other 
than the exponentiation operator or the unary operators depends upon the data type 
of its operands, as shown in Table 3—3. The data type of a unary operation is the - 
Same as its operand. 


PRIMARY 1 


INTEGER 


primaryl op primary2 
PRIMARY2 ; 


INTEGER REAL aa COMPLE 


u 

















integer 


REAL complex 





u 





O.P. d.p. 











COMPLEX complex complex 


NOTES: opis a nonunaty operator: +, -, *, .. 
D.P. means double precision. 
The letter u means that the result is not defined in standard FORTRAN and depends 
upon processor implementation of the operation. 


Table 3-3. Type Rules for Conventional Arithmetic 
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© Of special importance in Table 3-3 is integer division. Only the integer portion of 
the quotient is retained; the remainder is dropped without roundoff. Thus, the 
arithmetic expressions 0/4, 1/4, 2/4, and 3/4 are evaluated as the integer 0. 


3.2.4, Evaluation of Arithmetic Expressions 
Rules: 


(1) If the value of an arithmetic expression is not arithmetically defined, it cannot 
be evaluated. For example, the following arithmetic expressions must not be 
used: 0**0, X/0.0, and 0**(-3). 


(2) In general, arithmetic expressions are evaluated from left to right governed by 

| ' the priority of the operator (shown in Table 3—1). In the course of this evaluation, 
expressions in parentheses are evaluated before proceeding to the next evalu- 
ation, with innermost parenthetical expressions evaluated first. 


There is a permissible exception. If mathematical use of operators is associative, 
commutative, or both, the order of evaluation may be changed internally to take 
advantage of these qualities, provided that integrity of parenthesized expressions 
is‘not violated. (An operation is associative if A op B op C can be evaluated as 
A op (B opC) orA op B) op C with no algebraic change in the results; an operation 
is commutative if A op B can be evaluated as B op A with no algebraic change 
in the result.) The only associative FORTRAN arithmetic operators are + and *; 
the only commutative FORTRAN arithmetic operators are + and *. Thus, in some 
processors the expression A + B + C can sometimes be evaluated as A + (B +C) 
@ and, at other times (in the same processor) as (A t+ B) + C; the expression A*B 
can sometimes be evaluated as if written B*A. The associative and commutative 
laws do not apply to evaluation of integer terms containing division; evaluation of 
such terms proceeds from left to right. For instance: K*M/N is evaluated as (K*M)/N. 





(3) The evaluation of any function references in the expression must not alter the 
value of any other element within the expression or statement that contains the 
function reference (see Section 8). 


Examples: 


(1) Evaluate (-B-SQRT(B**2.0-4.0*A*C))/(2.0*A) where A, B, and C have the values 
of 1.0, -3.0, and -10.0 respectively. The order of evaluation is: 


(-B-SQRT(B**2.0-4.0*A*C))/(2.0*A) 
(a) (-B-SQRT(9.0-4.0*A*C))/(2.0*A) 
(b) (-B-SQRT(9.0-4.0*C))/(2.0*A) 
(c) (-B-SQRT(9.0-(-40.0)))/(2.0*A) 
(d) (-B-SQRT(49.0))/(2.0*A) S 
(e) (-(-3.0)-SQRT(49,0))/(2,0*A) < 
(£) (3.0-SQRT(49.0))/(2.0*A) 
(g) (3.0-7.0)/(2.0*A) 
(h) (-4.0)/(2.0*A) 
(i) -4.0/(2.0*A) 


G) -4.0/(2.0) 
(k) -4.0/2.0 
(1) -2.0 
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Note that in these steps, the successive evaluations are shown in terms of exact 
values. In actual practice, these values will be approximated. The degree of 
approximation will depend upon the processor implementation of real type 
arithmetic. 





(2) Evaluate the expression (N/2)*2-N. 


All possible results are shown in the following truth table. (Remember that, in 
integer division, the remainder is truncated.) 


POSITIVE/NEGATIVE ODD/EVEN RESULT 
odd 


positive -1 
-positive even 
zero 


negative 


negative 





Note that this expression can be used for odd/even testing of an integer value. 


(3) Intrinsic functions (see 8.3) or assignment statements can be used to get around 

some. of the restrictions imposed by the mixed type requirements of Table 3-3. 
It must be understood that some precision may be lost in this process and that 
it cannot always be done, since the range of a real or double precision value 

is always greater than the range of an integer value. (The restrictions on the 
use of an assignment statement for data type conversion are described in 4.2.) 
The following is an example of the need for such conversion in the calculation 
of the volume of a room where the length is originally given as an integer. 





Cl NEXT) LINE USE S| ION TRLIUN Sit Ci FUN CIT! ON, JF LIOAT 

V0.0 UMIE =F LiOlA TiC Lh EIN, G:TiH.) 1* Wil O.TIH,* HE, 11G.Hi TF, is | 
ci wwe.x(r| LIN ELS H,OlW 0 A Ni M N A Mi EiN 
Pe a ape | X LN, G T=, L, E: NIG: TH ‘ 


0 MiE = xX LIN D » Waee. » Bi * Gu { i n 
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A relational expression defines a relation between two arithmetic expressions. At 
execution time this relation is evaluated as either true or false. 


3.3.1. Relational Operators 


A relational expression is made up of two arithmetic expressions separated by one 


of the relational operators shown in Table 3—4. Blank characters may be used freely 
to improve readability. 


OPERATOR MEANING 


greater than © 


greater than or equal to 
less than 
fess than or equal to 


equal to 





not equa! to 


NOTE: The periods, as shown, are necessary parts of relationat 
operators. 


@ Table 3-4. Relational Operators 
Examples: 


If KO, K1, K2, KK1, and KK2 represent the values 0, 1, 2, -1, and -2, respectively, 
then the expression: 


K1L.EQ.1 is true. 

K2.GT.KK2 is true. 
KK1.NE.KO is true. 
KK1.LE.KK2 is false. 
K1**K1.NE.K2**K0O is false. 
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3.3.2. Type Rules for Relational Expressions 


Only the combinations indicated in Table 3—5 are permitted for the relational 
expression: 


expl op exp2 
where: opis any relational operator, and exp! and exp2 are arithmetic expressions. 


expl op exp2 


INTEGER no no 
REAL 

O.P. 

COMPLEX 


NOTE: O.P. means double precision. 








Table 3~5. Type Rules for Relational Expressions 


Where a real arithmetic expression appears with a double precision arithmetic 
expression, the relation is evaluated as if written as (expl-exp2) op ODO. 


3.3.3. Applications of Relational Expressions 


(1) In this sequence a set of data cards is read, each of which contains an integer 
value right-justified in the first 10 columns of the card. The program is to find 
the greatest value and print it out. It is assumed that none of these values will 
be - 999999999, so that a card containing - 999999999 indicates the end of the 
set containing at least one integer. 


feeee a VIN, AO ..1[GT. NIM BR NM BR [= NUFIO: : 
i | leo, 7,0) 2,0 


WARE TIE1C 35, 14:01) IN iM, BIR ‘ 
F,ORMAIT Ch tat | \ 


$:T,0,P : \ , 


t 
Ny DD: aaa: | H A l | ie ae a Co | pot 
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Note that the contents of the first data card must be read into NMBR to create 


a basis for comparison. Then the value of the next card is read into INFO. The 
fourth line is a logical IF statement with the relational expression INFO.EQ. 
—999999999. If this relation is true, the GO TO 30 statement causes a jump 
out of the loop and a printout of the contents of NMBR. If this relational ex- 
pression is false, GO TO 30 is disregarded and the next line, which is 
another logical IF statement, is executed. This statement tests the value of 
INFO. If the value of INFO is greater than the value in NMBR, the value of 
INFO is transferred to NMBR. Then control returns to 20 and a new value for 
INFO is read. Cards are read until a card containing - 999999999 is reached. 
Note that it is necessary to test for -999999999 before the attempt to update 
NMBR; otherwise, the NMBR printed out will be -999999999. 


(2) This sequence compares two complex numbers, KMPLX2 and KMPLX4. If their 
real portions are equal and the imaginary part of KMPLX2 is greater than the 
imaginary part of KMPLX4, control is passed to one set of statements; otherwise, 
control is passed to another set of statements (in this case, beginning with the 
statement label 50). 


This sequence introduces two intrinsic functions, REAL and AIMAG. The 
intrinsic function reference REAL (x) obtains the real portion of the complex 
variable x; the intrinsic function reference AIMAG (x) obtains the imaginary 
portion of the complex variable x. 


pea ta Pe Li 


U A a A Ai ML A Friese i it rae fee 
«| j ' | 
e 
. f 1 , : i 
F RIE AL KIM, PL, X,2 -N EL., RE AIL KM, PIL XxX 4 GO; .T|O, 5.0 
F.C AyE MAG CK MPLX2)), LIE AI MIA,G (, KMIP.L.X,4,)/), 6,0, |T,0, 5,0] 
! Cian oan RTO Fe doe ae oe 


Note the numerous parentheses used in the logical IF statements. The format 

of a logical IF requires one set of parentheses around the relational expression; 
the format of a function reference requires a set of parentheses around the list 
of arguments. In writing such statements it is not unusual to occasionally mis- 
place a parenthesis. It is a good idea to check such statements and ensure that 
the number of left parentheses is the sameas the number of right parentheses. 
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3.4. LOGICAL EXPRESSIONS @ 


A logical expression is a group of one or more logical elements and operators which, 
at execution time, is evaluated as either true or false. 


3.4.1. Logical Operators 
The three logical operators are shown in Table 3-6. 
NOT. 1 


logical negation 
-AND. 






logical conjunction 


logical disjunction 


NOTE: Parentheses can force the order of evaluation. 


Table 3—6. Logical Operators 


The meanings of the logical operators are more precisely illustrated by Table 3-7 
where the format of a logical expression is /el.op./e2 for the .AND. and .OR. 
operators, and .NOT. Je3 for the .NOT. operator. Each le is a logical element. 


AND 


le2 lel. AND. Je2 | lel.OR Jez 
ts ts IS 





Table 3~7. Truth Table for Logical Operators 
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3.4.2. Formation of Logical Expressions 


A logical expression is formed with logical elements and logical operators. A 
description of logical elements and their combination to form a logical expression 
is shown in Figure 3—2, and the text which follows the figure. Figure 3-2 also 
implies the order of evaluation: primary, factor, term, logical expression. 


logical 
expression 





logical logical expression.OR.logicai expression 


term 


logical term.AND.logical term 


logical 
factor 


logical -NOT.logical primary 


primary 









logical logical 
constant vatiable 


iogical 
attay 
element 


logical relational logical 
function expression expression 
reference in 
Darentheses 


Figure 3~2. Structure of Logical Expression 


A logical primary is a logical constant, a logical variable, a logical array element, 
a logical function reference, a relational expression, or a logical expression in 
parentheses. 
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Examples: 

-FALSE. logical constant 

ALPHA logical variable 

SWITCH (K2, K4, K6) logical array element reference 
» XXX(2) logical function reference 

A.LE.0.034 relational expression . 


A logical factor is a logical primary or a construction of the form: 
-NOT. logical primary. 


Examples: 

-TRUE. logical primary 
.NOT..TRUE. evaluated as false 
A.NE.B logical primary 
.NOT.A.EQ.B equivalent to A.NE.B 


A logical term is a logical factor or a construction of the form: logical term .AND. 
logical term. 





Examples: 

ALPHA logical factoz 

-TRUE..AND..FALSE always evaluated as false 
-TRUE..AND.ALPHA evaluated as true only if ALPHA is true 
A.GT.B.AND.A.GT.C.AND.D evaluated as true only if A is eects 


than B and greater than C,and D is true 
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A logical expression is a logical term or a construction of the form: 
logical expression .OR. logical expression. 


Examples: 


.TRUE..OR..FALSE. always evaluated as true 


-FALSE..OR.ALPHA evaluated as true only if ALPHA is 


true 


A.GT.B.OR.C evaluated as true if C is true and/or A 
is greater than B 


A.GT.(B.OR.C) an invalid expression because both 
operands of a relational expression must 
be arithmetic expressions 


A.GT.B.OR.A.GT.C evaluated as true if A is greater than C 
and/or A is greater than B 


A.AND..NOT.B.OR..NOT.A.AND.B called the exclusive OR function; evaluated 
as true only if either A or B is true, but 
not both 


(A.OR.B).AND..NOT.(A.AND.B) another way of writing the exclusive OR 
function 


Evaluation of Logical Expressions 


The order of evaluation of logical expressions is determined by the priority of the 
logical operator(s), as shown in Table 3—6, and the order of evaluation implied by 
Figure 3—2, without violating the integrity of elements in parentheses. When 
parentheses are present, innermost parenthetical expressions are evaluated first. 
Thus, parentheses can be introduced to force the order of evaluation. When two 
elements are combined by an operator (or an element follows .NOT.), the element(s) 
must be evaluated before the logical relation can be evaluated. : 


Examples: 


.NOT.A.AND..NOT.B evaluated as true only if both A and B 


are false 
NOT.(A.AND.B) evaluated as true if A and/or B is false 
A.AND.B.OR.C evaluated as true only if both A and B are 


true and/or C is true 


A.AND.(B.OR.C) evaluated as true only if A is true and 
either, or both, B or C is true 
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3.4.4, Applications of Logical Expressions : @ 
An application of logical expressions is shown in the following example. 
In this program, an indicator called IND is set to —1 if integers K2 and K3 are both 


negative, to +1 if both are not negative, or to 0 if one is negative and the other is 
not. 


ND = ’ 
1 iF K,2,. ,G,E,.,0,. ,A,N{D,. K,3,. .G,E,. ,0 1.N,D, = +,1 


IF. ( jK2,. LT,.,0,. AND. K3,.,L,7.,0,))] IND =, -1, | 
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4. ASSIGNMENT STATEMENTS 


GENERAL 


Execution of an assignment statement causes assignment of a value to a variable 
or array.element. This new value becomes its current value until the variable or 
array element is redefined. 


There are three assignment statements in FORTRAN: the arithmetic assignment 
statement, the logical assignment statement, and the GO TO assignment statement. 
Because the GO TO assignment statement can only be used with an assigned GO TO 
statement (which is a control statement), it is discussed in Section 5 in conjunction 
with the assigned GO TO statement. 


ARITHMETIC ASSIGNMENT STATEMENT 
Function: 


The arithmetic assignment statement evaluates an arithmetic expression and assigns 
this value to an arithmetic variable or an arithmetic array element. 


where: vis an arithmetic variable or arithmetic array element; 


eis any arithmetic expression. 





Operation: 

An arithmetic assignment statement is performed in up to three steps: 

(1) The arithmetic expression e is evaluated to yield a single numerical value. 

(2) The data type of this single numerical value is converted to the data type of v. 


(3) This converted single value replaces the contents of v. 







4 
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SECTION: 


With arithmetic assignment statements such as K = 3 or S = 3.0, no evaluation or 
conversion is required. Since the statement can be performed in three distinct steps, 
it is possible for the same variable or array element to appear in both v and e. For 
example, the sequence 


| le AY K,+,2 K .+,3 





assigns the value of 7.0 (its real type approximation) to the array element ARRAY(6). 


The arithmetic expression is evaluated and its value assigned to the variable or array 
element in accordance with the rules in Table 4-1. 


INTEGER DOUBLE COMPLEX 
PRECISION 


INTEGER 


REAL 


DOUBLE 
PRECISION 





COMPLEX 


Assign e to v without change. 

Truncate any fractional part of e and assign result to v as an integer. 
Transform e to real type value and assign this value to v. 

Transform e to double precision value and assign this value to v. 


Evatuate e by rules of 3.2.4 (or any more meaningful rules), transform to type of v, and assign 
to v. 


Prohibited combination. 





Table 4—]. Type Conversion by Arithmetic Assignment Statement 
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Rules: 


(1) Some arithmetic assignment statements should be avoided to save execution time 
or avoid possible inaccuracies. For example, the statements A = 1 and I = 10.0 
should be replaced, if possible, by the statements A = 1.0 and I = 10, respectively. 


(2) Some arithmetic operations can be avoided by using the appropriate form of 
constant, thus saving execution time. For example, the statement X = 2 E04 
is preferred to X = 2.0 * 10.0 ** 4.0. The two values may not be identical due 


to the approximations involved in the computer representation of real type data. 


(3) References to intrinsic functions are available (see Section 8) for handling complex 
type data and converting the real and imaginary parts as required. 


Examples: 


(1) This sequence shows how to compute and store the fractional remainder when the 
integer K2 is divided by the integer K3. 


fii. | lou P\REIC 1 Sit OIN, .K.2, 21, K 


ee 
rat 

uoyu 
Tas A 





(2) This operation can be source coded more conveniently with the basic external 
function DMOD (see Table 8—3) as follows: 


| 
fi. | fou ete PREC St ON, )K 12,21, Ki3 31, IRiEiM tei 
Pileseie er 


\ fae re boae 
Pisa | : si 
eee ce | 
ee ieee 7 roan 
sa | deem 01,0. K12.2., K.313 Cas SW Te 


(3) The next example introduces the intrinsic function FLOAT and shows how it can 
eliminate an arithmetic assignment statement used for type conversion. 


i= iN zg | i 


P| ome: THAW DTH HE GIT La 
ley eesxts! 


>< 


1 
A N Lh A HUN 7 UN 


ee got cha LiU MIE, = FLOAT ( L EING TH) * Wil OD TIH: *: HE UGH T 














4 






UP-7536 
Rev.1 








FUNDAMENTALS OF FORTRAN 


SECTION: 





4.3. LOGICAL ASSIGNMENT STATEMENT 





Function: 


This statement evaluates a logical expression and assigns this value (either true or 
false) to a logical variable or logical array element. 


v=-¢€é 


where: v is a logical variable or logical array element. 





eis any logical expression. 


Execution of this statement consists of two parts: the evaluation of the logical 
expression as either true or false, and the assignment of this logical value to the 
logical variable or array element. 


Examples: 


(1) The following sequence evaluates x as true if K is an integer less than or equal 
to 3. 





(2) With each execution of the following loop, SWITCH will alternate in value from 
true to false and can be used as an odd/even counter within the loop. 


0,GIC AL wiist,Gk 


Pio | swititom = FALSE ! es 
di iol [swurcuy es buowT § wu TcH ! 
Bisel a 

ian 2 ee 

bis. | leo, ror 1.0 

aaaes! 
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5S. CONTROL STATEMENTS 


5.1. GENERAL 


Control statements modify the normal sequence of execution. Some of these statements 
specify unconditional modification of the normal sequence; others contain a test that 
determines whether or not the sequence of execution shall be changed. 


Execution of a program or procedure starts with the first executable statement and 
continues sequentially until a control statement is encountered. 


The control statements of FORTRAN are: 
a GO TO statement 
gm IF statement 
m DO statement 
: = CONTINUE statement 
@ , a STOP statement 
a PAUSE statement 
@ CALL statement 
= RETURN statement 
- The CALL and RETURN statements are not described in this section because they are 
associated with external procedure subprograms. These statements are described in 
Section 8. The GO TO assignment statement, because it is associated with the assigned 
GO TO statement, is described in this section. 
5.2. GO TO STATEMENTS 
The GO TO statements are: 
@ unconditional GO TO statement 
@ computed GO TO statement 


m assigned GO TO statement 
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5.2.1. Unconditional GO TO Statement 
Function: 


To transfer control, unconditionally, to a statement specified by statement label. 


GO TO s/ 


where: s/ is the statement label of an executable statement within the same 
program unit. 





Rules: 


(1) si must be the statement label of an executable statement. 

(2) Any executable statement immediately following the unconditional GO TQ 
statement in the program unit must have a statement label, otherwise it can 
never be executed. 


Example: 


After the following sequence is executed, the variable K will have a value of 8. 


5.2.2. Computed GO TO Statement 
Function: 


To transfer control to one of several listed statement labels, as determined by a 
previously defined integer value. 


GO TO (sly ,8lp,...,8l,)) 1 


where: each s/ is a statement label of an executable statement in the same 


program unit, separated from the next s/ in the list by a comma. 


tis an integer variable representing a value, such that l<i<na. 
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Rules: 


(1) If the value of ris 1, the statement with sl, will be executed; if the value 
is 2, the statement with sl, will be executed; and so on. 


(2) The value of 7 must be defined before execution of the statement. 


(3) The parentheses around the list of statement labels and the comma before the 
integer variable are required as shown in the format. 


(4) There is no restriction on other uses for the integer variable 1. 


(5) There is no standard FORTRAN restriction on the maximum value of i ora, 
but it is possible that a particular processor may specify a maximun value. 


Examples: 


(1) After execution of this sequence, control is transferred to the statement with 
statement label 35. 





rere ° ee ne GEO OBC EC APT OO a aC oe 
Pipe td : bet a 
i zoel ft TE MAK 24,1 phe og peg leg op eds fy a pgp 


rdskeaicuale 0: 17,0 a PE VMbe st ha tot ee ot lets ep Cl eps ob 
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(2) Another application of the computed GO TO statement shows how it can be 2 
used to determine the course of processing in the main program. A data card & 
having a decimal digit punched in the first character position is read in. If 
this digit is 1 through 4, it indicates a processing sequence. Any other digit 
(or the blank character) is treated as an end of file indication. 





NTE GIER, CHO, CE 


| 


H1011 E Q:. 10 OLR H:0.!}'C E.. Tl. 4 $|T,0.P 


iG) 
© 
e 


a 
fe 
+4 





7) 
°o 
o~{ 
() 


m 
4 
o 


MS Pee Pe ah ee ae ee 
Pr) ar oe sy . 
4 
° 
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(3) The following example shows how the computed GO TO statement can be used 
to create a closed internal block (a group of statements that cannot be entered 
by normal sequential execution, but must be called for execution by control 
statements). 


Hoes SET OD RRO IOE TES © aN 
INTERNAL 
ener uiemeonm deal 
Qo, . TIO 10,0, ,.2,0,01,,3,0,0,,,4,6,0 K 
Remar He 


Pde eh CLS TIO, 12,0 
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5.2.3. Assigned GO TO Statement 





Function: 


To transfer control te one of several listed statement labels, as determined by an 
integer value previously defined by a GO TO assignment statement. 


GO TO i(sly, Sl, oe -) 


where: each s/is a statement label of an executable statement in the same program 
unit, separated from the next s/ in the list by a comma. 


ris an integer variable, followed by a comma, previously assigned a value 
by a GO TO assignment statement. This value is equal to an s/ in the list. 








(1) Prior to execution of the assigned GO TO statement, the integer variable i must 
have been assigned a value by a GO TO assignment statement. 


(2) Although standard FORTRAN places no restriction on the number of statement 
labels in the statement (except that there must be at least one), a particular 
processor may specify a limit. 


(3) Standard FORTRAN specifies that the value assigned to 7 must be an s/ in the 
list. If the value of ; is not an s/ in the list, some processors will treat this 
condition as if the value of i were in the list. The manual for the processor will 
specify how this condifion is handled. 





Example: 


Logically, the assigned GO TO statement can be used whenever a computed GO TO 
statement is used (see 5.2.2). The format requirements differ and the assigned GO TO 
statement requires at least one previous ASSIGN statement (GO TO assignment 
statement). Figure 5—1 shows how the assigned GO TO statement can create a multi- 
legged GO TO after a series of statements shared by different parts of the same pro- 
gram unit. 
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ASSIGN 7 TO K2 








Figure 5=1. 


GO TO K2, (18, 20, 7, 42) 


20 WRITE... 
e 
e 





18 C(1, 2,3) 5... 


Use of Assigned GO TO Statement 
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5.2.3.1. GO TO Assignment Statement 





Function: 


To assign a statement label to the integer variable that is used in an assigned 
GO TO statement. 


ASSIGN k TO i 


where: kis the statement label of an executable statement in the same 
program unit. 


ris an integer variable used in an assigned GO TO statement of the 
same program unit. 





Rules: 


(1) The ASSIGN statement must be executed before the assigned GO TO statement 
to which it applies. 


(2) The statement ASSIGN 10 TO J does not have the same meaning as the 
arithmetic statement J = 10. In particular, the sequence 





will not produce a meaningful result. If J is to be used as an arithmetic 
variable after its use in an assigned GO TO statement, it must be defined 
as such at some point after the assigned GO TO statement. The following 
sequence illustrates proper use of the ASSIGN statement: 


AS STIGN, 10 TO J 





_ 
“ 
— 
wn 


~K- 
it) 


Control is transferred to statement 10 after execution of statement J. The 
variable K is assigned a value of 14. 


Example: 


See Figure 5—1. 
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5.3. IF STATEMENTS 


The IF statements are the decision-making elements of FORTRAN. The test specified 
in the IF statement may modify the normal sequence of execution. FORTRAN provides 
two IF statements: the arithmetic IF statement and the logical IF statement. 

5.3.1. Arithmetic IF Statement 


Function: 


To act as a three-way branch, as determined by evaluation of an arithmetic expression. | 


IF (exp) sl, ,sly,Sl, 


where: exp is any arithmetic expression except complex type. 


each s/is a statement label of an executable statement in the same 
program unit. 





Rules: 


(1) If exp is negative, control is transferred to sly; if zero, to Sly; if positive, to 
sl. 
~ (2) The list must contain three statement labels; however, any two, or all three, 
may be the same. If all three are the same, the statement is, in effect, an 


@ unconditional GO TO. 


(3) An arithmetic IF statement must not branch to itself. For example, the following 
statement is illegal: 





- (4) The first executable statement following an arithmetic IF statement must have a 
statement label or it can never be executed. 


(S) If the arithmetic expression is real or double precision type and contains trunca- 
tions and roundoffs, caution must be exercised, especially when the expression 
is tested for the zero condition. 
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Examples: 
(1) If, in solving the quadratic equation ax” + bx + c, where all three coefficients @ 
are integer type numbers, the discriminant b“ — 4ac is negative, the result is 


two complex roots which are conjugates of each other; if zero, two equal real 
roots; if positive, two unequal real roots. The following sequence shows how an 
arithmetic IF statement can be used to select one of three root evaluation 
procedures based on the evaluation of the discriminant. 


NDIS,;C = NBt*2 j= 4 * NA * NEC 


tte cN py sc) 10, 20, 30 


pir oce dure for two complex roots 


|. 20] | pir 9 cedure for two equal real roots 


eens oir .o¢ ed vr ie fio.r, t we vurneiqucal real foots 


(2) Now consider another procedure for the problem in (1), but with the coefficients 
real type numbers. This sequence introduces the basic function reference ABS(x} 
which returns the absolute value of a real type argument, x. 





O1rsc,=B** 224.07 A*C 


|e AtB.$( 0,1 ,8.C,) - 1j&-.10 40 5j0., 5,0 
rr 


Ol! S$Cy=,0,. 0 } 


| sol tr on sc) 10 ,.20,, 30 
Pe OT two, ¢,omplje.x, 0,0¢ § 


two equa} real ef oaty;s | 


meee two, vine qujat, ..¢ail rojo ts ! ; ; . i 
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In this program, all absolute values of the discriminant less than 10719 ate 
treated as zero. This handles possible errors occurring in the value of the 
discriminant due to computer representation of real type values, truncation, 
and roundoff. The criteria for these limits vary with the nature of the problem 


and the manner in which a particular computer treats real type values and 
computations. 


(3) Paragraph 3.2.4, example (2), contains an arithmetic expression that could be 
- used for testing odd/even integers. The following example shows how this 
expression might be used. 


1 oF Cyt N72.) * .2.-N) GO 790 100),,200, 100 


10,0] | re ofustji nie, for odd, jt, 9,t ,e,9,e/F ,s 

neeee 

| 200] | rf ojwtine for even integ rs 

ae : 

aes : 


Logical IF Statement 


Function: 


To determine whether or not a single executable statement, written as part of the 
logical IF statement, shall be executed. 


IF (e) s 


where: eis any logical expression. 


s is any executable statement except a DO statement or another logical 
IF statement. 





Operation: 


Execution of the IF statement proceeds in two parts: evaluation of the logical ' 
expression and (possibly) execution of the statement. If the logical expression, e, 
is evaluated as true, the statement, s, is executed; if false, the statement s is 
ignored and control is passed in normal sequence to the next executable statement. 
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Rules: & 


(1) A logical IF statement cannot refer to itself. For example, the following statement 
is illegal. 





(2) Because execution of the statement proceeds in two parts, the same variable 
may appear in both the logical expression e and the statement s. For example: 





(3) Caution must. be exercised if the logical expression involves the comparison 
of real or double precision arithmetic expressions. 


Examples: 


(1) The following program determines how many months it takes for a deposit to 
double its original value for three different compound interest rates. The rate 
is the interest computed at the end of each three-month period, starting from 
the month of deposit. The three interest rates are 1%, 1.25%, and 1.5%. The 
months required for the three different rates are stored in the array MONTHS. @ 





DIMENSION MONTHS (3) 


ieee BS 
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(2) This program finds the values of y for y=3x 72x45 for values of x between 0.1 
and 0.2 at intervals of 0.001, that is, x = 0.100, 0.101, 0.102, ..., 0.199, 
0.200. This gives 101 values of y, stored in array Y. 


Di MENS! ON Y( 10,1) { 
eee eR: 
to] [ue net 
P| beer N 


To | [xe tyoat.c.%) 7183 


Ree eee ee ee 

PL hemp. tat. 10,1), 60, 7.0, 10 

ea BR | 
P| [eno 

eee PeraEs 


Because the integer type variable N is used for the counter and in the logical 

IF statement, the result will be exactly 101 values. Because a division operation 
is used to compute X each time around the loop, each sample of X will be as 

close to the desired value as the processor approximation to a real number 

permits. This avoids any pitfalls that might arise if X had been incremented 

each time around the loop by a real type value, by avoiding possible approximations 
due to successive truncations and roundoffs. 


DO STATEMENT 
Function: 


To initiate and control repeated execution of a set of executable statements. 


DO na FM, M7,M3 


or 


DO em re ) 


where: nis the label of an executable statement called the terminal statement, 
which follows (not necessarily immediately) the DO statement. 


iis an integer variable called the control variable. 
my is an integer constant or an integer variable called the initial parameter. 
m, is ‘an integer constant or an integer variable called the terminal parameter. 


m. is an integer constant or an integer variable called the incrementation 
parameter. In the second form, its value is implicitly 1. 
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Operation: 


For successive repetition of the same group of executable statements, the DO state- 
ment eliminates separate statements that set a variable to a starting value, increments 
this value after the group has been executed, and tests the new value to determine 
whether the group shall be executed again, as was done in the sample program of 

1.6. The operation of a DO statement is shown in the following examples. 


WITHOUT DO STATEMENT WITH DO STATEMENT 


=m, DO n Fm) ,M7,M3 


label first statement of set ( first statement of set 
DO 


range : 
n terminal statement of set 


n terminal statement of set 
=m 
IF (.LE.m,) GO TO label 





The steps in the execution of a DO statement are: 


(1) The control variable i is assigned a value represented by my. This value must 
be less than or equal to the value represented by My. 


(2) The range of the DO is executed. If this range contains a reference to a function 
or subroutine, the function or subroutine is considered part of that range when it 
is executed. 


(3) After execution of the terminal statement, the control variable is incremented by 
the value of m,. 

(4) This new value of the control variable is tested. If it is less than or equal to 
the value of m,, program control is transferred back to the first statement of the 
DO range, with the new value for the control variable; if it is greater than mM, 
the DO is satisfied, and the control variable becomes undefined. 


(5) If this DO range is nested within another DO range, and both have the same 
terminal statement, when the inner DO is satisfied, the control variable of the 
next outer DO is incremented and tested, and its DO range (which includes the 
inner DO) will be repeated until satisfied. This will continue for all nested DO 
statements sharing the same terminal statement until the outermost DO is satisfied. 
If there is no nesting of DO’s with the same terminal statement, after a DO 
statement is satisfied, its control variable becomes undefined and program control 
is passed to the first executable statement after the terminal statement. 
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Rules: 


(1) The terminal statement must be in the same program unit as the DO statement. 
It must not be a GO TO, arithmetic IF, RETURN, STOP, PAUSE, or a logical 
IF containing any of these forms. However, if the logic of a DO range indicates 
that such a statement is a terminal statement, such a statement can be followed 
by a CONTINUE statement (which has no logical function); the CONTINUE 
statement is then labeled and used as the terminal statement of the DO range. 
If the terminal statement is a CALL statement, the subroutine will be executed; 
after the RETURN of the subroutine is executed, the control variable is tested 
to determine whether the DO range shall be executed again. 


(2) At execution time the parameters of the DO statement must be defined as values 
greater than zero. 


(3) Because the contro! variable is tested at the end of the DO range execution, a 
DO statement will always be executed at least once when encountered. 


(4) No statement in a DO range may redefine the control variable or any parameter 
of the DO statement; however, the control variable may be referenced in the DO 


range, as in: 


DO na iFM ,m5,M 3 


K=r INT 


n terminal statement 


(5) If a control statement causes an exit from a DO range before the DO is satisfied, 
the control variable remains defined until redefined. 


(6) DO statements can be nested in outer DO statements with this restriction: the 
range of each nested DO must be completely contained in the range of its next 
outer DO and may share the same terminal statement. 
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VALID NEST INVALID NEST 
DO10 =... (The range of the second DO is not 
A contained within the range of the first 
DO.) 
DO 202... DO10 =... 
DO 20 =... DO 20 =... 
20 executable statement 10 executable statement 


. . 


10 executable statement 20 executable ‘statement 





=p A special type of nest is a completely nested nest. This is a nest of DO 
statements which satisfies both of the following conditions: 


® The outermost DO statement of the nest is not contained in the range of 
another DO statement. 


m The first occurring terminal statement within the nest physically follows 
(not necessarily immediately) the last DO statement within the nest. 


COMPLETELY NESTED NEST VALID NEST BUT NOT COMPLETELY 
NESTED NEST 


(The first occurring terminal statement 
does not follow the last DO.) 


sf es DO 10... 


DO 20... 


20 executable statement 
DO 30... 


20 executable statement 30 executable statement 


10 executable statement 10 executable statement 
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A DO range can have an extended range if it is a nonnested DO or the innermost 
nest of a completely nested nest and both of the following conditions are true: 


@ The sange contains a GO TO, or arithmetic IF, or a logical IF containing either 


of these two statements, that can pass control out of the range to another set 


of executable statements. 


@ This set of executable statements contains a contro! statement that could return 


control back to the DO range. 


Examples of extended ranges are: 


DO 10... 


* 


IF (. . .)20,20,30 


executable statement 


e 


executable statement 


terminal statement 


. 


GO TO 50 


executable statement 


extended 
range 


IF (...)GO TO(40,.. 
executable statement 





The extended range of a DO must not contain, in the 


statement that has an extended range. 


DO 10... 


IF (.. .)}GO TO 20 


executable statement 


terminal statement 


. 


executable statement 


extended 
range 


GO TO (30, 40),-K 
executable statement 


same program unit, a DO 
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(8) A GO TO, arithmetic IF, or a logical IF with either of these two forms cannot 
pass control into the range of a DO unless that GO TO or arithmetic IF is being & 
executed as part of the extended range of that DO. If more than one DO statement 
share the same terminal statement, only the innermost DO range may contain a 
GO TO or arithmetic IF that transfers control directly to the terminal statement. 

The following shows examples of legal and illegal transfers of control: 








VALID TRANSFERS OF CONTROL INVALID TRANSFERS OF CONTROL 


| z 





DO 






t 


extended 
range 








Examples: 


(1) The following sequence shows how to set all array elements of a one-dimensional 
array to zero. 


0,1 M EIN,S,!,0.N, ,T,ABLIE 1,0,0 


fed td ; | 77 
heictents dell (Be 110, Ki=.3,, 1, 0,0 
Eat 


TAB LIE(,K,) =,0,.,0 
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(2) The following sequence shows how to set all array elements of a three-dimensional 


array to zero. 


Po | onreiceR 
we 
OSs 2 


eee eee 
a 10, .K 
Perens 
ae alial SET. (iby, Kk 


SE T,(,3,,14,,,5,) 


(3) The following example is an illustration of processing by DO statements. 


ees a 
Peaee 5 


D,O, 110 K 
LIN T.GIR.S 


etenl 

rest 

L ataiee 
D.9 0 K 


rete 
f. 4.2.0] [in reie.s / 
eee 
edit 


A 


NiSi/ O.N 


LNT. GIRIS 15 


K,*K 


IN, TGIR,S, (OK, 2 


After the first DO (with terminal statement 10) has been satisfied, the array 
elements of INTGRS will contain, in succession, square numbers, l, 4, 9, ..., 
225. After the second DO has been satisfied, all odd-numbered array elements 
will contain the value 1; the even-numbered array elements are unchanged. 


(4) 


The following internal numerical sort program contains an extended range. 


Originally the array elements of array NMBRS contain a series of integers in 
random sequence. The program is to re-arrange the integers from !ow to high, 
with the lowest in the first array element. 
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The sorting method used is known as a ‘‘bubble’’ sort. The first array element r 
is compared with the second array element; if the first is greater than the second, & 
the two are interchanged. Then the second array element is compared with the 

third array element, and so on, until the next-to-last array element is compared 

with the last array element. This is the initial pass. The result is that the 

greatest number has ‘‘bubbled’’ through to the last position. A record is kept 

of the last positions interchanged. This record determines if another sorting pass 

is required and how many array elements must be compared in that next sort. 


Di! MEIN S| ON, ON MBRIS(,1,0,0)) 


po | four 





=,1 ' 

oe He 
|. 4,0| (0,0, 210, Mi=,1,,. 
se Bl 1. F,(,.NIM,B.R,$,(.M,),-,G, TI. .NM,B.R.S, (.M,+, 1 G0, .T,0, , 3/0 
Sener. ; 
ica scan J.=,K,O|U NT -.1 
Pd decay ee, ,0 $,T,0}P 
Pt | leo, te, 4,0 
P30 [Kou NiT.=.M 

N.S A,VjE,=.N.M,B_RS.(.M 


| | ntsMBR/S,( .M,) ,=,N.MjBLR,S|( (Mis .1,) 
| INM BORIS. (.M+.1,) 2.N,S ALIVE 


ee aPC er 
eae 


END 


Note the following: 


w The terminal statement cannot be a logical IF containing a GO TO; therefore, 
a CONTINUE statement is added and used as the terminal statement. The 
CONTINUE statement (also known as a ‘no-op’’) does not perform any logical 
operation and can be used to satisfy format rules of FORTRAN. 


mw The extended range is considered part of the DO range being executed; therefore, 
parameters of the DO range can be referenced in the extended range. 
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(5) A DO statement can be used to advantage in cases where its use is not readily 
apparent. One such case is the evaluation of polynomials which have the form: 


n n~1 
a,x + agx Faoeet ayx + ane] 
where each a and the value of x are known as execution time. 


For polynomials where n is no greater than 3 or 4, an arithmetic assignment 
statement can be used, but as the value of a increases, execution time can be 
reduced with a DO statement. 


Consider ayx3 + ay x2 + agx + ag. By successive factoring of x, this can be 
represented as x(x(x(ay)tao)+a3)ta4 which consists of three ‘‘multiply by x 
and add a constant’’ operations. With the DO statement, this can be evaluated 
by: 





For the general case, where n has the value N: 


acta vet a 

heehee Idle N,.E,Q,. 0 GO| 7.0, .2,0 
co tat el aee 10, tis it N 

ee eee eer eee ert 


5.5. CONTINUE STATEMENT 
Function: 


To act as a dummy executable statement. It is used primarily as the terminal 
statement of a DO range. It may be used wherever a dummy executable statement 
is required. 


CONTINUE ; 


Rules: 


(1) The CONTINUE statement does not perform any logical operation. 


(2) The sequence of statement execution is not changed by the CONTINUE state- 
ment. 
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Example: e 


The following sequence changes all negative values in array MAP to their corresponding 
positive values. Because the last logical operation is an arithmetic IF statement, a 
CONTINUE statement is added and made the terminal statement of the DO range. 


po | lo ews 


OLN, Mii 1.0.0 
path 
Pee | 
eee Gece 


Pt de MIA, Pi (K 2,01, 11,0;,)1)0 
2,0) 


M,A,P (}K,) .=, ,—.M,A;P,(.K 


| aol leon, TINUE 


5.6. PROGRAM CONTROL STATEMENTS 


A program control statement either temporarily halts execution of a program (the 
PAUSE statement) or terminates execution of a program (the STOP statement). 


5.6.1. PAUSE Statement 


Function: 





To temporarily halt execution of a program. 
PAUSE 


or 


PAUSE n 


where: nis a string of one to five octal digits (the digits 0 through 7). 





(1) A PAUSE of either form temporarily halts execution of the program. The method 
of resuming execution differs with each processor and is described in the pro- 
gramming manual for that processor. 


(2) The decision of resuming execution is not under program control, but is usually 
made by the operator under instructions from the programmer. 


(3) With the second form of the PAUSE statement, the digit string is displayed or 
accessible by other means. Use of these digits depends upon the particular 
processor used. 


(4) If execution is resumed without changing the state of the program, program 
control is passed to the next executable statement in normal sequence. 
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Examples: 


@ The following sequence checks a list of account numbers in array LIST. If any 
of these numbers is not greater than zero, control is passed to an error routine 
and execution of the program is temporarily halted. 


D,! M,EIN,S,1,ON, Lt ,S,7T/(,2,0,, 3,0,, 4,0) 


pon 
Eeacaad ; 
Fs | |o.0, 110, om 
peeens 


=,1,,,4,0 

DO, 140 L,=,1 3,0, 

ac | lo.8 110, K,=,1 2,0 
LF Lil s,T,( K,, by, M,) |. bE. ,0,) GO, jT,0, ,2,0 


Pol leon ti wwe 
mes e.rilrjovr r,o,u,t,i nile 


5.6.2. STOP Statement 
Function: . 
To terminate execution of the program. 
STOP 


or 


STOP a 


where: n is a string of one to five octal digits (the digits 0 through 7). 





Rules: 


(1) There must be at least one STOP statement in a program to terminate execution 
of the program. 


(2) Action that follows execution of the STOP statement depends upon the particular 
processor used. 


(3) The digits of the second form of STOP statement are not necessarily accessible, 
depending upon the particular processor used. 
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Example: 

The following program has read in a list of telephone numbers into array NUMBER. & 
Each telephone number contains seven decimal digits (no area code). To ensure 

that all values have been read in correctly, each number is checked to see that it 

is greater than zero and not greater than 9999999, If any such error is detected, 

control is passed to an error routine (which may print a message) and execution of 

the program is terminated. The second STOP statement is assumed to be the STOP 

statement encountered during normal execution and is included to show that a program 

may have more than one STOP statement. 


Ree arr es SFOs OOO 
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6. INPUT/OUTPUT AND FORMAT 
STATEMENTS 


6.1. GENERAL 


Input statements fetch data from input and auxiliary storage devices to be used in 


. the program. The input statement is the READ statement. Output statements store 


results obtained in the program on auxiliary storage devices or display the results 
on output devices. The output statement is the WRITE statement. 


The sample program deck of cards of Figure 1-8 illustrates use of the READ and 
WRITE statements. The same program can be used without change for different data 
decks. All that need be changed in a particular application is the deck of data cards. 


FORTRAN also gives the programmer some control over external devices with the 
BACKSPACE and REWIND statements for positioning such devices as magnetic 
tape, disc, and drum units. The ENDFILE statement can be used to demarcate files. 


These statements apply to the transfer of sequential files to and from the processor. 

A file is the entire set of data on the [/O device designated in an I/O statement. A 

file may be subdivided into records. In an 80-column card reader, each record is an 
80-column punched card; in a printer, each line is a record; in a magnetic tape unit, each 
record may have a different size (in characters) up to a maximum which is specified for 
each computer. Tape records are usually separated by a gap which contains no data. 
Each record may be subdivided into fields, the size of which is determined by the 
programmer. 


The term sequential file is used as opposed to a random access or direct access file. 
For example, with standard I/O statements it is not possible to read the fifth record 
of a file directly; it is necessary to indicate that the preceding four records are to be 
passed over. Once the fifth record is read, it may be impossible (as in the case of-a 
card reader) to go back and read the third record; however, this can be done where 
the BACKSPACE or REWIND statement is effective (as on magnetic tape). 


READ and WRITE statements may refer to a FORMAT statement which describes the 
characteristics of the data being transferred. Such statements are called formatted 
statements. FORTRAN alse provides for the transfer of information from one medium 
to another without change or conversion; these I/O statements are called unformatted 
statements. In addition, a READ or WRITE statement usually contains a list that identi- 
fies the items being transferred. 

The remainder of this section describes: 

mw elements of I/O statements 

gs FORMAT statement 

gm formatted READ. and WRITE statements 

a unformatted READ and WRITE statements 


@ auxiliary I/O statements 
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6.2. ELEMENTS OF READ AND WRITE STATEMENTS 


Each READ or WRITE statement may reference a FORMAT statement or specification 
and a logical unit, and may contain an I/O list. 


6.2.1, 


6.2.2. 


Logical Unit Number 


The logical unit number is an unsigned integer that designates the I/O device contain- 
ing the file being referenced. A file may be transferred from one medium to another. [f 
it becomes necessary to access the same file later on, it will have a different logical 
unit number than the one originally used. In previous examples, the integers 1 and 3 
were used to denote a punched card reader and a printer, but these numbers were only 
for use in examples. There is no standard convention for assigning numbers to logical 
units. This information must be obtained from programming manuals for a particular 
processor. 


Input/Output List 


The purpose of an input/output list is to identify transferred items so that they can be 
referenced in the program. A transfer initiated by a READ or WRITE statement is not 
complete unless all items in the input/output have been transferred. It is convenient 
to define an input/output list in terms of a simplé list and a DO-implied list. 


A simple list is a variable, array element, array name, or two simple lists separated 
by a comma. For example, 


V2,ARRAY ,MATRIX(5) 


is a simple list. Previous examples were restricted to simple lists in READ and WRITE 
statements. When an array name appears in a simple list, it refers to all elements of 
that array in the order described in 2.7.4. 


An //O list is a simple list, a simple list enclosed in parentheses, a DO-implied list, 
or two lists separated by a comma. When there is no I/O list in a READ or WRITE 
statement, the I/O list is said to be empty. For example, 


V2,(ARRAY ,MATRIX(S),(NAME)) 


is an I/O list. No I/O list may contain a constant except in a subscript expression or as 
a parameter of a DO-implied list. 


A DO-implied list is a list followed by a comma and a DO-implied specification, all 
enclosed in parentheses. A DO-implied specification has the format: 


i= My ,My,M3 
or 


i= my ,m9 
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pa The italicized parameters and the control variable are the same as those for the DO 
@& statement (see 5.4). For example, the DO-implied list (ARRAY(K), K=3,5) refers to 

: the array elements ARRAY(3), ARRAY(4), and ARRAY(5), in that order. This illustrates 
one of the advantages of a DO-implied list. It enables selected array elements to be 
referenced without the use of a DO statement. The DO-implied list (MATRIX(J), 
ARRAY(J), J=1,3) refers to the array elements MATRIX(1), ARRAY(1), MATRIX(2), 
ARRAY(2), MATRIX(3), and ARRAY(3), in that order. 


Examples of DO-implied lists: 
(1) (ARRAY(J,K), J=1,5,2), K=3,4) 


refers to ARRAY(1,3), ARRAY(3,3), ARRAY(5,3), ARRAY(1,4), ARRAY(3,4), and 
ARRAY(5,4), in that order. 


(2) (ARRAY(J,K), J=1,3) 


refers to ARRAY(1,K), ARRAY(2,K), and ARRAY(3,K), in that order, where K 
was defined previous to execution of the READ or WRITE statement containing 
the DO-implied list. 


(3) ((ARRAY(J,K,M), J=1,J2), K=1,K2), M=1,M2) 


refers to the array elements of ARRAY in their natural order if J2, K2, and M2 
were the declared dimensions. 


(4) An example of an I/O list containing a DO-implied list is: 
A,B, (C,ARRAY(K), K=3,5) 
2 tefers to A,B,C,ARRAY(3),C,ARRAY(4),C, and ARRAY(5), in that order. 


(5) The elements of an array can be referenced in any order. For example, 


((ARRAY(,K), K=1,L), J=1,M) 
interchanges the order of subscripts of ARRAY. 


(6) An example of a DO-implied list within a DO-implied list is: 
(ACJ), [=1,10,2), B(J,3), J=1,K) 
The order of reference is clearer if shown as follows: 
DO J=1,K 
[ DO [=1,10,2 


ALJ) 
B(J,3) 


Thus, the order of reference is: 


A(1,1),A(3,1),A(S5, 1), A(7,1),A(9,1), BQ1,3), 
A(1,2),A(3,2),A(5,2),A(7,2),A(9,2), B(2,3), 


A(l .K), (3, K), ACS, K),AC7 ,K),A(9,K), B(K,3) 
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6.3. FORMAT STATEMENT 


Function: , S$ 





To provide conversion and/or editing information between the internal representation 
and the external character strings in conjunction with a formatted READ and/or WRITE 
statement. 


Format: 


(440124929 + + + tt 2p_-1'n 99) 


where: each q is a series of one or more slashes (/) to act as a record demarcator and 
may be omitted. 
each t is a field descriptor or a group of field descriptors, 
each z is a field separator. . 
The part (q4¢ 24 i Be q) is called the format specification. 





(1) A FORMAT statement must have a label. It can be referenced by one or more 
formatted READ and/or WRITE statements. 


(2) If a formatted READ or WRITE statement references a format specification by 
array name, only the format specification (including the enclosing parentheses) 
must be the first item in the array. Any information in the array after the rightmost 
parenthesis is ignored. For example, . 





indicates that K2, K3, and K4 are to be interpreted as integers occupying five 
character positions, and that K5, K6, and K7 are to be interpreted as integers 
occupying eight character positions. The READ statement could also have been 
written as: 


DIM ,EIN,S,! ON ,1,N 7T.()2,) 


| | toate 1N.7.(,1,),/,5,4](,3,1,5,,,/,,,1,N,T/(, 2,),7,4,H,3,1,8,) |/ 
Raden | | 


|. | [reap (,1,,,1,.N,T,) K, 2), ,&,3,,,K,4,,,K,5,,)K,6,,,«,7 


where the DATA statement initializes the contents of array INT to the required 
format specification, using Hollerith constants. 
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The only restriction on the use of a format specification in an array is that the 
format specification cannot contain a Hollerith field descriptor of the form nH. 


Biank characters may be used freely in the FORMAT statement or a format 
specification, except with a Hollerith field descriptor. 


If there is an I/O list in the formatted READ or WRITE statement, the format 
specification must contain at least one field descriptor other than nH or nX. 


Format control (control by a FORMAT statement or specification) is initiated when 
execution of a formatted READ or WRITE is started. (Further details on format 
control are furnished in 6.6.) - 


The first character of a record to be printed is not printed; it is used as a form 
control character as shown in Table 6-1. 


CHARACTER VERTICAL SPACING BEFORE PRINTING 7 


ONE LINE 


TWO LINES 


TO FIRST LINE OF NEXT FORM 


NO ADVANCE 





Table 6—1. Form Control Characters 


There are many ways of assuring that the required form control character is the 
first character. The safest way is to use 1HB or 1X for the blank character and 
1HO, 1H1, 1H+, respectively, for the form control characters in Table 6—1. However 
the sequence : 


= 
Az 
" 


Roos tPF war 


FLOR MjJA,T (I 2,,,-,-)- 1) 


will have the same effect as if the FORMAT statement had been written 





In both cases, the printer will advance to the start of the next line and print the 
digit 1 in the first print position of the line. 
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However, if the FORMAT statement had been written : 





the printer would have advanced to the first line of the next form before printing the 
remaining items of the I/O list. 


(7) A format specification without any field descriptors is valid. For example, the 
format specification (/////) is valid. It causes the printer to advance five lines. 


Record Demarcator 


The term record demarcator refers to the one or more consecutive slashes (/) that 
appear anywhere in a format specification. If there are n slashes at the beginning 
or at the end of the specification, a records will be skipped; if there are n slashes 
anywhere else, processing of the current record is terminated, and n—1 records are 
skipped. 


For example, in reading punched cards, the sequence 





will obtain the value of K2 from the current card and the value of K3 from the next 
card. The sequence 





causes skipping of two punched cards before obtaining the value of K2 from the third 
card; skipping the fourth card; obtaining the value of K3 from the fifth card; then skipping 
two cards for processing of the next READ statement (if any) for the same file. 


Field Separators 


A field separator is either a comma or a series of one or more consecutive slashes. 
It is used to separate field descriptors in the list of format specifications. The 
slash(es) also acts as a record demarcator and ends processing of the current record. 


For example, in printing information, the sequence 





will cause printing of the values K2 and K3 on one line, and the value of K4 on the 
next line. 
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6.3.3. Field Descriptors 


Function: 


To indicate how items from/to an input/output device shall be represented internally/ 
externally and provide editing information in the form of spacing, form control (printer 
only), and Hollerith data. 






xPrFw.d 
xPrEw.d 
xPrGw.d 
xPrDw.d 
tlw 
rLw 
rAw 
aHhyhy ...A 
aX 









n 









where: the letters F, E, G, D, I, L, A, and H indicate the conversion and editing 
and are called conversion codes. 


w and n are unsigned integers greater than zero indicating the number of 
character positions in a field of the external medium. 


d is an unsigned integer constant indicating the number of.digits in the 
fractional part of the external character string (except for the G conversion 
code). 


r is an optional nonzero unsigned integer, the repeat count, indicating how 
many times to repeat a basic field descriptor that follows it. 





xP (optional) indicates scaling, x being an unsigned integer representing 
the scale factor. 






A represents a character from the processor character set. 


Rules: 


(1) The field width must be specified for all descriptors. 


(2) For descriptors of the form w.d, d must be specified even if it is 0, and w must 
be greater than or equal tod. 


(3) All output is right-justified in the output field width specified by w, preceded by 
leading blank characters (w permitting). 


(4) The number of characters produced on output cannot exceed the field width, w. 
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6.3.3.1. Blank Field Descriptor 





The field descriptor for blank characters is nX. On input, n characters of the input © 
record are skipped; on output n blank characters are inserted in the output record, ~ 
except for the carriage control character. 


6.3.3.2. Numeric Data 
There are five conversion codes for handling numeric data: 
m the rlw code for handling integer type data 
@ the xPrFw.d code | 
mw the xPrEw.d code for real type and complex type data 
@ the xPrGw.d code 
gw the xPrDw.d code for handling double precision type data 
-On all numeric conversions, a blank character in the specified field is treated as 
a 0. A blank field is treated as the integer zero. Depending upon the processor and 
the field width specified, positive values may be preceded by a + or no character 


position; a negative value requires a character position in the field width for the 
minus sign. Leading blank or zero characters are not significant. 


6.3.3.2.1. Integer Type Conversion 


The basic field descriptor Iw indicates that the external field occupies w positions 
as an optionally signed integer and is represented internally as an integer type 


item. @ 


On input, the extemal field may contain only a sign followed by digits, with 
blank characters anywhere in the string. No other characters are permissible. If 
the sign is plus, it may be omitted; if the sign is minus, it is required. On 
output, the external field consists of blank characters (if necessary) and a minus 
sign or an cptional plus sign followed by the magnitude of the intemal integer 
type value. ° 





Examples: 
(1) An input card contains 


Position 1 8 


BbodH1DS2 
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and the following statements are used: 





The first five columns will be read to obtain the value for NAME2, which is 
1. Then the next three columns will be read to obtain the value for NAME4, 
which is 2. 


If the following statements are used instead: 


P| PRE ADE G110.) NALMEL2L UNAM, Ed, NAIM EL6 


FO .RMjA,T 14. t2,, 11, 3 


the first four columns will define the value of NAME2 as 0; the next two 
columns will define the value of NAME4, as 10; the next three columns will 
define the value of NAME6 as 20. Note that trailing blank characters are 
significant and treated as zeros. 


If the following statements are used instead: 


Pe | fee A.C 1.1.0.) NL AMLEL? 
a 


| |e ,0,RMIJA,T 194) 


the first nine columns will define the value of NAME7 as 10020. 
(2) A punched card contains 
1 2 3 4 


8 
Position 1 0 0 0 0 _0 
( 6bbb 1234665 -12346-5b 123456 -12365+-B515253465b ... B 





and the following statements are used: 


| | [ot MEIN, S,1,0N LAR ALY,(.2,), 
veve LEVOTESETOETE 


poor 1 10 K, 1 LA RiAi Yi. Ki 2h, (Ki 3 Ki 4 


FO .R.MIA.T 3. X U5 3,X% 6 ma:) 115 11,0 { 
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The first three columns are skipped; 

the next five columns read in the value for K1, which is 1234; 

the next three columns are skipped; 

the next six columns define the value of IARAY(1), which is -12340; 

the next eight columns define the value of IARAY(2), which is -12340, 
the next five columns define the value of K2, which is —123; 

the next 10 columns define the value of K3, which is 102034; 

the next five columns define the value of K4, which is.0. 

The list is now satisfied and any remaining items on the card are ignored. 





(3) After the following sequence is executed, 





the printed line is: 


106b5- 10 starting in print position 1. 


If the output list items in the WRITE statement had been 


WR TE 3 1,0, ) MAMIE, 4... NAME 2 7 


an error condition would result. The first print character is a minus character, 
which is not a legitimate form control character. If this minus sign is interpreted 
as a form control character (in some processor implementations), there is not 
enough field width provided for printing —10. In this case the FORMAT state- 
ment requires a change such as: 


ee tol eommat Me BR ee a Ps 


to ensure that NAME4 will be printed correctly. 
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6.3.3.2.2. Input of Real Type Data 


For input numerical data to be represented internally as real type data, three 
conversion codes are available: the F, E, and G codes. Operation of these three 
basic field descriptors is identical for input data (their operation differs on output 
data). 


The basic form of the extemal input field consists of an optional sign (optionally 
preceded by blank characters) followed by a string of digits and blank characters 
with one optional decimal point anywhere in the string. This basic form may be 
followed by an exponent, having any of the following forms: 


@ a signed integer constant 
w the letter E followed by an optionally signed integer constant 
gm the letter D followed by an optionally signed integer constant 


A decimal point in the basic form overrides the decimal point specified by the 
d designator. 


Example: 


A punched card contains 


1 2 3 4 5 6 
Position 1 0 0 0 0 0 
3456789 345.6789 3456789+3 -67.89E-3 -6789D10 -67.89D-2 


and the following statements are used: 


Al is defined as 3456789.0; 

A2 as 345.6789; 

A3 as 3456789.0; 

A4 as - 0.06789; 

A5 as - 678900000000.0; 

A6 as 0.6789 (note that the minus sign is lost). 


11 
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6.3.3.2.3. Output of Real Type Data 


On output of numerical data that is internally represented as real type data, three 
conversion codes are available: the F, E, and G codes. The output form is different 
for the F and E codes; the G code results in form similar to either the F or E code. 


On output, the Fw.d results in an output field consisting of blanks (if necessary), 
a minus sign or optional plus sign (depending upon the processor), followed by a 
string of digits containing a decimal point and a fractional part rounded to d 
fractional digits. 





On output, the Ew.d results in an output field of the form: 


0.x; ae ae xgE+tyyy9 
or 
(0.x, ecw XqgtV¥ 17273 
where: the choice of form depends upon the processor. 


X, -. + Xq are the d most significant rounded digits of the value to be 
transmitted, 


y represents a digit of the decimal exponent, 


the plus sign following the letter E may be represented by a blank 
character (depending upon the processor); 





t is either no character position or a minus sign; the 0 may be replaced 
by no character position (depending upon the processor). 


On output with the Gw.d basic field descriptor, the form of the external field 
depends upon the absolute value of the internal real type value. If N is the 
absolute magnitude, the effect conversion is as follows: 


MAGNITUDE OF N EFFECTIVE CONVERSION 
O1<N<l F(w- 4) .d,4X 
1<N< 10 F(w-4) .(d-1) , 4X 
10d-2<n<10¢-1 F(w-4) .1, 4X 
10%1<Nn <10¢ F(w-4) .0, 4X 

otherwise Ew.d 


For effective F conversion, the absolute value of the real type item in storage 
must be equal to or greater than 0.1 and equal to or less than 104. 
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g four positions for the decimal exponent; 
® one position for the decimal point; 
B@ one position for the 0; 


m@ one position for the sign. 


Therefore, this general rule pertains to these basic field descriptors: w must be at 
least 7 greater than d on output. 


Examples: 
(1) The successive elements of array table are 


.001234 

.001234 
—12.345 
~—12.345 


and the following statements are used for printing: 


@ Po | ye T203,..1 0) TABLE | 
| tol {ror aiact TX Fi 15) 17, E15) 7, Fi U5) 21 E51 7 


The printed line will be (subject to processor options): 


1 3 4 6 
Position 1. 5 0 5 ot 
666 bbH56.00 12340b6b0.1234000E-02bbbbbbb5b-1 2.355-0.1234500E +02 


(2) Examples for the basic field descriptor Gw.d on output for different values 
of the transmitted item are: 


ITEM CODE EFFECTIVE PRINTED 
VALUE DESCRIPTOR DESCRIPTOR ITEM 
—.0123456 G10.3 E10.3 -0.123E- 01 
- 23.456789 G11.2 F7.0,4X Bbb- 23 BESS 
- 23.456789 G11.3 F7.1,4X Bb- 23.55555 
- 23456789 G11.5 F7.3,4X - 23.4575bbB 
—123.456730 G11.2 E11.2 BB- 0.12E+03 
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6.3.3.2.4. Double Precision Type Conversion @ 


On input, the basic field descriptor Dw.d is used for items to be represented 
internally as double precision type data; on output, it is used for items that are 
represented internally as double precision type data. 


The form of the external input field is the same as that for real type conversion 
(see 6.3.3.2.2). The form of the external output field is the same as that for the 
Ew.d field descriptor (see 6.3.3.2.3) except that the letter D may replace the 
letter E. As with real type numbers on input, a decimal point in the input value 
overrides the d specification. 


The advantage of the D descriptor is that it can store and output more significant 
digits than the E descriptor can. 


Rule: 
The internal field must be explicitly declared as double precision type. 


Example: 


The following short program illustrates use of the D conversion code: 


DOUBLE PREC IS!IO.N, X,,,X%,2 


Pere Eee 1 1,0 xX 
igs 


F,O.R MjA,T Di1,5,-45 

X,2, ,=( ,2,.,0 * x 

WR! TIE, (,3,,,2,0,) X,2 
| 2,0! |F,0,R,MIA,T,(.D,1,7,-,5,) 
estes 


$s ,7,0,P \ 


E.N,D 
| 
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& 6.3.3.2.5. Complex Type Conversion 


A complex data type item is represented internally as two consecutive real type 
items. Therefore, for each of the pair (the first is the real part and the second is 
the imaginary part), a real type conversion code is required.. 


Requirements for the external input field are the same as for real type conversion. 
Each successive pair of items is interpreted as a complex type item. 


The external output field requires a separate real type conversion code for each 
part of the complex value. 


Rule: 


The internal field requires explicit type declaration. 


Examples: 
(1) Input 


A punched card contains 


P Position 1 0 0 0 0 0 
-123.566 -123-3 45.6-02 








and the following statements are used: 


é C,OM,PiL,E,X, .C,P,X,1,,,CiP,X,2 


(iets ee 
a 
m 
> 
Oo 


1,,,1,0,) CP} X,1,,,A,,,C,P, X,2 
1.0] [F,O,RMIA.T, (,F,1,0.., 0,1, F1,0,-,1,,,F,1,0,-, 14, £ 1,0; ,0, 5,6 1,0) .,0,) 





~123.56- 0.0123; is read into CPX1, 
0.456 into A; 
1.23+0.01237 into CPX2. 


| Since the input list is satisfied, the last G10.0 of the FORMAT statement is 
disregarded, 
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(2) Output 





After execution of the following program 


COMPILE xX, CP.X,2,, CiPX.4 


ae CPX 2 (,-,1,2,3,., 5/6 -,3,4,5,£,-,3,) 
Tt leeax4 7,8,9,. 5,6), .-:8,9,0, .,1 
eet 


See 


uo 


P| ls 7 bs 
st [FORMAT E43. ,5,,:21X% Fi t2.- 517, 5X, F511 3 / Fy 8-24 2X Fi hy I 
cent 


$,T,0.P 
END 


The first print line will be: 
—0.12356E+03bb—0.34500E +03bbbtb—. 456 


The second print line will be: 
6789.5655—890.1 


Note that the first character of each print line (a blank character) is interpreted 
as the form control character and is not printed. 


6.3.3.3. Logical Type Conversion 


On input, the basic field descriptor Lw is used for items to be represented internally 
as logical type data; on output, for items that are represented internally as logical 
type data. 





Rules: 
(1) The internal field must be explicitly declared as logical type. 
(2) On input, the external field consists of optional blank characters, followed by 


either a T for true, or F for false, followed by optional characters which are 
ignored. : 


(3) On output, the external field consists of w- 1 blank characters followed by 
a ToranF. 
Example: 


An input card contains 


1 1 
Position 1 0 9 
TON FOFF 


and the program contains the sequence: 
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| | feye.6.tyepat S.W,T, CH, 2).,5,¥,7, CGH, 4 
ee ee 1.45 S,W,T.C/H, 2, S$, T.C,H, 4 


erie L,1,0,,,b,9 


| wiki tie 3 S$.W,T!C.H SW 


The printed line will be: 
BHDSDHDHTHSDHDODOF 


6.3.3.4. Hollerith Field Descriptors 


Hollerith information may be transmitted by the field descriptors nH and Aw, 
as follows: 


(1) On input, the nH reads Hollerith data into the n characters following the nH 
field descriptor in the format specification. On output, it writes the last 
previously defined n characters of the format specification following the nH 
field deseriptor. 


(2) The Aw field descriptor reads or writes w Hollerith characters into or from 
an element of the I/O list. 


Rules: 


(1) The nH field descriptor must not be used in a format specification if a 
READ or WRITE statement references the format specification by array name. 


(2) For the Aw field descriptor, let g be the maximum number of characters that 
can be represented in a single storage unit (see 2.5.1). If wis greater than 
or equal to g on input, the rightmost g characters will be accepted and any 
remaining leftmost characters will be lost. If w is less than g on input, the 
w characters will be accepted left-justified internally with remaining storage 
positions filled with blank characters. If w is greater than g on output, the 
external field will contain the g characters right-justified with the remainder 
of the field filled with blank characters..If w is less than or equal-to g, the 
external field will contain only the leftmost w characters from the internal 

a representation. 


[-xamples: 


For explanatory purposes, it has been assumed that no more than five characters 

can be stored internally in an integer, real, or logical data type item. The number of 
characters that can be stored in a single storage unit varies with processor implementa- 
tion; there is not standard capacity. 
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(1) In this program, the aH field descriptor causes printing of specified characters. 





K,=,1 
are Se 
| fuse 2,5 ,+.K 
UL iweretre 3..,1,9 SiU.M 
Pavel’ mo ealiay : O.TA =iul ) 
Peewee ETT 
| den, 





The printed line will be: ee 
TOTAL = 3 


(2) Consider the following short program which reads in one data card containing: 


1 2 3 4 5 
Position 1 oe 0 0 0 0 0 
CURRENT 1.2 56.78 2.3 34.56 : 


i 


|. | lot Me1N, $10.6 RS.T.RIC1,0,),,,¥,0,b,7,${(,1,0)) », ApM.P.S,({1,0)) 
| leew. L110 RiSITRGK 1 VO TIS GK Keeyhi2 


E.O.R MjAT 610 HBV O LIT AG ESB FT Cl. 0, Fit OO FINO. Oe Fi. 0 





| |p.o 210, ,K; =)3 
| aM.P.SiCK)=.VOL.T.S$ (1K) /.RS.7.RCK) 
Pa | leet 31..330; AIM.P =: 
RE S$, 7,0,P 
END 


iy 


ua 


_ The printed output is: 
CURREN 1555500047 .325655615.03 


The READ and WRITE statements are processed as: 


(a) After the card was read in, the FORMAT statement became, in effect, 
10 FORMA T(LOHBCURREN Tbb,F10.0,F10.0,F10.0,F10.0). 
On printing, the first B was interpreted as a carriage control character 
and was not printed. 


(b) It was not necessary to define all elements of the arrays since these 
were not used. 


(c) Although the WRITE statement required only two I/O list items, its 
FORMAT statement contained four field descriptors. As soon as an [I/O 
list is satisfied (all items are accounted for), format control is terminated. 
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(3) 


Peeea 
Pest 
ceewed 
ieee 
19 | 
REeee i 
ieewe | 
Peer 
Peas i 
Reeees 
ets 
Deane 
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The next example shows how the coding of repeated field descriptors can be 
simplified. 


The program of the previous example is simplified and made more I/O oriented by: 

(a) reading in the format specification; 

(b) reading in the number of array elements to be processed in the DO-implied. 
lists; 


(c) simplifying coding of repeated field descriptors; 


(d) reading in the numbers for input/output devices. 


The first input record from the card reader is organized as follows: 


(a) In columns 1 through 8, (8F10.0) is used as the format specification for 
input data for arrays RSTR and VOLTS. This specification is named 
INFORM. 


(b) In columns 11 through 31 (Ai,15X,A5,A2/(1X,F10.0)) is used as the 
output format specification, OUTFRM. 


(c) In columns 39,40, a two-digit unsigned integer denotes NIN, the input 
device for data for arrays RSTR and VOLTS. 


(d) In columns 44,45, a two-digit unsigned integer denotes NOUT, the output 
device for the computed data in array AMPS. 


(e) In columns 49,50 a two-digit unsigned integer (not greater than 10) indicates 
NMBR, the number of array elements required in each of the arrays. 


(f) In columns 51 through 57, CURRENT denotes output heading. 


(Statements are continued for illustrative purposes only.) 


Di M,EIN,S, 1, OLN, RLS. TRIO 1, 90 V,0,L,T,S$1/(,1,0 ALM.P.S,(/ 1,0 
{NL FjO RM( 2 OU TIF,.RM,( (5 TI TILE (2 

D A,T A; CN TRL 7/14 0)/ 

R,E,A,D 1,.11,9 UNL FIO,RM,, 0,0 TF, RM TiN, iN, O,U,T 
NM BIR Tel Tk 

FOR MIA,T.(.2,A,5,, ,5,A:51, 13, (3%, 11 2 Ai5,,1A,2 

RE,ADICIN IN INGE .O, RIM R,S,T, RC iK V,O,L,7T:5,(K 
K,= li, NUM, BR 

1 N.M.BLR 

VOLTS, (jM,),./,R,8,T,R(.M,) 


D,O, ,2;0, .M 
AMPS} (.M,) 


a“ 


WLR TIE, (.N.O.U.T,, 0.0, TIF.RM CN TRL TITLE 
A MIP.S,( J J,=,1,,|NM BR 
$,T,0, P ! 


| len. jes i 
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Up to three of the remaining input records can contain, as before, the input i 
values for voltages and resistors (RSTR and VOLTS) in fields of 10 positions @ 
per value. The input device for these records is designated by the programmer. 

The output contains a header and thé calculated output values (AMPS), one 

value to a record. The output device for these records is also designaied by 

the programmer. 


The repeat specifications and the use of seemingly redundant parentheses within 
a format specification are described in 6.3.3.5 and 6.6, but their particular applications 
in this program are described in the following paragraphs. 


In lines 4 and 5, the I/O list requires12 items: 2 for INFORM, S for OUTFRM, 1 
each for NIN, NOUT, and NMBR, and 2 for TITLE. 


In line 6, the two items required for INFORM are covered by 2A5, which could 
also be written as A5,A5; and the five items required for OUTFRM are covered 
by SAS. However, a group repeat designator is. used in this line. This indicates 
successive repetition of the group within the parentheses. The group repeat 
designator 3(3X,12) could also be written as 3X,12,3X,12,3X,I2. 


Line 7 refers to the format specification in INFORM, which is (8F10.0). The 8 means 
that there can be 8 fields to a record with the basic descriptor F10.0. However, as 
many as 20 values may be required: 10 for RSTR and 10 for VOLTS. If the DO-list 
indicates that more than 8 fields are required, the first 8 fields will be read in, a 
new record automatically started, and the same specification (8F10.0) used, and so 
on, until all values required for the I/O list are read in. 





The WRITE statement refers to the format specification in OUTFRM, which is (Al, @ 
15X,45,A2/(1X,F10.0)). The Al refers to the carriage control character called 
CNTRL, that was initialized to 1HO. It is not permitted to use an nH field descriptor 
directly in a format specification that is stored in an array; therefore, this indirect 
method must be used. (An alternative method would be to read it in.) The slash in 
this format specification indicates the end of a record. The field descriptors that 
follow apply to the next record which permits only one F10.0 field descriptor to a 
record, applying to the output results in AMPS. If AMPS contains more than one 
value (indicated by NMBR), it will automatically start a new record and use the 
1X,F10.0 for each new value and record. If the 1X,F10.0 was not‘enclosed within 

its own parentheses, format control wceuld have been started at Al, instead of 

after the slash each time a new record was required. 
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Repeat Specifications 


Repetition of a field descriptor (except nH and nX) in a format specification is 
accomplished by using the repeat count r (see 6.3.3) immediately before the basic 
field descriptor. For example: 


5L3 is the same as L3,L3,L3,L3,L3. 


Repetition of a group of one or more field descriptors and/or field separators is 
indicated by enclosing the group in parentheses and optionally preceding the 
group with an integer constant. If no group repeat count is indicated, its value 
is implicitly 1. This form, with or without a group repeat count, is called a 
basic group. For example: 


2(2X,215,F 10.0) is the same as 2X,15,15,F10.0,2X,1S,15,F 10.0. 
2(//) is the same as ////. , 
2(4HHALT) is the same as 4HHALT,4HHALT. 


A further grouping may be formed by enclosing field descriptors, field separators, 
and/or basic groups within parentheses. Again, a group repeat count may be 
either specified or implicitly 1. The parentheses around this further grouping are 
called second level parentheses. Parentheses are permitted in a FORMAT state- 
ment only to a second level; therefore, this further grouping may not be contained 
in another grouping by parentheses, 


The first left parenthesis and the last right parenthesis required of a format 
specification are not considered group delimiters. 


Example: 





The elements of the I/O list are read in as: 


M(1) with specification [2 
M(2) with specification [2 
N2 with specification 12 
skip three character positions 

A(1) with specification F6.3 
B(1) with specification F6.3 
skip two character positions 

c(i) with specification F6.4 
A(2) with specification F6.3 
B(2) with specification F6.3 
skip two character positions 

C(2) with specification F6.4 


skip four character positions 
N with specification [2 
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6.3.3.6. Scale Factor 


Input and output using the F, E, G, and D conversion codes can be scaled up or @ 
down by a power of 10 if the conversion code has the form xP immediately pre- 

ceding a repeat specification (if any), where x is an integer constant that may be 

preceded by a minus sign. The effect of x, the scale factor, is to multiply the 

corresponding I/O list item by a power of 10. 


Rules: 


(1) When format control is initiated (see 6.6), a scale factor of zero is established. 
Once a scale factor has been established, it applies to all subsequently inter- 
preted F, E, G, and D field descriptors until another scale factor is encountered 
and then established. 





(2) The effect of the scale factor is temporarily suspended for input with F, E, G, 
and D conversions that contain an exponent in the external field. 


(3) For input with F, E, G, and D conversion and no exponent in the external field, 
the internally represented value is the external value divided by 10*, 


(4) For output with F conversion, the external value is the internal value multiplied 
by 10%, 


(5) For output with E or D conversion, the basic constant part of the output is 
multiplied by 10* and the exponent is decreased by x. 


(6) If the effective use of E conversion is required for output with G conversion 


(see 6.3.3.2.3), the scale factor has the same effect as for E conversion. If & 
effective use of the F conversion code is required, the scale factor has no 
effect. 

Example: 


Two identical punched cards have all values right-justified in their fields of 
10 positions each, as follows: 


columns 1 through 10: 100.21 
-- columns 11 through 20: 100.21 
columns 21 through 30: 100. 21D-3 
columns 31 through 40: 100.21 
columns 41 through 50: 567890 
columns 51 through 60: 1234SE-2 
columns 61 through 70: 1234.5 
- columns 71 through 80: 654321.321 


with the following statements: 


Y | DI MENS LON AC7),, B7) 

| | [DOUBLE PRECISION AS, B8 

eee: 

tN ate ge eel 
| | READ 1, 10) A, AS, 8, B88 


| tol [FORMAT (F 10.0, 2)P2£10.0,2610.09,0P2F10.5,2PD10.¢4) 
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A(1) is read in as 100.21 because the scale factor is established as 0 (unless 
stated otherwise) at the start. 

A(2) is read in as 1.0021 because the scale factor of 2 is effective. 

A(3) is read in as 0.10021 because the scale factor of 2 is temporarily suspended 
by the exponent in the extemal field. 

A(4) is read in as 1.0021 because the scale factor of 2 is still established and 
effective. 

A(5) is read in as 5678.9 because the scale factor of 2 is still established and 
effective. 

A(6) is read in as 0.0012345 with the scale factor now established as 0. 

A(7) is read in as 1234.5 because the scale factor is still established as 0. 

A8 is read in as 6543.21321 because a scale factor of 2 is now established and 
effective. 

Because a new record start is required by the list, the next card is read, but there 
has been no termination of format control. 

B(1) is read in as 1.0021 because the scale factor of 2 is still established and 
effective. 

B(2), . . . ,B(7), and B8 are read in with the same values as the corresponding A’s. 


FORMATTED READ STATEMENT 
Function: 


To initiate input of data from a specified input device and to scan and interpret this 
data in accordance with a format specification. 


READ (u,f) I/O list 


or. 


READ (u, f) 


where: u is either an integer constant or an integer variable that identifies an input 
unit. 
f is either the statement label of a FORMAT statement or the name of an array 
that contains the format specification. : 
I/O list is described in 6.2.2. 





Rules: 


(1) The number of records to be read depends upon the I/O list and the format speci- 
fication. 


(2) Format control is started and terminated in accordance with the rules of 6.6. 


(3) There are no standard conventions for assignment of integers to input devices; the 
integer for a particular input device depends upon the computer being used. 
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Examples: 


(1) 








Regardless of the input device represented by INPUT, one value is read in for A, 
then seven values for ARRAY. A new record is automatically started, one value 
is read in for B, followed by seven values for BRRAY. 


(3) 








The first 10 characters of the record from INPUT replace HEADINGDDBB in the 
format specification; the same FORMAT statement can then possibly be used in 
a WRITE statement to supply the header. @ 





(4) 


DI MENS 1 ON ARRAY (1) 


| AT AL ARRAY CT) Ji HL) / 
eenee se 


pris | | 
| | REA DC UN, ARRAY) 


3 The format specification for the READ statement is obtained from the array named 
ARRAY. The effect of the READ statement is to advance the input device by two 
records. 











Pees | INT 
| | [wri tec Nout, 10) A,B,C | 


|, yo] TF,0,R,MIA,T CT HAT, 15,X1,,1,5,H NAME, OF, PROB LEM, /j1,6,%,,,3,F,1,0,. 41) 
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6.5. FORMATTED WRITE STATEMENT 


Function: 


To indicate the output of data to a specified output device in accordance with a format 
specification. 


WRITE (u,f) [/O list 


or 


WRITE (u, f) 


where: u is either an integer constant or an integer variable that identifies an input 
unit. ; 
fis either the statement label of a FORMAT statement or the name of an 
array that contains the format specification. 
[/O list is described in 6.2.2. 





Rules: 


(1) The number of records to be transmitted depends upon the I/O list and the format 
specification. 


(2) Format control is started and terminated in accordance with the rules of 6.6. 


(3) The first character of a print record is not printed; it is interpreted as a form 
control character (see Table 6-1). In standard FORTRAN the printed line may 
start in position 1; however, many printers indicate the presence of the control 
character by printing a blank character, so that printing may start only with posi- 
tion 2 of a print line. This characteristic should be checked for a particular 
processor. 


(4) There are no standard conventions for assignment of integers to particular output 
devices, this depends upon the processor being used. (In examples, 3 is used to 


indicate a printer.) 


Examples: 


(1) 


PROBLEM) 





The printer starts at the first line of the next form, spaces 15 positions, and prints 
NAME OF PROBLEM. 


(2) 
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The output device, if magnetic tape (as an example), will contain the character 

1, followed by 15 blank characters, followed by the characters NAME OF PROBLEM, 
all in the first record. The second record will contain 16 blank characters followed 
by the values of A, B, and C, each in accordance with the F10.4 conversion code. 
These records may be sent to a printer later on. 





6.6. FORMAT CONTROL 


The following rules describe the relation between format control and the I/O list 
(if any) of the formatted READ and WRITE statements. i 


(1) Format control is initiated with the start.of each execution of a formatted READ 
or WRITE statement. 


(2) When the format control encounters an I, F, E, G, D, A, or L basic descriptor in 
a format specification, it determines if there is a corresponding element in the 
1/O list. If there is no such element, format control is terminated; if there is a 
corresponding element, it is converted and transmitted, and format control proceeds. 


(3) If format control proceeds to the outermost right parenthesis of the format speci- 
fication, it determines whether another I/O list item is to be transmitted. If not, 
format control is terminated. If another list item is to be transmitted, format control 
starts a new record and control is transferred to that group repeat specification 
(which may be an implicit 1) terminated by the last preceding right parenthesis 
or, if none exists, to the first left parenthesis of the format specification. This 
action does not affect the scale factor that has been established. 





For example, if there are list elements to be transmitted, a new record is started & 
and format control continues at 


Ice atye se? ae ch to gee Pea esse 


where: xis a group repeat specification (which may be an implicit 1). 
If the format specification does not contain any inner group parentheses, 
format control continues at . 


id 


Co ned 


> (4) There is no corresponding I/O list element for an X or H basic descriptor. An 
I/O list requires at least one field descriptor other than nX or aH. 


Examples: 


(1) 





The eight values of ARAY2 are obtained from one record. 
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(2) 


| td eeaD (.N,P,U,T,,,1,0 GAR AY,2 K, *K, =, 1,,,8,) 
| vol le.oR MAT. (.F1,0,..5,) 


The eight values of ARAY2 are obtained from eight successive records. 


(3) 


WAC UCTARCU UR MERRITT Cr a nin on TT 
eer 


FLOR MIAT, (.2,F,.1,0,.,5,/1/3:F,1,0,.,5 


The first two values for ARAY2 are obtained from the first record; the next three 
values from the second record; the next two values from the third record; and the 
last value from the fourth record. 


(4) 


| ReaD P.U,T,,,1,0,) (,ARAY,2,(,K,)],, ,K,=,1,,,8,) 
| tol [FORMIALT 2,F,1,0,.,5,/1(,3,F,1,0,..5 


The first two values for ARAY2 are obtained from the first record; the next three 
values from the second record; and the last three values from the third record. 
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6.7. UNFORMATTED WRITE AND READ STATEMENTS 
Function: 


To transmit the exact binary configuration of I/O list elements to and from an external 
device. 


WRITE (u) [/O list 


and 


READ (u) [/O list 


or 


READ (u) 


where: wis an unsigned integer constant or an integer variable designating an input 
or output device. 


I/O list is described in 6.2.2. 





(1) An unformatted WRITE statement transmits the exact binary internal representation 
to an external device. An I/O list is required with an unformatted WRITE. 


(2) An unformatted READ (of either form) can only be performed for records created by 
an unformatted WRITE. 





(3) The unformatted READ with I/O list transmits items until the list is satisfied, 
provided that the record contains at least as many items as required by the I/O 
list. 


(4) An unformatted READ without I/O list can be used for positioning the file to the 
next logical record. ; . 


(5) A formatted READ should not be used for a record created by an unformatted WRITE. 


Example: 


Execution time of unformatted READ or WRITE statements is much faster than the 
corresponding formatted statements. They are generally used to write out large lists 
on temporary files and read them back when required by the program so that main 
memory is free for other processing. For example, 
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DJ! M.EIN,S,1,0,N, JA. RR, ALY, (1 0,0,0 


nama LCeace TO Pee? A[R,R,A,Y 


6.8. 


6.8.1. 





R,E,A,D)(,M,T,A,P_E,), ARIR, AY 


where: MTAPE represents a temporary file on magnetic tape, drum, or disc. The binary 
contents (1’s and 0’s) are written and then read back just as they appeared in the 
1000 storage units (see Table 2—4) reserved for ARRAY. 


AUXILIARY INPUT/OUTPUT STATEMENTS 
The auxiliary input/output statements are: 

ws REWIND statement | 

a BACKSPACE statement 

g ENDFILE statement 


REWIND Statement 
Function: 


To cause an input/output device to be positioned at its initial point. 


REWIND u 


where: u is either an integer constant or an integer variable identifying an input/ 
output device. 





Rule: 


The REWIND statement is applicable to such I/O devices as magnetic tape, disc, 
and drum units. Its implementation on devices such as card readers and printers 
depends upon the computer being used. 
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Example: 


Poe 10] feo M ATC 1) 
eee s : 


een eee: 
| | dee Abi M.T.ALPLE 1,0,) 





ft dk E,A,DI( .M,T,A,P,E,, ,1t,0,) 


anes eee" 7 : 


| tee tiW.D, MAPLE 


> 


See seaare | 
Td le E,A.DI(.M.T,A,P,E,,,1,0,)}. 


The first READ references the first record on MTAPE, the second READ references & 
the second record, and the third READ references the first record. 


6.8.2, BACKSPACE Statement 
Function: 


To backspace one record on a specified input/output device. 


BACKSPACE u 


where: u is either an integer constant or an integer variable identifying an input/ 
output device. 
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Rules: 


& (1) The BACKSPACE statement is applicable to such I/O devices as magnetic tape, 
disc, and drum units. Its implementation on devices such as card readers and 
printers depends upon the computer being used. 


(2) If the unit identified by u is already at its initial point, the BACKSPACE state- 
ment has no effect. 


Example: 


Lo iol [Rog wat 


Femees 
Fotis ale 
eau lee M,T,A.P.E,,,1,0 
Beers f . 

ene 


Reeee it 3 
poe MEAIB 1,0,) 
rode 


fi | [BAC KIS PACE, MTLAPLE 


ee MT_A,P_E., 1,9 


The first READ references the first record, the second READ references the second 
record, and the third READ references the second record on MTAPE. 
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6.8.3. ENDFILE Statement 





Function: 


To record an endfile record on a specified input/output device. 


ENDFILE u 


where: uw is either an integer constant or an integer variable identifying an input/ 
output device. 





(1) Execution of this statement causes creation of a unique endfile record. The form 
of this record depends upon the processor being used. 


(2) When such an endfile record is encountered during execution of a READ statement, 
the action taken depends upon the processor being used. 


Example: 


| ol ror MA TA....,..) 
poop 


poem MT AP, E), 1,0 
Seeenerrree LE, MT A,PLE 
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7. SPECIFICATION STATEMENTS 


GENERAL 


Specification statements are nonexecutable and must precede all other statements in 
a program unit except the FUNCTION, SUBROUTINE, BLOCK DATA (all in sebprogian 
units), and FORMAT statements (see Table 2-2). 


The specification statements are: 


type-statements 
DIMENSION statement 
EQUIVALENCE statement 
COMMON statement 
EXTERNAL statement 


This section describes all but the COMMON and EXTERNAL statements. Because of 
their association with procedure subprograms, these statements are described in 
Section 8. 


TYPE-STATEMENTS 


Function: 


To explicitly declare the data type of a symbolic name and/or declare an array. 


INTEGER vy,v5, ... 


‘and 


REAL V15¥9, - 


and 


DOUBLE PRECISION v,,v5, . 


and 


COMPLEX V1+V9> see 


and 


LOGICAL vj,v5, ... 


where: each vis a variable, an array name, a function name, or an array declarator, 


separated from the next v by a comma. 
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Rules: 


(1) Explicit type declaration of a symbolic name applies to all appearances of that 
symbolic name in the same program unit. 


(2) Any symbolic name of data type double precision, complex, or logical must be 
typed explicitly, since these are not governed by the rules for implied typing. 


(3) Explicit type declaration of an array refers to each of its array elements. In a 
main program unit, the.array declarator may contain only integer constants in the 
subscript. In a subprogram unit, the array declarator may contain integer constants 
and/or integer variables. 


(4) Type-statements of the forms shown may not be used for explicit type declaration 
of a function name in a function subprogram if that name appears in the FUNCTION 
statement. The FUNCTION statement provides for explicit type declaration of 
the function name of the function subprogram. 


Examples: 
(1) 
TNT E/G ER, MAT REX1(2,, 3,,,4,),.X,QiR,T 
This statement types the array MATRIX and the variable XQRT as integer type. 


(This is redundant, since MATRIX is implied integer type.) It also declares the 
dimensions of array MATRIX. Such a statement can appear in any program unit. 


(2) 





This statement can only appear in a procedure subprogram because it.contains an 
adjustable dimension. The function reference or subroutine call that must precede 
such a statement defines the value for N at execution time. 
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& 7.3. DIMENSION STATEMENT 


Function: 


To deciare one or more array(s). 


DIMENSION vy,¥5, ... 


where: each v is an array declarator separated from the next by a comma. 





Rules: 


(1) If any of the arrays is adjustable, it can be declared only in a procedure subprogram. 


(2) An array may be declared in a DIMENSION statement and explicitly typed in a 
type-statement, although a type-statement can accomplish both these functions. 


Examples: 


(L) The statements 


een ae NS 140 ALR LRALY 1, (35, 141,55 ALRLRALY (2 5,16 


LNT EIGER, ,A.RR,A,Y,1 


e can also be written as: 


|. | out M£INS,1,0.N, AR,RAIY,2,(,5,,,6,) 
earaeere BARRY TIC 3 di 5 


(2) The statement 








can only appear in a procedure subprogram because of the adjustable array. 
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7.4. EQUIVALENCE STATEMENT 





Function: 


To permit sharing of the same storage space by two or more entities of the same 
program unit. 


EQUIVALENCE (k,), (ko), ...- 


where: each k, enclosed in parentheses and separated from the next k by a comma, 
is a list having the form: 


ay» 89) ese 28n 


each a is either a variable or an array element (not a dummy argument) with 
only constants as subscript expressions; 
m is equal to or greater than 2. 





Operation: 


The following series of statements causes sharing of storage units: 


DIMENSION V2(7T,.2),K3,2) 


|... | Je,0.m, Pyte,x ¥,1¢,1 21,,3,.),,¢ 
| | jp ouse PRECISION 0(2,2) 
| | fuoer eat, ot 


|. | fe,0,0,1 1¥,A,t,8,N,C.8 CY 1)..,.4,.2.0,1) 7)... 60,02). KO). 2) 4) Cd) 





SHARED SHARED 

1 V1(1,2, 1) 
vi¢(1,1,2)]. 

v1(1,2,3) 


From Figure 7-1, it can be seen that: 


V1(1,2,2) 





Figure 7—1. Effect of EQUIVALENCE Statement 


mg Sharing is accomplished on the basis of storage units. For example, V2(4) shares 
the same storage space as the second storage unit of V1(1,2,1). (See Table 2—4 
for storage unit attribute of the different data types.) 


we An array can be referred tc as a one-dimensional array in the EQUIVALENCE 
statement regardless of the number of dimensions in the array declarator. © 
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(1) Each entity in a given list, k, is assigned the same storage or part of the same 
storage. The sequence of items in a list is unimportant. 


(2) The number of subscript expressions for an array element must be either the 
same number as in the array declarator or it must be 1. However, the number of 
storage units reserved by the array dectarator must not be exceeded by the array 
element reference. 


(3) If a two-storage entity is equivalenced to a one-storage entity, the one-storage 
entity will share storage with the first storage unit of the two-storage entity. 


(4) If one array element is equivalenced to an element of another array, both arrays 
are equivalenced. (See 2.7.4 for the order of the array elements.) 


(5) A dummy argument of a subprogram must not appear in an EQUIVALENCE state- 
ment. 7 


(6) When one entity of a list, k, has its value defined, all its associaced entities of 
the same list are defined. For example, when D(2,1) in Figure 7-1 is assigned a 
value, the contents of K(1,2), K(2,2), L(4), L(5), and C are defined. 


(7) The programmer must avoid contradictions when referring to the same array more 
than once in an EQUIVALENCE statement. For example, the following statement 





causes equivalence of A(3) and D({1), which is a logical contradiction. 


(8) Special considerations apply when a COMMON statement is involved: 


w When two entities are equivalenced, both entities must not appear in COMMON 
statements of the same program unit. 


we Although it is always possible to equivalence an array past its end in an 
EQUIVALENCE statement (as was done in Figure 7-1), it is not always 
possible to extend an array ahead of its beginning (as was done to array DP 
in Figure 7-1) if one of the arrays is listed in a COMMON statement. 














7 
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Examples: ; & 


(1) The primary intent of the EQUIVALENCE statement is conservation of storage. 
For example: 


Paces 





DI MEIN;S ! ON | NT GR (5,00 RiL.(45,0,0 
Lo | pou pie PREC SION, P0504 
—> | le eut iv at eN.C.£ DiP.(.) 1LN.LTG.RIC,1 DP 5,3 
os Deearee © 


Peet 
P| foo 10 


UN. T.GiR(K,) 2 


4 

— 
bral 
io 


reef . 
P| tbo. 210, 211.115.9040 
| 20) 


RL CK = 


a at Wikre M.P,RiLLN,T., 1310 1 HLT.G,Ri, ik oe fg. etnies 
FLOR MjAT 

cewale 

| too, 410, sky =,1,.,5,0,0 
OP, (,K}) = 





The EQUIVALENCE statement cuts down storage requirements by 1000 storage 
units. 


(2) The EQUIVALENCE statement can force the storage of non-COMMON variables 
or arrays. In this example, the variables R, T, B, and A are stored consecutively 
in memory: 





When the statement 


R,E,A.D/( M,C, A,R,0.,.F,) V,0,K 


is encountered, a value for R, T, B, or A is read in depending on whether the 
value of K is 1, 2, 3, or 4, respectively. 
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(3) An EQUIVALENCE statement can make allowances for errors in spelling that 
occur in a program unit. For instance, if an array was declared as CHIEF(3,4,5) 
and the name CHEIF was occasionally used to refer to the array, it is only 
necessary to write 





rather than correct every statement containing CHEIF. 
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8. PROCEDURES AND PROCEDURE 
SUBPROGRAMS 


8.1. GENERAL 


Up to now, programs that contained only one program unit, the main program, have 

been considered, with little discussion of procedures. This section describes procedures 
that can be defined in a program unit and expands the concept of procedure subprograms, 
showing how the programmer can use these and, where they are not supplied with the 
compiler, create his own external procedure subprograms. 


The four categories of procedures and procedure subprograms are: 

® statement functions 

@ intrinsic functions 

@ extemal functions 

@ extemal subroutines ‘ 

A procedure may be a single executable statement or it may be a fixed series of 
statements. The operation(s) to be performed by a given procedure is always the same. 
Each time a procedure is invoked, the values to be used in it can be defined by the 
programmer if required and the procedure can return values. The advantage of using 
procedures is that the series of operations for a given procedure need be defined only 
once, eliminating the necessity of redefining the procedure each time it is invoked. 
There are two methods for passing values to and from procedures: (1) by use of argu- 
ments, (2) by use of the COMMON statement. Either, or both, may be used, except 
that a function reference requires at least one argument. 

A function procedure is invoked by a function reference of the form: 


function name (ay ,a9,.--,4n) 


where function name is the symbolic name that identifies the function procedure, and 
each ais an actual argument. 
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A function reference appears as a primary either in an arithmetic or in a logical expres- 
sion. For each actual argument in the function reference there must be a dummy argument in @ 
the definition of the function procedure that will be replaced by the value of its actual 

argument during execution. Some function procedures, namely intrinsic and basic ex- 

ternal functions, are furnished with the processor and need not be defined by the 

programmer. Statement functions and function subprograms are written by the programmer. 

An example of a basic function reference is: 





The function reference is COS(THETA), and THETA is the actual argument that must 
be defined before the function reference is encountered. During execution of the 
arithmetic assignment statement, THETA is evaluated and a value for COS(THETA) 
is then returned to the arithmetic expression so that it can be computed and a value 
can be assigned to C. A function reference must have at least one actual argument, 
but may have many. For example, a function subprogram is called AREA and computes 
the area of a triangle, given the three sides, and is referenced in an arithmetic 
expression as Y = AREA(X,Y,Z). This subprogram can be expanded to also return a 
value for the perimeter and reference it by Y = AREA(X,Y,Z,PER), where the value 
of PER is used in any statement that follows the function reference. A function 
reference cannot redefine any other values in the same expression containing the 


function reference. & 


A subroutine call invokes a procedure external to the calling program unit, but it 

does not necessarily return any values to the calling program unit and it need not have 
any actual arguments. After the task of the subroutine (defined by the programmer) is 
completed, control is returned to the next executable statement of the calling program 
unit. For example, if the same list of records is read in many times in a program, the 
procedure can be defined once in an external subroutine subprogram and be called 
upon many times in the program by CALL READIN(ARRAY), where ARRAY is the 
actual argument. : 


8.1.1. Statement Functions and Intrinsic Functions 


Statement functions and intrinsic functions are function procedures requiring few 

: computer (machine coded) instructions for execution. Therefore, when the program 
is submitted in machine code form for execution, it is feasible to insert these 
computer instructions inline for each reference. For example, if the intrinsic 
function ABS is referenced many times in a program, the instructions for this 
procedure will be repeated inline in the machine coded executable program as 
shown in Figure 8-1. 
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normal 
instruction 
sequence 


cena for ABS 


instructions for ABS 


instructions for ABS 


Figure 8~]. Inline Coding of Statement Functions and Intrinsic Functions 


8.1.2. External Procedure Subprograms 


External procedure subprograms (function and subroutine subprograms) are function 

procedures that generally require more than one FORTRAN statement for definition 

and many computer instructions for execution. Therefore, the set of instructions is 

recorded only once in the machine coded program and placed out-of-line so that it 

can be entered only by a machine coded control instruction. After execution; control 

is returned to the referencing program by another machine coded control instruction. 
For instance, if the cosine function COS is referenced many times in a program, the 

@ machine coded program will appear as in Figure 8—2. ; 


last instruction to be executed 


l data 


cosine procedure 





Figure 8~2. Out-of-Line Coding of External Procedure Subprograms 
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Communication Between Program Units os 


Values may be transmitted to a procedure or a procedure subprogram through the ® 
actual arguments of the reference or call and returned through a function name and 

actual arguments. The COMMON statement is another method of transmitting 

information between different program units. The EXTERNAL statement identifies 

a symbolic name used in the list of actual arguments as the name of an external 

subprogram. 


‘Valid Forms of Arguments 


Table 8--1 lists, for reference purposes, the valid forms of procedure arguments. 
These are explained, in more detail, in the applicable paragraphs. 


FORM OF ARGUMENT(S) 
PEEVES ACTUAL DUMMY 


statement function arithmetic expression, logical variable 
expression 


basic extemal arithmetic expression of required 
function type 


function subprogram@) arithmetic expression, logical variable, array 


: expression, atfay name name, name of 
external procedure 





subroutine sibprderam) adthmetic expression, logical variable, array 
expression, array name, name name, name of 
of extemal procedure, Hollerith external procedure 
constant 


intrinsic function arithmetic expression of 
required type 


@ if an actual argument corresponds to a dummy argument that is defined or rede- 


fined in the subprogram, the actual argument must be a variable, array element, 
or array name. 


Table 8—]. Forms of Argument 











8.2. 


8.2.1. 












8 


SECTION: 





FUNDAMENTALS OF FORTRAN 


STATEMENT FUNCTION 


There are two types of statement functions: arithmetic statement function and logical 
statement function. 


Arithmetic Statement Function 
Function: 


To define an arithmetic procedure with one FORTRAN statement. This procedure 
may be used in its program unit as many times as required. 


function name(a,,a » «-,@ )= Hmited arithmetic expression 
n 


where: function name is a symbolic name identifying the procedure. 


the a’s comprise the list of one or more dummy arguments, a list of variables 
enclosed in parentheses and separated by commas, 


limited arithmetic expression is similar to an arithmetic expression except that it 
may not contain an array element as a primary, and it may not reference, in any 
way, the statement function of which it is a part. 


the maximum value of n is defined for each processor; there is no standard 
value. 





Rules: 


(1) An arithmetic statement function reference is a primary in an arithmetic expression. 
The actuai arguments are arithmetic expressions and must correspond in number; order, 
and type with the dummy arguments. 


For example, the statement function 





The value for the actual argument R is substituted for its dummy argument A, S for B, 
and T for C. A value is returned to the statement function reference so that the arith- 
metic expression Y — AVRGE(R,S,T) can be evaluated and assigned to Z. The state- 
ment function reference 
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(3) 


(4) 


(S) 


(6) 


(7) 
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is incorrect because the data type of N is inconsistent with the data type of C. The 
statement function reference 





is incorrect because there are only two actual arguments, and three are required by 
the statement function definition. 


All arithmetic statement functions must appear between the specification statements 
(if any) and the first executable statement of the program unit (see Table 2—2). The 
limited arithmetic expression may use references to previously defined statement 
functions as primaries. 


The function name may not be used as a variable within the same program unit; it is 
used as part of a statement function reference. 


Each dummy argument in the list must appear at least once in the limited arithmetic 
expression. 


The arithmetic statement function definition is a nonexecutable statement and cannot 
have a statement label. 


The function name is governed by the mles for imptied and explicit typing of symbolic 
names; it need not have the same data type as its arguments. Evaluation of the limited 
arithmetic expression is subject to the rules for mixed type evaluation of arithmetic 
expressions, and assignment of a value to the function reference is subject to the 
tules of Table 4-1. 


Variable names used as dummy arguments are purely local to the arithmetic statement 


_ definition. They may be used elsewhere in the same program unit for any purpose. 
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Examples: 


(1) 


peers G? NS. 1 ,0,N) ,ARIRIALY 
| | lsusmsja Ay Bi =| ,Ai*),* 12; ,+, ,8,*|* + C.*,*)2 


Ps | [bo 10, at3 
10! layer aly N =, JN 


$,UM.S.Q CAR, RAY (1), AR RA VIC 2) ARRAY, (3) 1), 


ime 
N 
" 


A value of 14.0 is assigned to Z. Note that the mixed data types in the limited 
arithmetic expression conform to the rules listed in Tables 3-1 and 3-2. 


(2) The same statement function can be defined with only variables in its limited arithme- 
tic expression, but these variables must be defined prior to the reference. For example: 


DI MEIN, S,1,0.N, JA, RR ALY, (3 


reaene MS 1O.(,A.,,8,, .C =] ,A** KT, +, BI** K,2, 1+, .C.*.*1K 3 
id galehile 


eS irae 
Shien g | Hees = 
wae Me | 
a an N =, JN : 
ret 


Z, =, |S U.MS,.Q( ARR Y,(.1,),,,4,R,R.A, %1(.2.), ARRAY. ()3.).) 


Nm TN I 


(3) A statement function definition may contain references to previously defined statement 
functions. For example: 


P| we ura exuty 

| sum sia(.a,.B) =. 44 412 +, ,8,%,*,2, .+] ,S,U,M,(,A,, ,B,) *)*,2 
Foeee 

Peres : 


-,0 
$UMS,Q(,Z,+,2,.0,,.)4,. 0), +, 2 


ae 
QO ft 
0 i 


The variable C will be evaluated as the real type approximation of 5.07 + 4.07 + 9.02 
+ 3.0 or 125.0. 
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8.2.2. Logical Statement Function 
Function: 


To define a logical procedure with one FORTRAN statement that may be referenced 
in its program unit as many times as required. 


function name(a,,a9, +++ 4Q) = limited logical expression 


where: function name is a symbolic name, explicitly typed as logical, that 
identifies the procedure. 


the a’s comprise the list of one or more dummy arguments, a list of 
variables enclosed in parentheses and separated by commas. 


limtted logical expression is similat to a logical expression except 

that it may not contain an array element as a primary, and it may not 
contain a reference, directly or indirectly, to the statement function 

of which it is a part. 


the maximum value of n is defined for each processor; there is no 
standard value. 





Rules: 


(1) A logical statement function reference is a primary in a logical expression. 
The actual arguments may be arithmetic or logical expressions and must 
correspond in number, order, and type with their corresponding dummy arguments. 


(2) All logical statement functions must appear between the specification state- 
ments and the first executable statement of a program unit. 


(3) The function name may not be used as a variable within the same program 
unit; it may only appear as part of a statement function reference. 


(4) Each dummy argument in the list must appear at least once in the limited 
logical expression. 


(5) The logical statement function definition is a nonexecutable statement and 
cannot have a statement label. 


(6) Evaluation of the limited logical expression is subject to the rules for evalu- 
ation of a logical expression (see 3.4), and it assigns a value of either true 
or false to the statement function reference. 


‘(7) Variable names used as dummy arguments are purely local to the statement 
function and may be used elsewhere in the same program unit for any purpose. 
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Examples: 


q1) 


a | Ieee. teat N,0,0,0 
es | bi op.ILK = Ki/i 21) 812) pay Ki) ye ING EL 10 — 
an 
eer 
beh eee 1.0.0.0 G10, ,T.0 


In the logical IF, if I is an odd integer, the GO TO will be executed (see 3.2.4, 
example 2). 


y y 
e 
~ O 
o 
18 OH 
* Oo 
it 
’ 


‘ 
In the logical IF, if M is a positive odd integer, the GO TO will be executed. 


8.3. INTRINSIC FUNCTIONS 


Intrinsic functions (built-in functions) are provided with the processor and are not 
written or modified by the programmer. A list of standard FORTRAN intrinsic functions 
is presented in Table 8—2. 
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_. “| NUMBER TYPE OF @ 
INTRINSIC FUNCTION DEFINITION OF SYMBOLIC 
ARGUMENTS NAME ARGUMENT | FUNCTION 


absolute value ja| 1 ABS real real 








IABS integer integer 
DABS double double 
truncation sign ofa times largest 1 AINT real real 
integer<| e| INT real integer 
IDINT double integer 
remaindering* a (mod a5) , 2 AMOD real real 
MOD integer integer 
choosing largest value | max (¢,,.a5, wat) 22 AMAX0O integer real 
AMAX1 real real 
MAX0 integer integer 
MAX1 real integer 
DMAX1 double double 
choosing smallest value} min (a, Qo; sae) 22 AMINO integer real 
AMIN1 real real 
MINO integer integer 
MIN1 real integer 
DMIN1 double double 
float conversion from integer 1 FLOAT integer real 
to real 
fix conversion from real to 1 IF IX real integer 
integer 
transfer of sign sign of a, times 1a4| 2 SIGN real real 
ISIGN integer integer 
DSIGN double double 
Positive difference ie min(a, 99) . 2 DIM real real 
{DIM integer integer 
obtain most significant 1 SNGL double real 
part of double precision 
argument 
obtain real part of 1 REAL complex real 
complex argument 
Obtain imaginary part 1 AIMAG complex real 
of complex argument 
express single precision 1 DOBLE real double 
argument in double 
Precision form 
express two real argu- | a, +a5\/-1 2 CMPLX real complex 
ments in complex form 
obtain conjugate of a 1 CONJG complex complex 


compiex argument 


*The function MOD or AMOD (a, 85) is defined as a —[xla5, 
not exceed the magnitude of 9 /ay and whose sign is the same as a, /a9 (see example 3). @& 


where[x] is the greatest integer whose magnitude does 





Table 8=2. Intrinsic Functions 
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(1) An intrinsic function is referenced as a primary in an arithmetic or relational expression, 
by name and list of actual arguments. The actual arguments must agree in type, number, 
and order with the specifications of Table 8—2 and may be any expression of the speci- 
fied type. Note, however, that the number of arguments for the MAX and MIN intrinsic 
functions is variable. 


(2) The intrinsic functions AMOD, MOD, SIGN, ISIGN, and DSIGN are not defined when the 
value of the second argument is zero. 


(3) It is not necessary to declare the type of an intrinsic function in the program unit that 
contains a reference to an intrinsic function. The data type is already known to the 
processor. 


(4) Fora valid intrinsic function reference in a program unit, the symbolic name: 


(a) must appear as specified in Table 8-2, followed by the list of actual arguments, 
in parentheses, also specified in Table 8-2; 


(b) must not appear in an EXTERNAL statement or be used as a variable or array 
name in a program unit where it appears as a reference; 


(c) must not appear in any type declaration different from the implied type decliara- 
tion of Table 8—2. 


(5) If a particular intrinsic function is not referenced in a program unit, its symbolic 
name may be used for any valid purpose in that program unit. 


Examples: 


(1) This series of statements reads in the six complex type elements of array VCTR, 
then prints the conjugate of each array element, using intrinsic functions CMPLX, 
REAL, and AIMAG. 2 


pf foowet et VC TLR. 61) 
pe eb (geting) 


|, , { {po 2,0 Kel, 6 
[zl Wve eick PLAX Cf 
haa | ees di ee 
3.0! lromwast 
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This example uses the MAX0 intrinsic function. 


eee oe 
P| Tt CM A.X.0,0,1, 0.5.5.6.) 1-5 4, 0) 1,0,0,,,2;0,0,,,3,0,0 


(3) 


(4) 





In this case, the maximum value is 40, causing the arithmetic IF to jump to the 
statement labeled 200. 


The remaindering function (also called the modulo or residue function) is useful in 
modulo-arithmetic. This intrinsic function divides one argument by the other and 
retains only the remainder prefixed by the sign of the quotient. For example, the 
reference MOD (9,8) returns a value of 1, and the reference AMOD (1.22, 1.1) 
returns a value of 0.12. In the statement 





a jump is made to the statement labeled 20 only if J is zero or an exact multiple of 
K (K can be multiplied by an optionally signed integer to obtain J). 


The IFIX and FLOAT can be used to calculate the values of arithmetic expressions 
where certain combinations of data types are not defined in standard FORTRAN (see 
Table 3—3). For example: 


The FLOAT function uses the real type form of K for the multiplication. The IFIX 
converts the real type result of the multiplication to integer type form so the sub- 
traction can be performed. : 





12 
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8.4. RETURN STATEMENT 
Function: 


To return program control from a function or subroutine subprogram to the program 
unit that referenced or called the subprogram. 


RETURN 


Rules: 


(1) There must be at least one RETURN statement in every function and subroutine 
subprogram. 


(2) When executed, the RETURN statement terminates further execution of the 
subprogram that contains it. 


Examples: 
See 8.5.2.2 and 8.6.3. 

8.5. EXTERNAL FUNCTIONS 
External functions are of two types: external subprograms written by the programmer 
in FORTRAN, and basic functions. Basic functions are external function procedures 
supplied with the processor and stored in auxiliary storage in non-FORTRAN repre- 
sentation. These basic functions do not require compilation and can usually also be 
utilized by non-FORTRAN users of the processor. 


8.5.1. Basic External Functions 


_ Table 8-3 lists standard basic external functions. 
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NUMBER ween & 
BASIC EXTERNAL of SYMBOLIC ; 
FUNCTION DEFINITION ARGUMENTS**| NAME ARGUMENT FUNCTION 





exponential e@ 1 EXP real real 
1 DEXP double double 
1 CEXP complex complex 
natural logarithm tog _(@) 1 ALOG real real 
: 1 DLOG double double 
3 1 CLOG complex complex 
common logarithm log, g(a) 1 ALOG10 real real 
— 1 DLOG10 double double 
trigonometric sine sin(a) 1 SIN real real 
1 DSIN double double 
1 CSIN complex complex 
trigonometric cosine cos (a) 1 cos real real 
1 ocos double double 
1 ccos complex complex 
hyperbolic tangent tanh (a) 1 TANH real teal 
square root (ay! /2 1 SQRT real “real 
1 OSQRT doubie double 
1 CSQRT complex complex 
arctangent arctan (a) 1 ATAN real real 
1 DATAN double double 
arctan (a, (Ag) 2 ATAN2 real real 
2 DATAN2 double double : 





remaindering* @, (mod a5) 2 DMOD double double 


modutus Vireat part)“+ (imaginary part) i CABS complex real 


*The function OMOD (a, ,a, ) is defined as a, ~ (xa, where [x] is the largest integer whose magnitude does not 
exceed the magnitude of ay /& and whose sign is the.same as the sign of a, /ap. 


**All angles are expressed in radians. 


Table 8-3. Basic External Functions 


Rules: 
(1) A basic external function is referenced as a primary in an arithmetic expression. 
(2) The form of the reference is: 


function name (arg, ) 


or 


function name (arg, arg) 
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where: function name corresponds to one of the function names listed in 
Table 8-3. 


each a is a dummy argument, the list is enclosed in parentheses, and the 
a’s are separated by commas. 


(3) It is not necessary to explicitly type a basic function name in the referencing 
program unit. 


Examples: 


(1) The following sequence finds the sine of 30 degrees (which is 0.5): 


| | jpouste PRECISION, PI. Y 
Ps | a taL Pit /. 3.51 4.159 D. 0 
See A 

eel 

ee Lee es OS mPa t/a 10 

| weer M, PRINT 


(2) This example shows a basic external function reference containing another 
fuaction reference in the list of actual arguments: 


(3) This example shows a basic external function reference in a logical expression: 





8.5.2. Function Subprograms 


Function subprograms permit the programmer to create external arithmetic and logical 
procedures and refer to these many times in-the body of the program using function 
references. The first statement of the external function subprogram definition must 
be a FUNCTION statement, and the last line must be an END line; between these 
there must be at least one RETURN statement and a definition of the function name. 
Each reference to a subprogram in an arithmetic or logical expression will have a 
value returned to it at the point of reference. Additional values may be returned for 
us@ in statements that follow the reference. 
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FUNCTION Statement 





Function: 


To identify an external function subprogram. 


t FUNCTION function name (41,45, +++ aq) 


where: tis explicit type declaration of the value to be returned in the function 
teference and may be omitted, in which case the implied type of 
function name determines the data type of the value. 
function name is a symbolic name identifying the function subprogram; 
each a is a dummy argument; the list is enclosed in parentheses, and the 
a’s are separated by commas. ; 
n must be at least 1, and its maximum value depends upon the processor 
being used. 





Rules: 


(1) If t, the explicit type declaration, is present, it must be INTEGER, REAL, 
DOUBLE PRECISION, COMPLEX, or LOGICAL. 


(2) Each dummy argument must be an external procedure name (see EXTERNAL 
statement), a variable, or an array name. 


(3) The function name used in the FUNCTION statement must not appear in any . 
other nonexecutable statement of the function subprogram. 


(4) The symbolic names of the dummy arguments must not appear in an 
EQUIVALENCE, COMMON, or DATA statement of the function subprogram. 


For examples, see 8.5.2.2. 

Function Subprogram Definition 

Function: 

To define the procedure that computes a result(s) returned to the referencing 


program unit. 


FUNCTION statement 


RETURN 


END 
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Rules: 


6 (1) 


(2) 


(3) 


(4) 


(S) 
(6) 


(7) 


All rules applicable to the FUNCTION statement are applicable to the external 
function definition. 


The function name of the subprogram must appear as a variable at least once 

in the subprogram. During every execution of the subprogram, this variable must 
be defined before it may be referenced or redefined. The value of the variable 
at the time of execution of any RETURN statement in this subprogram is the 
value of the function and is the value returned to the function reference. 


The subprogram may define and redefine one or more of its arguments so as to 
effectively return results in addition to the value of the function. 


The function subprogram may contain any statements except BLOCK DATA, 
SUBROUTINE, another FUNCTION statement, or any statement that directly 
or indirectly references the subprogram being defined. 


The function subprogram must ¢ontain at least one RETURN statement. 


If a function reference causes a dummy argument in the referenced function to 
become associated with another dummy argument’ in the same function or with 
an entity in common (see COMMON statement), a definition of either within the 
function is prohibited. An example of such a function reference is: 

Y = ADD(A,A) 


The following rules apply to arguments involving arrays or array elements: 


(a) If an actual argument is an array element, its dummy argument must be 
either a variable or an array name. 


(b) If an actual argument is an array name, its dummy argument must be an 
array name and that array must be declared in the subprogram with a 
size (in elements) that does not exceed the actual argument array. 


(c) If the actual argument is the xth element of an array containing z elements 
and the dummy argument is an array name, that array must be declared in 
the subprogram with a size that does not exceed z — x tl elements. 


(d) A dummy array declarator may use one, two, or three subscript expressions, 
tegardless of how the actual array was declared. Each subscript expression 
of the dummy array declarator may be either an integer constant or an integer 
variable. If any subscript expression of the dummy array declarator is an 
integer variable, that array is called an adjustable array. 


An adjustable array declarator must have each of its integer variable sub- 
script expressions listed as dummy arguments and each must be defined by 
its actual argument. These variables must not be redefined in the subprogram. 
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Examples: @ 


(1) This subprogram calculates the sum of the squares of all positive odd integers 
from 1 through the given positive integer K: 


SM, =] 0 
po | teswriem ve, at 
Lack Saas ilbrg 1/0, .N,=,1 K 
gems te Pd Wit = - KS .W TIC.H 
Pea ol lice KIS, W.T.C.H. .. .E,Q. H SM, ,=, SIM, ,+, N.*.* 42 
| | leet URN ; 
Po | leven 


Thus, if this externa! function subprogram is referenced by 


fe apes a ad 


elie ER M,Q,D,S,Q e 
Perren 


Peeeee 
fees Gece 
asd 





5, +, ,5.M,(,N,3,) t 


the actual argument N3, which has a value of 5, will be substituted for its 
dummy argument K in the subprogram. The value of SM is 1° + 3“ + 5° = 35. 
When the RETURN statement is reached, this value is substituted for the 
function reference, so that the value assigned to M is 40. 


u 
uy 


(2) This function subprogram, which has more than one RETURN statement, 
calculates the absolute sum of all array element values in a 100-element 


array: 
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© | | de UN.ciT,1,0.N S,MA,R, AVY, ,(,Z,R,4,Y,) 
Perea her creer 


D.O 0, Ki=114, 131010 
S MARIA, Y=,S MA RAY +/Z, RAY, (K 
Ree alee S MARA, Y¥,.,G,E, 10). 10 R, & TURIN, 
po | Ia RAYic, 1aiSMAARALY 
i | tee, 7.upRN 
jovi | leno 


@ Serere 0a = S M,A.RA Y,(, AIR, RA, Y,) +, 8S MJA,R,A,Y,( BRAY 


(3) This function subprogram also returns a value through one of its dummy argu- 
ments. Given the three sides of a triangle, it calculates the area and returns a 
value for the perimeter. 


eee ALR.E.A.(\A,, By, Ci, :P.E.R 


PER, j=, A, +, ,By +, IC 
S. =, IPJER/ 2. 10 

fcaee tela REAL =, $.Q.2,7.(S,*1(,$,-,A) .*,( 8 .-,Bi) .* (48 ,=,C,) ,) l 
R,E,T,UIR.N 


D 
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The referencing program contains: 


|| ieee L. = ARE A(X, Y,,.2,, PER! ) 


Note that, in the referencing program, one value for AREA will be returned in 
the function reference. Another value will be returned as the value of PERIM. 
This value must not be used in the same statement as the function reference, 
but can be used in statements that are executed after the function reference. 

If a value is to be returned to an actual argument in the function reference, that 
actual argument must not be a constant; this would be an attempt to redefine a 
constant. , 


The following subprogram illustrates the use of an adjustable array. The 
subprogram calculates the sum of the elements in an array, but each time 
the subprogram is referenced, the dummy array may have different dimensions. 


rT | fT J 
wy 
io 
= 
> 
" u 
- = 
oy = 
Es 
> 
cs 
> 
Fo] 
7 
> 
=< 
Fad 
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The referencing program unit may contain: 


| lo MeNS.15.N $F. T.110.4,, 51236 SjE 7.2.01 ,1 SE, 7,3,0,8,.,1,9,) 


Renee i 
i | lara =; 18.1.6 $,E.7,3.,14:°. 5) 916 4,. 0, *, 5] .,0, *,6,-,9)) 
feel | 


P| davies 
haart 


$,1,GM,A,()S,E,73,, ,8,°%,1,0,); ./, ,€,8,-,0,*%,1,.;0,) 


Hi 
i 
~ 
wo 


Note that this last reference uses an array element as the actual argument. 
This array element is the 119th element (see Table 2—6) in an array that can 
contain a maximum of 120 elements. Hence, the function reference is obtaining 
the sum of the last two array elements of array SET1. 


& (5) This example uses the Newton-Raphson method to find the cube root of a 
number, y, as follows: 


(a) Make an initial guess at the cube root of y. (In this program, y is divided 
by 3 if its absolute value is greater than or equal to 1; otherwise, it is 


multiplied by 3.) This value is called xy: 


(b) Substitute x1 in the expression 


1 
3 (2x1 + y/x 7) 


to arrive at a new value, Xoe 


(c) Compare x9 with X4- If the comparison satisfies a criterion set up by the 
programmer, x, is the cube root of y; if not, X» becomes x,, and this new 
Xy is substituted in the above expression to arrive at a new value for x5. 
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Steps (b) and (c) are repeated until the comparison of X» with % is satisfied. 


0 BR. OO.Ti= 0,* 


ree ie C}B,.R,0,0,.T,. .£ Q,. xX 
Po | e.c.pik.o.0.7 
mace 


G,O, ,Tj90, .1,0 
END 


A reference to this subprogram could be: 


Using this subprogram, the following is a list of successive CBROOT’s calculated 
in finding the cube root of 1.0, 0.9, 0.1, and -0.027. The solution is repeated 









; 8 
SECTION: 


R,E, TURIN 





several times because, for this particular processor, the arithmetic operations 


are performed internally with greater precision than was required for the output. 


cube root of 1.000 
3.2222223 
2.1802527 
1.5236256 
1.1593397 
1.0208964 
1.0004248 
1.0000002 
1.0000000 


cube root of .900 
1.8411523 
1.3159345 
1.05053145 
.9721879 
.9655354 
-9654894 
-9654894 
.9654894 





cube root of .100 
.5703704 
.4827094 
.4648626 
.4641599 
.4641589 
.4641589 


cube root of -.027 
~1.4257421 
-.9549222 
-.6464846 
-.4525238 


* -.3456326 


-.3057596 
-.3001078 
-.3000000 
-.3000000 
~.3000000 
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Such procedures can be refined by more accurate first estimates, the use of 
double precision, and different exits from the loop. For example, the precision 
of the result could be limited to a satisfactory value with the statement: 





where: DELTA is an additional argument. 


8.5.2.3. References to Function Subprograms 


Function: 


To obtain a single value for use in an expression, by reference to a function 
subprogram. 


name(ay ao, +117) 


where: name is the programmer-written symbolic name of the function. 


each a is an actual argument; the list is enclosed in parentheses, and 
the a’s are separated by commas. 
a has a minimum value of 1. 





Rules: 


(1) 


(2) 


(3) 


(4) 


(S) 


(6) 


Each actual argument may be a logical or arithmetic expression, an array 
name, or the name of an external procedure (see EXTERNAL statement). 


If a value is to be returned through an actual argument, that actual argument 
need not have its value defined prior to the reference, and it must not be a 
constant. . 


The actual arguments must agree in order, number, and type with the dummy 
arguments of the subprogram. 


If the actual argument is an array element, its dummy argument may be either 
an array name or a variable; if the actual argument is an array name, its 
dummy argument must be an array name. 


If an actual argument is the name of an external procedure, its dummy argument 
must be used in the subprogram as the name of a function or subroutine external — qe 
to the procedure. 


The operation of a RETURN statement in the function subprogram must 
intervene between each reference to the same function subprogram. 


Examples: 


See FUNCTION, EXTERNAL, and COMMON statements. 
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8.6. SUBROUTINE SUBPROGRAMS @ 


A subroutine subprogram is an external subprogram defined by FORTRAN statements, 
starting with the SUBROUTINE statement, and is invoked by a CALL statement. 


8.6.1. CALL Statement 


Function: 


To invoke a subroutine subprogram. 


CALL name 
or 


CALL name (4 ,ay, +) 


where: name is the symbolic name of the subroutine subprogram. 
each a is an actual argument; the list is enclosed in parentheses, and the 
a’s are separated by commas. 





Rules: 


(1) Each actual argument may be an expression (logical or arithmetic), an array . 
name, a Hollerith constant, or the name of an external procedure (see EXTERNAL 
statement). @ 





(2) The actual arguments must agree in order, number, and (except for a-Hollerith 
constant) type with the corresponding dummy arguments of the subroutine. 


(3). If an actual argument corresponds to a dummy argument that is defined or rede- 
fined in the subroutine, the actual argument must be a variable, array element, 
or array name. 


(4) If an actual argument is an external function name or subroutine name, the 
corresponding dummy argument must be used as an external function name or 
subroutine name, respectively. 


(5) If an actual argument is an array element, its corresponding dummy argument 
must be either a variable or array name; if an actual argument is an array name, 
its corresponding dummy argument must be an array name. 


(6) If a subroutine reference causes a dummy argument in the referenced subroutine 
to become associated with another dummy argument in the same subroutine or 
with an entity in common (see COMMON statement), a definition of either entity 
is prohibited in the subroutine. An example of such a reference is: CALL ADD(A, A). 


(7) Between any two successive calls on the same subroutine there must be the 
operation of a RETURN statement in the subroutine. 
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8.6.2. SUBROUTINE Statement 


8.6.3. 


Function: 


To identify a subroutine subprogram. 


SUBROUTINE name 


or 


SUBROUTINE name (a1 ,a9, +144) 


where: name is a symbolic name that identifies the subprogram. 
each a is a dummy argument; the list is enclosed in parentheses, and the 
a’s are separated by commas. 





Each dummy argument is a variable, an array name, or an external procedure name 
(see EXTERNAL statement). 


Subroutine Definition 
Function: 
To completely define the step-by-step procedure for a subroutine. 


SUBROUTINE statement 


RETURN 





Rules: : 

(1) The first statement in the subroutine must be the SUBROUTINE statement, and 
the last line must be the END line. The subroutine definition must contain at 
least one RETURN statement. 


(2) The symbolic name of the subroutine must not appear in any statement of the 
subroutine except the first. 


(3) The symbolic names of the dummy arguments may not appear in an EQUIVALENCE, 
COMMON, or DATA statement of the subroutine. 
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(4) The subroutine may contain any statements except BLOCK DATA, FUNCTION, 
another SUBROUTINE statement, or any statement that directly or indirectly calls & 








this subroutine. 


(5) The subroutine may define or redefine one or more of its arguments so as to return 
results through its arguments. If the subroutine returns a result through a dummy 
argument, its actual argument must not be a constant. 


(6) Adjustable arrays may be used in the subroutine, in which case any adjustable 
dimensions may be passed as actual and dummy arguments for use in the array 


declarator of the subroutine. 


(7) The RETURN statement returns program control to the next executable statement 
following the CALL statement in the calling program unit. 


Examples: 


(1) This subprogram clears the elements of any array to 0.0. 


Ir acketcdels B.RiO.U, TINE LiEIARi (AIR RIAL, NI) ! pili bout 
Di MEN Si ON, AR RALY(N 


Foon WO, K =,1 N 
Pee 


ALRiALAL FT “ 


Peete 


A calling program unit contains: 





P45 ewes | ead LE ALR A,RALY, 14:7, 51 * 16 


papa CLEA, RC AR ALY (.1,,,1,,,1,),,,1{0,0,) : ‘ 
ae cial | | mats 
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The first three subroutine calls will clear all of the referenced artays; the last 
subroutine call will clear only the first 100 elements of the array ARAY. 


The following subroutine subprogram, UPDATE, computes the current value of 
a savings bank account given a starting value, V, the number of interest periods, 
M, and the interest rate per interest period, R. 


Pi 0] |v.2,¥, +182 
Po | ee tue ! 
Sen 


| au 


This subroutine subprogram is invoked in the following program. The program 
computes the current value of an account based on the assumption that the 
interest is computed and accrued at the end of each three-month period from 
the month of deposit. 
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The first card of the data deck is organized as follows: 





columns 1 and 2 — the current month (1-12), right-justified 
column 3 - blank character 
columns 4 through 7 — current year 


Each succeeding data card is organized as follows: 


column 1 — blank character 

columns 2 through 7 — account number, right-justified 

column 8 — blank character | 

columns 9 vsueh 48 — name of depositor, left-justified 

column 49 — blank character 

columns 50, 51 — month of last update or (if no previous update) 


month of deposit, right-justified 


hyphen (-) or minus character 


column 52 


year of last update or (if never updated) year 


of deposit @ 


column 57 — $ (dollar sign) 


columns 53 through 56 





columns 58 through 66 - dddddd..dd, the tast updated value or original 
deposit, where each d is a decimal digit 


columns 67, 68, 69 — blank characters 


column 70 — months (0, 1, 2, or 3) from last update or driginal 
deposit to end of next interest period. (If original 
deposit and no update, it must be 3.) 


The last card of the data deck contains all blank characters, except for column 
70, which contains a digit greater than 3. 


The subroutine is invoked by the following program which prints the updated 
record and punches an updated card record concurrently (assuming output device 
5 is a card punch). 
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ae err CY BAIR, -, VE AR) +0 .CiM,O.N, TH, MON, TIH,) + (3 .- RMN TINS 
Mo =, IMT E 


dol THEOL TDI Te EL CMON TH CVE AR WY Ab L EL RMN Tu 
P| wou tie Te 
P| tous ts 


(3) The following simple example shows how a Hollerith constant might be used in 
a CALL statement. The printed line will print LOW, AVRGE, or HIGH, depending 
upon the income read from an input card. 


Ot .M, EIN,S,1,0.N, NAM 


| tol feo, RM ALT 8.4.5, 1,0, 11,0 
pat ea PRINT MAMIE. Sst b8 
r [toe 
30] lh eC i COME 2.0000), 4.0,,.4,0).,5,0, | 
| 40] lepaststy JP RIN TUNA WME, , 5 4,A,¥,R,GE 
ss 


Pare eer PRINT CNAME. CLASS 
| | fot MN, 5,1 0,6 NAM, E}(8,) 


| lwttye.0.3.,1,0,) NAIM E,, CLASS 
|. tol feo, MALT G1 xy, 8, AS i, 15. %..1, 9,4 SiS, ,0,F, ,1,N,C,O,ME, .1,8, .,,A)5,), | 


R RiN 
P| fewo 
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8.7. EXTERNAL STATEMENT 





Function: 


To indicate that a symbolic name (with no arguments) in a list of actual-arguments is 
the name of a subroutine or external function. 


EXTERNAL v,,v., ...,v 
1 n 


where: each v is the symbolic name of an external procedure. 





(1) If an external procedure name is used as an actual argument in 2 function reference 
or subroutine call, the name must appear in an EXTERNAL statement of the program 
unit that contains the function reference or subroutine call. 


(2) If a symbolic name appears in an EXTERNAL statement of a program unit, it can 
appear in a type-statement of the same program unit only if that name is also being 
used as the function name in a function reference within the same program unit. 


(3) If an actual argument is the name of an external procedure, its corresponding dummy 
argument must be used as the name of an external procedure. 


(4) If a symbolic name is being used as the name of an intrinsic function or a statement 
function within a program unit, it must not appear in an EXTERNAL statement of é 
the same program unit. 

(5) The EXTERNAL statement is a specification statement; its position in a program 
unit is shown in Table 2-2. 


Examples: 


(1) The following subroutine can be used to print a table of either sine values or 
cosine values for tenths of a degree, from 0.0 degrees to 90.0 degrees. 


rt] 
(oJ 
(Cae S) 
h 
1 
ld 
yy 
+ 
> 


TABLE 4 ¢ K) A(TABLE2¢(K) (5°, ON EDGR) 


| fwett163,.,2,8) 
| tweet recs 3.0.), ¢( TIABLE2.(.%.), Tag 4 Xl: 2.9 
. 3.0l [por wa.r VX Fi 4 itt 0% ,0, 1,4. 47 

ps | every ipa 

| | leno 
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The calling program contains: 


The first call causes printing of double precision sine values; the second, double 
precision cosine values. 


(2) In the following program unit: 





The use of the EXTERNAL statement is not warranted. It is not the function name 


ABC that is being passed as an argument, but the value of the function reference 
‘ ABC(D). 





31 


UP-7536 
Rev. 1 











8 32 


SECTION: 






FUNDAMENTALS OF FORTRAN wages 
8.8. COMMON STATEMENT 
Function: 


To permit communication between program units without the use of arguments. 


COMMON /Xy [ay /%q/8q/ +6 /X,/ 8, 


where: each x is either a symbolic name, called a block name, or empty (no 
character or one or more blank characters), 
each a is a nonempty list of variables, array names, and/or array declarators, 
with each item in the list separated from the next by a comma. 
If oT is empty, the first pair of slashes before a is optional. 





Operation: 


The COMMON statements of a program unit cause the processor to set aside locations 
for the listed variables and arrays in an area of main storage called the common area. 
These locations can be shared by the different program units that make up a program. 
This common area is composed of blank common and labeled common blocks, If an x 
in the COMMON statement is empty, its associated variables and arrays are stored in 
the blank common block. If an x is a block name, its associated variables and arrays 

2 are stored in that block of the common azea identified by the block name. The total 
size (in storage units) of blank common and each labeled common block is determined 
by COMMON and EQUIVALENCE statements. The order of the locations in blank 
common and in each labeled common block is determined by the order of their appearance 
in COMMON statements. 


In the following sequence: 


| | [common AT, a2, a3, ARC2,2)/EBLT/ BROS), BIS C3, CR( 2) 


Fg tah 7 ter 2 7.0.02 2012) EBL, ZEB 3) 
P| [common S741, WR 3) LBL 27 XR. 2 
| {Ot mE H,S, 10M, X03) 
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eS . the common area will consist of: 


BLANK COMMON 
aI 2] 8 beet]? ste] " bee 


BLOCK LBL1 <j 





BLOCK LBL2 , ~~ 





An EQUIVALENCE statement can lengthen a block (blank or labeled); the only 
lengthening permitted is that which extends a block past the last assignment for 
that block made directly by a COMMON statement. For example, if the following 
sequence is added to the previous sequence, 


| te.e.u.t iv aA ene, OIR Cy 124, 37 Qi Ri 2 Xi 2, 1T R22 
eee 


0.) MEIN, S10 QR8 TRO4 


the labeled common block LBL2 is lengthened by two units. 





EQUIVALENCE 
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are illegal, since this is an attempt to extend a common block ahead of its. first 
assignment for that block, DR. & 
In each program unit, the size (in storage units) of a labeled block with the same 

block name must be the same, if it is present in a COMMON statement. The size 

of blank common need not be the same in all storage units. The total size of blank 

common is equal to the largest number of storage units assigned to blank common 

by a program unit of the program. However, association of a common variable or 

array is by storage unit in a labeled block of the same name or in blank common. 

The nth storage of blank common in one storage unit is shared with the ath storage 

unit of blank common in any other program unit. The ath storage unit of a labeled 

common block with a given name in one program unit is shared with the ath storage 

unit of a labeled common block bearing the same block name in any other program 

unit. 


For example, if one program unit contains 






COMMON, /,TABL ETWN/, A,B, C., 0 /,/)K,1,,,K,2,,,K,3, | 


and a second program unit contains 


oe Oi:Mi: MIOUN M AB i Vv d 


the correspondence would be: 





BLANK COMMON TABLE 1 


TTEMS() 
e———e TEMS (2) 











This means that the variable A of the first program unit is the same value as the 
array element VALUES(1) in the second program unit. Note that the second program 
unit cannot access K3 of the first program unit because of the difference in blank 
common size. However, KE of the first program unit could correspond to a blank 
common item in still another program unit of the same program. 


Since correspondence is by storage unit: 
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If program unit 1 contains 





and program unit 2 contains 


the following correspondence occurs: 















PROGRAM UNIT 1 PROGRAM UNIT 2 
real part of Al ARAY(1 
imaginary part of Al 


EO 
real part of A2 ARAY(3) 
imaginary part of A2 a ARAY(4) 






Rules: 


(1) A COMMON statement is a specification statement; its order within a program 
unit is shown in Table 2-2. 


(2) An array declarator in a COMMON statement must not contain dummy arguments. 


(3) A block name is not related to any variable or array in the same or any other 
program unit; it has no data type or value associated with it. Once a symbolic 
name is used as a block name, no other name can be used to identify the same 
block throughout the entire program. However, the same symbolic name can be 
used as an array name, variable, or statement function name in the same program 
unit or any other program unit of the program. 


(4) The same block name can occur more than once in a COMMON statement or in 
more than one COMMON statement of a program unit. All entities of the same 
block (blank or labeled) are stored consecutively, in the order of their appearance. 


(5) An EQUIVALENCE statement may extend.a block only past its end, not ahead 
of its beginning. : 


(6) The size, in storage units (see Table 2—4), of a common block for a program unit 
is the sum of the storage required for the elements introduced through COMMON 
and EQUIVALENCE statements. The sizes of labeled common blocks with the same 
block name must be the same in all program units of a program; the sizes of 
blank common in all the program units need not be the same. 
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(7) It is incorrect to equivalence two entities of a COMMON statement to each other, ‘ 
either directly or indirectly. : @ 


(8) Any program unit requiring access to a common block must have a COMMON 
statement. 


(9) In a subprogram, a symbolic name appearing in a COMMON statement may not 
identify an adjustable array. 


(10) Correspondence between different entities of the same common block in different 
program units is determined by order and by storage unit(s). 


(11) In any procedure subprogram, symbolic names of dummy arguments must not appear 
in a COMMON statement. 


(12) An item in blank common must not appear in a DATA statement; an item in 
labeled common may appear in a DATA statement. 


Examples: 


(1) A subroutine subprogram contains: 





ee O,M, M10, ARRA T.A,8,L,£,1 5 TA Bb, E2 4445 TA Bib, E3.( 3 0.)) 


The CALL CLARAY statement in program unit 2 causes all the arrays in the 
labeled common block ARRAYS to have their elements cleared to 0.0. The DO 
range in program unit 2 causes clearing of all elements in the array TABLE1. 
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(@) If program unit 1 contains 


-and program unit 2 will only access MTRX(3,4,5), it must provide space in its 
COMMON statement for those items in blank common that precede MTRX, even 
though they will not be used in program unit 2. For example: 
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9. INITIALIZATION 


GENERAL 


This section describes the DATA statement (data initialization statement) and the 
block data subprogram (specification subprogram). 


DATA STATEMENT 
Function: 


To initialize the values of specified variables and/or array elements at compilation 
time. 


DATA kj /dy/,ko/dy/, --1k,/d,/ 


where: each k is a list of variables and/or array elements separated by commas. 
each d is a list of constants and optionally signed constants separated 
by commas, any of which may be preceded by j*, where j is an integer 
constant indicating repetition of a constant j times. 





Rules: 
(1) Dummy arguments may not appear in a DATA statement. 


(2) Each subscript expression in an array element reference must be an integer 
constant. . 


(3) A Hollerith constant may be used in a list d. 


(4) There must be a one-to-one correspondence between the list-specified items and 
the constants. 


(5S) An initially defined variable or array element must not be in blank common. An 
item in a labeled common block can only be.initialized if the DATA statement 
is in a block data subprogram. 


(6) No item can be initialized more than once in the entire program. 


(7) Array names may not appear in DATA statements. To initialize an entire array, 
each array element must be listed separately. 


(8) A program unit may contain as many DATA statements as required. 
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Example: 





ID | MENiS,! ON, IM .3,) 
DATA, [A.B *,0,.,0 M.( 2 M,(43 Mi(1,)1/,2.*, 0, IT. 71 I ‘ 


The values will be initialized as follows: 


A to 0.0 

B to 0.0 

MATRIX(1) to all blank characters 

MATRIX(2) to 0 

MATRIX(3) to 0 ; - 


9.3. BLOCK DATA SUBPROGRAM 
Function: 


To initialize values of labeled common blocks, at compilation time. 


BLOCK DATA 


data initialization and specification 


statements (except EXTERNAL) 





END 





Rules: 


(1) The first statement of a block data subprogram must be the BLOCK DATA 
statement. 


(2) The body of the block data subprogram consists of one or more DATA statements, 
one or more COMMON statements, and all other required specification statements 
except the EXTERNAL statement. These specification statements are type- 
statements, DIMENSION, and EQUIVALENCE statements. 

(3) All specification statements must precede the DATA statement(s). 


(4) The last line of the subprogram must be an end line. 


(5) The block data subprogram, although independently compiled, is a specification 
subprogram and should not be confused with an external procedure. 


(6) There may be as many block data subprograms as required. 
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(7) If any entity of a common block is initialized in a block data subprogram, all 
items of that block must have their required specification statements, even 
though some of these items do not appear in a DATA statement. 


(8) Initial values may be entered into more than one common block in a single block 
data subprogram. 


(9) A block data subprogram cannot be referenced in a program. 
(10) No executable instructions are generated by the processor for a block data sub- 
program. Therefore, execution time of a program is decreased and less storage 


Space is required for execution. 


Example: 


| tetiacki Dis 

irc. B,2|, ,B,R4 

| love : QIN: . E 

| |i, 0 1,q LR | 

evo MON BLK 1/ sb REZ 7 BIL KP BLD BRAC el ty Ly 
ane A R R (2 2,7,-[T,RU, E,. Bi 2: /j— 1, Ei /1 2). 17,7 812, 8, 1,8,2 18 0,0 / 
BF 

i 


9.3.1. BLOCK DATA Statement 
Function: 


To identify a block data (specification) subprogram. 
BLOCK DATA 


Rule: 


The BLOCK DATA statement must (and can only) appear as the first line of a 
specification subprogram. 
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APPENDIX A. DIFFERENCES BETWEEN 
ANSI FORTRAN AND 
ANSI BASIC FORTRAN 


ANSI (American National Standards Institute, Inc.) Basic FORTRAN lacks some of the 
features found in ANSI FORTRAN. The following is a list of the features which Basic 
FORTRAN lacks, and the sections of this manual affected by the differences: 


(1) There is no DATA statement or specification (block data) Subprogram (2.2.1, 2.2.2, 
Section 9). 


(2) There is no explicit type declaration of any kind; there are no type-statements (2.2.2, 
2.7.2.1, 7.2). 


(3) There is no EXTERNAL statement (2.2.2, 7.1, 8.7). 
(4) There are no ASSIGN or assigned GO TO statements (2.2.2, 4.1, 5.2.3, 5.2.3.1). 
(5) There is no logical assignment statement (2.2.2, 4.1, 4.3). 


(6) All DIMENSION statements must precede all COMMON statements; all COMMON state- 
ments must precede all EQUIVALENCE statements (2.2.2). 


(7) There is no $ character in the Basic FORTRAN character set (2.3). 


@ (8) There is a maximum of five continuation lines instead of 19 (2.4.3). 





(9) There is a maximum of four digits in a statement label instead of five (2.4.4). 


(10) A real constant may not be written as an integer constant followed by a decimal 
exponent (2.6.2). 


(11) Hollerith data is sea ite only in a FORMAT statement; there are no Holter 
constants and no A field descriptor (2.5, 6.3.3, 6.3.3.4). 


(12) There is a maximum of five characters in a symbolic name instead of six (2.7). 
(13) There are no double precision, complex, or logical data types (2.7.2). 


(14) There is no provision for type declaration in a FUNCTION statement (2.7.2.1, 7.2, 
8.5.2.1). 


(15) Arrays are limited to two dimensions instead of three (2.7.4). 


(16) There are no adjustable arrays and no array declaration in COMMON statement (2.7.4.1, 
8.8). 


(17) There are no relational or logical expressions (Section 3). 


(18) There are no double precision or complex type arithmetic expressions (3.2.3). 
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(19) 
(20) 
(21) 
(22) 
(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 
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There is no logical IF statement (5.3). 

There is no provision for extended range in a Do loop (5.4). 

There is a maximum of four instead of five octal digits in a PAUSE statement (5.6.1). 
There is no form control character for formatted output records (6.3). 


There are no D, G, or L field descriptors (6.3.3). 


There is no provision in FORMAT statement for scale factor (6.3.3), data exponent 
on input for F field descriptor (6.3.3.2.2), or second level of parentheses (6.3.3.5, 6.6). 


In numeric fields, blanks are permitted only to the left of the first nonblank character 
and between the sign of the field and the next nonblank character (6.3.3.2). 


In formatted READ and WRITE statements, f must be the statement label of the FORMAT 
statement; no array name is permitted (6.4, 6.5). 


There are no logical statement functions (8.2). 


Basic FORTRAN provides only the following intrinsic functions: ABS, IABS, FLOAT, 
IFIX, SIGN, and ISIGN (8.3). 


Basic FORTRAN provides only the following basic external functions: EXP, ALOG, 
SIN, COS, TANH, SQRT, and ATAN (8.5.1). 


Function subprograms may not define or redefine any of their arguments or alter any 
entity in common or entity associated with common by an EQUIVALENCE statement 


(8.5.2, 8.8). 


There is no provision for labeled common blocks (8.8). 
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Term 


A 


ABS Basic External Function, 
example of 
Actua! Arguments 


Addition, 
operator 
order of evatuation in 
arithmetic expression 


Adjustable Array, 
in OIMENSION statement 
in function Subprogram 
in Subroutine subprogram 


Advance to Next Form, 
form control character for 


AIMAG Intrinsic function, 
example of 
.AND. 


Arguments, 
basic function actual 
function Subprogram actual 
function subprogram dummy 
intrinsic function actual 
statement function actual 
statement function dummy 
subroutine Subprogram actual 
subroutine subprogram dummy 


Reference 


see external 
functions, basic 
8.5.1 


see arguments 


3.2.1 

3.2.1, 3.2.4 
2.7.4.1 

7.3 

8.5.2.2 

8.6.3 

6.3 

see intrinsic 
functions 

3.3.3 

3.4.1 

8.1, 8.1.3, 8.1.4 

8.5.1 

85.2.2, 8.5.2.3 

8.5.2.1, 8.5.2.2 

8.3 

8.2.1, 8.2.2 

8.2.1, 8.2.2 

8.6.1 

8.6.2, 8.6.3 
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&13 


3-5 

3~1, 3-5 
2-25 

7-3 


&-17 
8-25 


6-5 


3-9 


3-10 


8-1, 8-4, 8-4 
8-13 

8-16, 8-23 
8-16, 8-16 
8-9 

8-5, 8-8 
8-5, 8~8 
824 

825, 8-25 






SECTION: 


Index 


INDEX 


Term 


Arithmetic, 
assignment statement 


data types 
expression 
expression, limited 
\F statement 

1/0 data 


operators 
Statement function 


Arithmetic Assignment Statement , 
in sample program 


Arithmetic IF Statement 


Array, 


COMMON statement declaration of 


declaration of 
declaration of adjustable 
DIMENSION statement 

deciaration of 
elements 


* example of adjustable 
function subprogram 
arguments involving 
location of elements in 
uSing name of 


Reference 


see arithmetic 
assignment 
statement 

2.5.1, 3.2.3 

3.2 

8.2.1 

see arithmetic. 
IF statement 

see numeric 
data 

3.2.1 

8.2.1 


5.3.1 


‘274 
8.3 
2.7.4.1 
8.5.2.2 


73 

see array 
element 

8.5.2.2 


8.5.2.2 
2.7.4.3 
2.7.4.1 





PAGE: 


Page 


211, 3— 
3-1 
&5 


2-23 
8-32 
2-24 
8-17 


7-3 


8-20 


8~17 
2-29 
2~25 
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Term 


Array Element. 
EQUIVALENCE statement for 
function Subprogram arguments 
involving 
location of 


number of subscript expressions 


required by 
Subscript expression in 
reference to 


Assembly Language 
ASSIGN Statement 


Assignment Statements, 
arithmetic 
GO TO 
logical 


Assigned GO TO Statement 


BACKSPACE Statement 
Basic External Functions 


Basic Field Descriptor in 
FORMAT Statement 


Bit. Definition of 


Blank Characters. 
in comment lines 
imend line 
field descriptor 
in FORMAT statement 
in FORTRAN coded lines 
in Hollerith data 
in! O numeric data 
in statement labels 
in symbolic names 


Blank Common 


BLOCK DATA statement 


Biock Data Subprogram. 
order of statements in 


Block name 


Reference 


2.7.4.2 
74 


8.5.2.2 
2.7.4.3 


7A 


2.7.4.2 


1.2.2 


5.2.3.1 


4.2 
5.2.3.1 
43 


5.2.3 


6.3.3 


1.2.1 


1.6, 2.4.1 
2.4.2 
6.3.3.1 
6.3 

2.4 

2.4.3 
6.3.3.2 
2.4.4 

206 


see COMMON 
statement 


9.3.1 


2.2.1, 9.3 
Table 2-2 


see COMMON 
statement 
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2-26 
7-5 


8-17 
2-29 


2-26 


1-3 


5~8 


4-1 
5-8 
4-4 


6-30 


2-3 


1-12, <-7 
2~8 

6-8 

6-5 

2~7 

2-9 

6-8 

2-10 
2-18 


2-2. 3-2 
2-6 








Term 


Cc 


CALL statement, 
within DO range 


Central Processing Unit, 
representation of fixed-point 
numbers 
representation of floating- 
point nunbers 


Characters, 
alphanumeric 
form control 
FORTRAN set of 
prccessor set of 
specia! 


Closed Internat Block, 
example of 


Comment Line, 
sample 


COMMON Statement 
Compilation Process 
Compiler 
Completely Nested Nest 
Compiex Type. 

constant 

field descriptor 

memory requirements 
Computed GO TO Statement 
Computer. 

decimal 

hexadecimal 

octal 
Constant(s). 

basic real 

integer 


real 


Continuation Line(s) of Statement 


Index 
SECTION: 


Reference 


8.6.1 
see DO 
statement 


15.1.3 
1.5.1.3.1 


1.5.1.3.2 


23 
Table 6-1 
2.3 
2.3 
2.3 


2.4.1 
1.6 


3.8 
1.4 
130 


see DO 
statement 


2.5.1, 2.5.1.4 
2.6.4 
6.3.3.2.3 
Table 2—4 


5.2.2 


1.2.1 
1.2.1 
1.2.1 


2.6 

2.6.2 
2.6.1 
2.6.2 














8-24 


1-8 


1-8 


2-7 
65 
2-7 
2-7 
2-7 


2~7 
1-9 


8-32 


14 


1-4 


2-11, 2-13 
2-16 
6-12 
2—11 


5-2 


1-2 
1-2 
1-2 


2-13 
2-14 
2-14 
2-14 
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Term 


CONTINUE Statement, 
as terminal statement 


Control Statements 


Control Variable, 
in DO statement 
in 00-implied list 


Conversion Codes in FORMAT 
Statement 


CPU 


Cube Root, 
function subprogram for 
catculation of 


DATA Statement, 
in sample program 


Data Types, 
Conversion of 
declaration of Hollerith 
explicit declaration of 
implied deciaration of 
in arithmetic expressions 
in relational expressions 
memory requirements for 
of function name 
of symbolic names 


Decimai Computer 
Decimal Exponent 
Deciaration, 
array 
DIMENSION statement array 
explicit type 
Hollerith value 
implied type 
DIMENSION Statement 


OMOD Basic External Function, 


example using 





Reference 


5.5 
see DO 
statement 


5.1 


5.4 
6.2.2 


6.3.3 


see central 
processing unit 


8.5.2.2. 


9.2 
1.6 


2.5.1 
42 
2.7.2.3 
2.7.2.1 
2.7.2.2 
3.2.3 
3.3.2 
Table 2—4 
8.5.2.1. 
2.7.2 


1.2.1 


2.6.2 


2.7.4.1 

7.3 

2.7.2.1, 7.2 
2.7.2.3 
2.7.2.2 


7.3 
see external 


functions, basic 
4.2 
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5-21 


5-13 
6-3 


6-7 


8-21 


9-1 
1~12 


2-11 
4-2 

2-22 
2-20 
2-20 
3—4 

3-8 

2-11 
8-16 
2-20 


1-2 

2-15 
2-24 

7-3 

2-20, 7-1 
2-22 


2-20 


7-3 


4~3 


Term 


00-implied, | 
list 
specification 


DO Statement 
Doubie Precision Exponent 


Double Precision Type, 
constant 
field descriptor for 1/0 data 
memory requirements 


Dummy Arguments 


ENDFILE Statement 


End Line, 
in function subprogram definition 
in sample program 
in subroutine subprogram 
definition 


EQUIVALENCE Statement - 
Executable Statements 


Execution, 
hait, temporary 


process 

sequence 

sequence modified by control 
statements 

termination 


Explicit Dectaration of 
data type 


Exponentiation, 
type rules for 


Expressions, 
arithmetic 
limited arithmetic 
limited fogical 
logical 
relational 
subscript 


SECTION: 


Index 


Reference 


6.2.2 
6.2.2 


5.4 
2.6.3 


2.5.1, 2.5.1.3 
2.6.3 
6.3.3.2.4 
Tabie 2-4 


see arguments 


6.8.3 


2.4.2 
8.5.2.2 
1.6 


8.6.3 
7.4 


2.2.2 


see PAUSE 
statement 

14 

2.2.1 


5.1 
see STOP 
statement 


2.7.2.1 


3.2.3 


3.1 
3.2 
8.2.1 
8.2.2 
3.4 
3.3 
2.7.4.2 
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6-2 
5-13 
2-16 
2-11, 2- 
2-16 


6-14 
2-11 


6-32 
2-8 

8-16 
1-13 
8-25 


7-4 


2-5 


1-4 
2-1 


2-20 


3-1 
3-1 
8-5 


3-10 
3~-7 
2-26 
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’ Term 


Extended range 


External Function(s), 
basic 
subprograms 


External Procedure(s), 
name as an argument 


subprograms, out-of-line 
machine coding of 


EXTERNAL Statement 


FE 


Factor, 
in arithmetic expression 
in logical expression 
. FALSE. 
field descriptor(s) in FORMAT 
Statement, 
blank 
repetition of 
Fietd Separators) 
Field Width ii Field Descriptors 


File 


Fixed-Point Representation of 
Numbers 


FLOAT Intrinsic Function, 
example of 


Floating-Point Representation of 
Numbers 


Flowchart. sample 
Form Contro! Characters 


Form, Typical FORTRAN 
programming 


FUNDAMENTALS OF FORTRAN 


Reference 


see 00 
statement 


8.5 

8.5.1 

8.5.2 

2.2.1 

see EXTERNAL 
statement 


8.1.2 


8.7 


3.2.2 
3.4.2 


6.3.3 
6.3.3.1 
6.3.3.5 


6.3.3 


6.1 


1.5.1.3.1, 2.5.1.1 


see intrinsic 


functions 


3.2.4, 4.2 


1.5.1.3.2, 2.5.1.2 


1.6 


Table 6-1 


2.4 
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813 
8-13 
8—15 


2-2 


8-3 


8-30 


3-2 
3-11 


2-17 


6-7 
6-8 
6-21 


6-7 


6-1 


1-8, 2-12 


3-6, 4-3 


1-9, 2-12 


1-10 


6-5 













Index 


SECTION: 





Term 


Format Control and 1/0 list, 
relation between 


FORMAT Statement, 
in sample program 


Function Reference, 
as primary in arithmetic 
expression 
as primary in logical 
expression 


FUNCTION Statement 


Function Subprogram(s) 
definition of 
order of statements in 
references to 


GO TO Statements, 
assigned 
computed 
unconditional 


i 


Hexadecimai Computer 


Hollerith, 

constants 

constant in CALL statement, 
example of 

constants in DATA statement, 
example of 

data 

in DATA statement 

in input, output statement 

field descriptors in FORMAT 
statement 

integer or logical type 

representation by real 


IF Statement(s) 
arithmetic 


logical 


Reference 


6.6 


6.3 
1.6 


8.1 
3.2.2 
3.4.2 
8.5.2.1 


“2.2.1, 8.5.2 
8.5.2.1 
Table 2-2 
8.5.2.3 


5.2 

5.2.3 
5.2.2 
5.2.1 


1.2.1 


2.6.6 | 
3.6.3 


6.3 
2.5 
9.2 
6.3.3 


6.3.3.4 
2.7.2 
2.5.1, 2.5.1.6 


5.3 

see arithmetic IF 
statement 

see logical IF 
statement 
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6—26 


111 
8-1 
3-2 
3-11 
8-16 
2-2, 8-15 
8~16 


2-6 
8-23 


1-2 
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6-18 
2-20 
2-11, 2~-13 














UP-7536 
Rev. 1 









Term 


Imptied Type Declaration 
tncrementation Parameter, 
in DO-implied tist 
in DO statement 


Initialization Statement 


tnitial Line of Statement 
Initial Parameter, 
in DO-implied list 


in DO statement 


Input Device(s), 
togical unit number of 


Input List 


Input Statement(s), 
auxiliary 


Integer Type, 
constants 
conversion of !/O data 
memory requirements 


Intrinsic Functions, 
inline machine coding of 


1/0 List 
K 


“Keyword, 
use as symbolic name 


L 


Label, Statement 


Labeled Common Block, 


initialization of 


Language(s), Programming, 
assembly 
FORTRAN 
machine 


Reference 


2.7.2.2 


see DO statement 


6.2.2 
5.4 


see DATA 
statement 


2.4.3 


6.2.2 
5.4 


15.11 
6.2.1 


6.2.2 


6.1 
6.8 


2.5.1, 2.5.1.1 
2.6.1 

6.3 

Table 2—4 


8.3 
8.1.1 


6.2.2 


1.6 
2.7 


see statement 
label 


see COMMON 
statement 

see block data 
subprogram 


1.2 

1.2.2 
1.2.3 
1.2.1 


FUNDAMENTALS OF FORTRAN 


Page 


2—20 


6-2 


6-2 


 §-13 


1-8 
6-2 


6-2 


6-1 
6-29 


2-11, 2-12 
2-14 


2-11 


89 
8-3 


6-2 


1-13 
2-18 


1-2 
1-3 


1-4 
1-2 









Term 


Logical Assignment Statement 


Logical Expression, 


limited 


Logical IF Statement, 
in sample program 


Logical Operators 
Logical Statement Function 
Logical Type, 
constant 
field descriptor 
memory requirements 
Logical Unit Number in 1.0 Statement 
Loop, 


example of nested 
in sample program 


Mi 


Main Program, 
order of statements in 


N 


Name(s), 
symbolic 


Newton-Raphson Method, 

in computation of cube root 
Nonexecutable Statements 
.NOT. 


Numeric Data, 
‘in tO data 


Oo 


Octal Computer 


Operating System 


SECTION: 


Index 


Reference 


4.3 


3.4 
8.2.2 


5.3.2 
1.6 


3.4.1 

6.2.2 

2.5.1, 2.5.1.5 
2.6.5 

6.3.3.3 
Table 2—4 
6.2.1 


2.7.4.2 
1.6 


2.2.1 
Table 2-2 


_ see symbolic 


name 


8.5.2.2 


2.2.2 


3.4.1 


6.3.3.2 


1.2.1 


15, 1.5.2 
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3-10 
8-8 


5-11 
1-11 


3-10 
8-8 
2-11, 2 
2-17 
6-16 
2-11 
6-2 


2-28 
1-11 


21 
2-6 


&21 


2-5 


3-10 


6-8 


1-2 


1-7, 1-9 
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Term 


Operators, 
arithmetic 


OR. 


Qutput Device(s), 


logical unit number of 


Output List 


Output Statement(s), 
auxiliary 


P 


Parameter(s), 
of DG-implied list 
of DO statement 


PAUSE Statement 


Polynomial Evaluation, 


by DO statement 


Primary, 


in arithmetic expression 
in logical expression 


Procedure(s), 
external 


Procedure Subprograms, 
order of statements in 


Processor, Definition of 


Program, 
compilation of 
execution of 
FORTRAN source 
main 
obyect 
sample executable 
source 
units of 


Program Unit(s), 


communication between 


definition of 
execution of 
organization of 


FUNDAMENTALS OF FORTRAN 


Reference 


3.2.1 


3.4.1 


15.1.5 
6.2.1 


6.2.2 


6.1 
6.8 


6.2.2 
5.4 


§.6.1 


5.4 


3.2.2 
3.4.2 


8.1 
2.2.1 


2.2.1 
Table 2-2 


1.2 


14 

see execution 

13 

2.2.1 

13 

1.6 

1.2 

see program units 


8.1.3 
2.2.1 
2.2.1 
2.2.2 


Page 


3-1 


3-10 


1-9 
6-2 


6—2 


6-2 
5-13 


5-22 


5-21 


3-2 
3-11 


2-2 


2-1 
2-6 


1-2 


1-4 


1-4 
2-1 
1-4 
1-9 
1-2 


2-1 
2-3 
2-5 








Term 


Programming Form, 
typical FORTRAN 


Punched Cards 
for sampie program 


Quadratic Equation, 
solution of 


READ Statement, 
formatted 
sample program 
unformatted 


REAL Intrinsic Function, 
example of 


Real Type, 
constant 
data input 
data output 
memory requirements 


Record(s), 
demarcator 


Relational , 
expression 


operators 


Repeat Count, 
in FORMAT statement 


RETURN statement. 
multiple use of 


REWIND Statement 


Scate Factor, 
with field descriptors 


Sequential File 


Index 


SECTION: 


Reference 


2.4 


1.6 


5.3.1 


6.4 
1.6 
6.7 


see intrinsic 
functions 
3.3.3 


2.5.1, 2.5.1.2 
2.6.2 
6.3.3.2.2 
6.3.3.2.3 
Table 2-4 
6.1 

6.3.1 


3.3 


3.3.1 


6.3.3, 6.3.3.5 


8.4 
8.5.2.2 


6.8.1 


6.3.3. 6.3.3.6 
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1-9 
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3-9 
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2-14 
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61 
6-6 





6~7, 6-21 


8~13 
8-16 


$—29 


6-7. 6-Z2 











UP-7536 
Rev. 1 









Term 
Simple List in 1/0 List 
Software, Definition of 


Sort, Internal 
example 


Specification Statements 


Specification Subprogram 


Statement Label(s), 
in sample program 


Statement Functions, 
arithmetic 
inline machine coding of 
logical 


Statements, 
continuation tine of 
executable 


initial line of 
fabels of 


list of FORTRAN 
nonexecutabie 


order of 
terminal 


STOP Statement(s), 
-in sample program 


Storage, 
auxiliary 
EQUIVALENCE statement 
conservation of main 
main 


Subprograms, 
description of 
block data 


function 


subroutine 


SUBROUTINE Statement 


FUNDAMENTALS OF FORTRAN 


Reference 
6.2.2 


1.5.2 


5.4 
71 


see block data 
subprogram 


2.4.4 
1.6 


8.2 

8.2.1 
8.1.1 
8.2.2 


2.4.3 

2.4.3 

see executable 
statements 

2.4.3 

see statement 
labels 

Table 2-1 

see nonexecutabie 
statements 

Table 2-2 

see DO statement 


5.6.2. 
1.6 


1.5.1.4 


74 
1.5.1.2 


2.2.1 

see block data 
subprograms 

see function 
subprograms 

see subroutine 
subprograms 


8.6.2 


Page 


6-2 


5-20 


T~1 


2-10 
1-9 


8-2 
8-5 
8-2 
8-8 
2-9 
2-9 


2-9 


2-6 


5-23 
1-9 


8-25 


Term 


Subroutine Subprogram(s), 
definition of 
invoking of 
order of statements in 


Subscript; 
expression 


Symbolic Names, 
uniqueness of 


T 


Term, 
in arithmetic expression 
in logical expression 
Terminal Parameter 
Terminal Statement 
Termination, 


of compilation 
of execution 


TRUE. 


Types of Data 


Type-Statement(s), 
examples of 


UW 


Unconditional GO TO Statement 


V 


Variable(s), 
control 


defining values of 
‘type declaration of 


Ww 


WRITE Statement, 
formatted 
unformatted 


SECTION: 


Index 


Reference Page 
2.2.1, 8.6 2-2, &—2¢ 
8.6.3 8-25 
8.6.1 8-24 
Table 2-2 2-6 
2.7.4.1 2-24 
2.7.4.2 2-26 
2.7 2-18 
2.7.1 2-19 
3.2.2 3-2 
3.4.2 3-11 


see DO statement 


see DO statement 


see end fine 
see STOP statement 


2.6.5 2-17 


see data types 


7.2 7-1 
2.7.2.1 2-20 

5.2.1 5-2 
2.7.3 2-23 
see DO statement, 
DO-implied list 

2.7.3 2-23 
2.7.2.2, 2.7.2.3 2-20, 2- 

6.5 6-25 

6.7 6-28 
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